信号处理Matlab试验toSt.doc

例1： %两信号叠加 clear all; f1=10; f2=30; t=(0:0.001:1); y1=sin(2*pi*f1*t); y2=3*cos(2*pi*f2*t); y=y1+y2; figure(1) subplot(5,1,1),plot(t,y1); axis([min(t),max(t),min(y1),max(y1)]); subplot(5,1,2),plot(t,y2); axis([min(t),max(t),min(y2),max(y2)]); subplot(5,1,3),plot(t,y); axis([min(t),max(t),min(y),max(y)]); %--------------------------------------------- %加噪声干扰 z=randn(size(t)); yr=y+z; subplot(5,1,4),plot(t,z); axis([min(t),max(t),min(z),max(z)]); subplot(5,1,5),plot(t,yr); axis([min(t),max(t),min(yr),max(yr)]); 试验内容： 改变采样频率，观察y1、y2波形的变化； 分别对y1、y2和y加噪声，观察波形变化。 例:2： %方波 clear all fs=1000; dlt=1/fs; len=2048; t=(0:len-1)*dlt; f1=5; f2=10; f3=20; s1=sin(2*pi*f1*t); s2=sin(2*pi*3*f1*t); s3=sin(2*pi*5*f1*t); s4=sin(2*pi*7*f1*t); s5=sin(2*pi*9*f1*t); s6=sin(2*pi*11*f1*t) s7=sin(2*pi*13*f1*t) s=(s1+s2/3+s3/5+s4/7+s5/9+s6/11+s7/13)*4/pi; figure(1) subplot(6,1,1),plot(t,s1); axis([min(t),max(t),min(s1),max(s1)]) subplot(6,1,2),plot(t,s2); axis([min(t),max(t),min(s2),max(s2)]) subplot(6,1,3),plot(t,s3); axis([min(t),max(t),min(s3),max(s3)]) subplot(6,1,4),plot(t,s); axis([min(t),max(t),min(s),max(s)]) subplot(6,1,6) y=do_fft(s); dlf=fs/len; ff=(0:len/2-1)*dlf; plot(ff,y); axis([0,100,min(y),max(y)]) 试验内容： 改变各信号（s1-s7）的振幅和初相位，观察波形变化； 改变信号（s1-s7）的频率，观察傅里叶变换得到的谱线变化； 例3： %窗函数和一个随机信号相乘（窗函数采样） clear all %t=linspace(0,100); fs=1000; dlt=1/fs; len=100; t=(0:len-1)*dlt; x=randn(size(t)); w1=gausswin(100,6); %高斯窗函数 w2=hanning(len); %汉宁窗 w3=hamming(100); %海明窗 w4=kaiser(100,6) s=w2; y1=x.*s; figure(1) subplot(3,1,1),plot(t,x); axis([min(t),max(t),min(x),max(x)]) subplot(3,1,2),plot(t,w1); axis([min(t),max(t),min(w2),max(w2)]) subplot(3,1,3),plot(t,y1); axis([min(t),max(t),min(y1),max(y1)]); 试验内容： 改变窗函数类型，观察波形变化； 改变窗函数长度，观察波形变化； 例4： %离散傅立叶变换和其逆变换 clear all f1=50; f2=120; t=(0:0.001:1); y1=0.5*sin(2*pi*f1*t); y2=2*sin(2*pi*f2*t); y=y1+y2; plot(y(1:50)); s=fft(y,512); Pyy=s.*conj(s)/512; f=1000*(0:155)/512; figure,plot(f,Pyy(1:156)) %---------