【精选资料】matlab信号分析

一、用MatLab完成三种不同信号的产生、合成或分解，时域波形分析(峰值，峰峰值，有效值，平均t=linspace(0,5*pi,65000);%在从0到2*pi取15000个点y1=20*sin(2*pi*t);%任给一正弦信号subplot(3,1,1);plot(t,y1);axis([02*pi-2525]);%画时域波形图max1=max(y1)%峰值vpp1=max(y1)-min(y1)%峰峰值avg1=mean(y1)%平均值rms1=norm(y1)/sqrt(length(y1))%均方根值（有效值）title('正弦波')y2=20*square(2*pi*t,50);%任给一方波信号subplot(3,1,2);plot(t,y2);axis([02*pi-2525]);%画时域波形图max2=max(y2)%峰值vpp2=max(y2)-min(y2)%峰峰值avg2=mean(y2)%平均值rms2=norm(y2)/sqrt(length(y2))%均方根值（有效值）title('方波')y3=20*sawtooth(t,0.5);%任给一三角波信号subplot(3,1,3);plot(t,y3);axis([05*pi-2525]);%画时域波形图max3=max(y3)%峰值vpp3=max(y3)-min(y3)%峰峰值avg3=mean(y3)%平均值rms3=norm(y3)/sqrt(length(y3));%均方根值（有效值）title('三角波')figuref1=4*cos(t)/(pi*pi);f2=4*cos(3*t)/(pi*pi*9);f3=4*cos(5*t)/(pi*pi*25);f4=4*cos(7*t)/(pi*pi*49);f=0.5*(-sawtooth(t,0.5)+1);%三角波的分解subplot(5,1,1),plot(t,f);title('三角波')%画三角波subplot(5,1,2),plot(t,f1);title('三角波的基波')%画基波subplot(5,1,3),plot(t,f2);title('三角波三次谐波')%画三次谐波subplot(5,1,4),plot(t,f3);title('三角波五次谐波')%画五次谐波subplot(5,1,5),plot(t,f4);title('三角波七次谐波')%画七次谐波legend('三角波的分解')figuref=0.5*(-sawtooth(t,0.5)+1);%三角波的合成f1=0.5+4*cos(t)/(pi*pi);f2=0.5+4*cos(t)/(pi*pi)+4*cos(3*t)/(pi*pi*9);f3=0.5+4*cos(t)/(pi*pi)+4*cos(3*t)/(pi*pi*9)+4*cos(5*t)/(pi*pi*25);f4=0.5+4*cos(t)/(pi*pi)+4*cos(3*t)/(pi*pi*9)+4*cos(5*t)/(pi*pi*25)+4*cos(7*t)/(pi*pi*49);subplot(2,2,1),plot(t,f1);title('三角波基波')subplot(2,2,2),plot(t,f2);title('三角波基波+三次谐波')subplot(2,2,3),plot(t,f3);title('三角波+三次谐波+五次谐波')subplot(2,2,4),plot(t,f4);title('三角波+三次谐波+五次谐波+七次谐波')legend('三角波的合成')值，均方根值等）二、用MatLab进行三种不同信号频谱分析（可选择功率谱，幅频相频谱，实频虚频）和相关分析（可选择自相关，互相关）Fs=5120;%设定采样频率dt=1.0/Fs;N=1024;%采样点数T=dt*N;%设置常量t=linspace(0,T,N);%在从0到T取N个点y=50*sin(2*pi*50*t);%任给一正弦信号subplot(3,3,1);plot(t,y);title('50sin(2*pi*50*t)的时域波形图');%画时域波形图subplot(3,3,2);yy=fft(y,N);%对信号进行快速傅里叶变换P=2*abs(yy)/Nf=linspace(0,Fs,N);%取幅值的频率plot(f,P);xlabel('频率(Hz)');ylabel('幅值');title('50sin(2*pi*50*t)的幅值谱');%画幅值谱axis([0500055])subplot(3,3,3);anglex=angle(yy)*180/pi;%求幅角plot(f,anglex);xlabel('频率(Hz)');ylabel('相位（度）')title('50sin(2*pi*50*t)的相位谱');%画相位谱axis([0100-200200])subplot(3,3,4);plot(f,real(yy));xlabel('频率(Hz)');ylabel('实部');%画实部谱title('50sin(2*pi*50*t)的实部谱');axis([0500-500500])subplot(3,3,5);plot(f,imag(yy));xlabel('频率(Hz)');ylabel('虚部');%画虚部谱title('50sin(2*pi*50*t)的虚部谱');axis([0500-300005000])subplot(3,3,6);plot(f,P.^2);xlabel('频率(Hz)');ylabel('功率');%画功率谱title('50sin(2*pi*50*t)的功率谱');axis([0500-2003000])y1=50*cos(2*pi*50*t);%任给一同频率余弦信号subplot(3,3,7);plot(t,y1);%画时域波形图title('50*cos(2*pi*50*t)的时域波形图');subplot(3,3,8);s=xcorr(y);%正弦信号的自相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s);xlabel('频率(Hz)');%画自相关函数图像title('正弦信号的自相关函数');subplot(3,3,9);s1=xcorr(y,y1);%正余弦信号的互相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s1);xlabel('频率(Hz)');%画互相关函数图像title('正余弦信号的互相关函数');figurey2=50*sinc(2*pi*50*t);%任给一sinc函数信号subplot(4,2,1);plot(t,y2);%画时域波形图title('50sinc(2*pi*50*t)的时域波形图');axis([00.18-2050])subplot(4,2,2);Pyy=fft(y2,N);%对信号进行快速傅里叶变换P=abs(Pyy)/Nf=linspace(0,Fs,N);%取复数形式的幅值的频率plot(f,P);%画幅值谱xlabel('频率(Hz)');ylabel('幅值');title('50sinc(2*pi*50*t)的幅值谱');axis([050001])subplot(4,2,3);anglex=angle(Pyy)*180/pi;%求幅角plot(f,anglex);%画相位谱xlabel('频率(Hz)');ylabel('相位（度）')title('50sinc(2*pi*50*t)的相位谱');axis([01000-10010])subplot(4,2,4);plot(f,real(Pyy));xlabel('频率(Hz)');ylabel('实部');%画实部谱title('50sinc(2*pi*50*t)的实部谱');axis([0500-100550])subplot(4,2,5);plot(f,imag(Pyy));%画虚部谱xlabel('频率(Hz)');ylabel('虚部');title('50sin(2*pi*50*t)的虚部谱');axis([0500-700100])subplot(4,2,6);plot(f,P.^2);%画功率谱xlabel('频率(Hz)');ylabel('功率');title('50sinc(2*pi*50*t)的功率谱');axis([050001])subplot(4,2,7);s=xcorr(y2);%sinc信号的自相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s);%画自相关函数图像xlabel('频率(Hz)');title('sinc信号的自相关函数');subplot(4,2,8);s1=xcorr(y,y1);%sinc信号与余弦信号的自相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s1);%画互相关函数图像xlabel('频率(Hz)');title('sinc信号与余弦信号的自相关函数');figurey3=20*square(2*pi*50*t,50);%任给一方波信号subplot(4,2,1);plot(t,y3);%画时域波形图title('20*square(2*pi*50*t,50)的时域波形图');axis([00.1-2525])subplot(4,2,2);Pyy=fft(y3,N);%对信号进行快速傅里叶变换P=abs(Pyy)/N;f=linspace(0,Fs,N);%取幅值的频率plot(f,P);%画幅值谱xlabel('频率(Hz)');ylabel('幅值');title('20*square(2*pi*50*t,50)的幅值谱')axis([02000-315])subplot(4,2,3);anglex=angle(Pyy)*180/3.14;%求幅角plot(f,anglex);%画相位谱xlabel('频率(Hz)');ylabel('相位（度）')title('20*square(2*pi*50*t,50)的相位谱');axis([02000-200200])subplot(4,2,4);plot(f,real(Pyy));%画实部谱xlabel('频率(Hz)');ylabel('实部');title('20*square(2*pi*50*t,50)的实部谱');axis([02000-400500])subplot(4,2,5);plot(f,imag(Pyy));%画虚部谱xlabel('频率(Hz)');ylabel('虚部');title('20*square(2*pi*50*t,50)的虚部谱');axis([02000-150001000])subplot(4,2,6);plot(f,P.^2);%画功率谱xlabel('频率(Hz)');ylabel('功率');title('20*square(2*pi*50*t,50)的功率谱');axis([02000-30200])subplot(4,2,7);s=xcorr(y3);%方波信号的自相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s);%画自相关函数图像xlabel('频率(Hz)');title('方波信号的自相关函数');subplot(4,2,8);s1=xcorr(y,y3);%方波与正弦信号的互相关函数x1=linspace(0,2*T,2*N-1);plot(x1,s1);xlabel('频率(Hz)');title('方波与正弦信号的互相关函数');%画互相关函数图像-高氯酸对阿胶进行湿法消化后,用导数火焰原子吸收光谱技术测定阿胶中的铜、“中药三大宝,人参、鹿茸和阿胶。”阿胶的药用已有两千多年的悠久历史,历代宫①马作峰.论疲劳源于肝脏[J].广西中医药,2008,31(1):31.①史丽萍,马东明,解丽