随机信号大作业

1、 随机信号分析大作业 班级：130204xx 姓名：xxx 学号：13xxxxxxxxx 1、 设有随机初相信号X(t)=5cos（t+），其中相位是在区间（0,2）上均匀分布的随机变量。试用Matlab编程产生其三个样本函数。 程序如下：clc;clear all;y=2*pi*rand(1,3);t=0:pi/20:8*pi;x1=5*cos(t+y(1);figure(1)plot(t,x1,'r-');hold on;x2=5*cos(t+y(2);plot(t,x2,'b');x3=5*cos(t+y(3);plot(t,x3,'g')

2、;xlabel('t');ylabel('X(t)');grid on;%axis tight;title('随机初相信号的三个样本函数');实现结果2.利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。(1)分析复合信号的功率谱密度，幅度分布特性；程序如下clc;clear all;fs=100; %采样信号频率fc=2; %正弦信号频率t=0:1/fs:2;n=2*fs+1; %信噪比x=sin(2*pi*fc*t);y=awgn(x,2);FY=fft(y); %求复合信号的傅里叶变换FY1=fftshift(FY);f=(0:n

3、-1)*fs/n-fs/2;plot(f,abs(FY1),'r-');title('复合信号的幅频分布特性');xlabel('w');ylabel('F(jw)');grid on;C,a=xcorr(y,'unbiased'); %C为复合信号的自相关函数FY2=fft(C); % FY2是对自相关函数求傅里叶变换，即为功率谱密度FY3=fftshift(FY2);figure(2);plot(a*0.25,abs(FY3),'b-');title('复合信号的功率谱密度');

4、xlabel('w');ylabel('G(w)');grid on;实现结果(2)分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性；程序如下clc;clear all;fc=2;fs=100;R=2;C=0.5;t=0:1/fs:2;n=2*fs+1;T=R*C;UC0=4;w=2;Zc=1/j/w/C;x=sin(2*pi*fc*t);y=awgn(x,2);absH=abs(Zc/(R+Zc);phiH=angle(Zc/(R+Zc);ucst=absH*sin(2*pi*fc*t+phiH)+y-x; %复合信号通过RC电路后的信号FY=f

5、ft(ucst); %求复合信号通过RC电路后的傅里叶变换，即求幅频特性FY1=fftshift(FY); %使正负半轴图像对对称f=(0:n-1)*fs/n-fs/2;plot(f,abs(FY1),'r-');title('通过RC电路后信号的幅频特性曲线');xlabel('w');ylabel('F(jw)');grid on; C1,a=xcorr(ucst,'unbiased'); %求复合信号通过RC电路后的自相关函数FY2=fft(C1); %对自相关函数求傅里叶变换，即求功率谱密度 FY3=fft

6、shift(FY2); figure(2);plot(a*0.25,abs(FY3),'b-');title('复合信号通过RC电路后的功率谱密度');xlabel('w');ylabel('G(w)');grid on;实现结果（3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。程序如下clc;clear all;fc=2;fs=100;t=0:1/fs:2;n=2*fs+1;x=sin(2*pi*fc*t);y=awgn(x,2); %求得的复合信号b=fir1(32,0.24);y0=filter(b,1,

7、y); %对复合信号进行一维低通滤波；FY=fft(y0); %对滤波后的信号进行傅里叶变换，得频域函数FY1=fftshift(FY);f=(0:n-1)*fs/n-fs/2;plot(f,abs(FY1),'r-');grid on;xlabel('w');ylabel('F(jw)');title('复合信号通过低通滤波器后的幅频分布特性');C2,a=xcorr(y0,'unbiased'); %复合信号通过低通滤波器后的自相关函数FY2=fft(C2); %对自相关函数求傅里叶变换，得到功率谱密度函数FY3=fftshift(FY2);figure(2);plot(a*0.25,ab