2.31上机题:设有随机初相信号X(t)=5cos(t+φ) ,其相位φ是在区间[0,2pi]上均匀分布的随机变量。试用MATLAB编程产生其三个样本函数。
代码如下:
clc,clear
close all
rng('default')%随机种子
t=0:0.001:10;
phi=rand(1,3)*2*pi;
x(1,:)=5*cos(t+phi(1));
x(2,:)=5*cos(t+phi(2));
x(3,:)=5*cos(t+phi(3));
figure()
hold on;
plot(t,x(1,:),'linewidth',2);
plot(t,x(2,:),'linewidth',2);
plot(t,x(3,:),'linewidth',2);
xlabel('时间');
ylabel('幅度');
grid on;
set(gca,'FontWeight','bold','FontSize',10);
运行结果图:
3.24上机题:利用matlab程序设计一正弦型信号加高斯白噪声的复合信号。
(1)分析复合信号的功率谱密度、幅度分布特性;
(2)分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性;
(3)分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。
代码如下:
clc;clear
close all
%%产生一个正弦信号
fs=100;%采样频率为100HZ
fc=10;%信号频率为10HZ
t=0:1/fs:2;
x=sin(2*pi*fc*t);
y=1*randn(size(t));
s=x+y;
[a,b]=xcorr(s,'coeff');
subplot(5,2,1);
plot(t,x,'r')
title('x=sin(2*pi*fc*t)');
ylabel('x')
xlabel('t')
grid;
subplot(5,2,2);
plot(t,y,'r');
title('y=1*randn(size(t))');
ylabel('y');
xlabel('t');
grid;
subplot(5,2,3);
plot(t,s);
title('复合函数s');
ylabel('s');
xlabel('t');
grid;
subplot(5,2,4);
plot(b*(1/fs),a);
title('复合函数s的自相关函数');
ylabel('b');
xlabel('a');
grid;
Fs=fft(s);
Fs1=fftshift(Fs);
f=(0:200)*fs/201-fs/2;
subplot(5,2,5);
plot(f,abs(Fs1),'r');
title('复合函数s的频谱图');
ylabel('F(jw)');
xlabel('w');
grid;
p=Fs1.*conj(Fs1)/length(Fs1);
subplot(5,2,6);
plot(f,p,'r');
title('s的功率谱密度图');
ylabel('G(w)');
xlabel('w');
grid;
b=10;
s2=conv2(s,b*pi^-b*t);
Fs2=fftshift(fft(s2));
f=(0:400)*fs/201-fs/2;
subplot(5,2,7);
plot(f,abs(Fs2),'r');
title('s通过RC后函数频谱图');
ylabel('Fs2(jw)');
xlabel('w');
grid;
p1=Fs2.*conj(Fs2)/length(Fs2);
subplot(5,2,8);
plot(f,p1,'r');
title('s通过RC后的功率谱密度图');
ylabel('Gs2(w)');
xlabel('w');
grid;
s3=conv2(s,sin(10*t)/(pi*t));
Fs3=fftshift(fft(s3));
f=(0:200)*fs/201-fs/2;
subplot(5,2,9);
plot(f,abs(Fs3),'r');
title('s通过理想滤波器后函数频谱图');
ylabel('Fs3(jw)');
xlabel('w');
grid;
p2=Fs3.*conj(Fs3)/length(Fs3);
subplot(5,2,10);
plot(f,p2,'r');
title('s通过理想滤波器后的功率谱密度图');
ylabel('Gs3(w)');
xlabel('w');
grid;
运行结果图: