实验2matlab基础及信号频谱分析.doc
实验2matlab基础及信号频谱分析实验日期:一.实验目的1.熟练掌握MATLAB基础知识。2.掌握应用MATLAB软件进行通信系统仿真的方法。3.掌握应用MATLAB求解音频信号的时域波形和幅频特性曲线。4.分析频谱特性曲线。二.仿真框图本实验主要工作是掌握MATLAB基础知识并应用MATLAB对通信系统进行仿真了解频谱特性曲线,本实验的仿真框图如图1:图1:本实验仿真框图实验任务1:框(1)、(2)和(3)组成,对音频“F215.wav”分别绘制时域波形、幅频特性曲线实验任务2:框(4)、(5)和(6)组成,对连续信号用fs=100Hz频率对其采样,分别绘制N=128、1024点波形图以及幅频特性图三.仿真要求及思考题1.任务1绘制时域波形、幅频特性曲线,观察幅频特性曲线2.任务2频谱特性曲线及分析信号x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)(1)问信号x中有无直流信号,交流分量的幅度、频率、相位分别是多少?(2)若用fs=100Hz频率对其采样,请分别绘制N=128、1024点波形图以及幅频特性图(3)根据幅频特性,分析频率组成以及幅度只分析fs/2以下的频谱(4)分析128点和1024点在频谱分辨率方面的区别音频文件(1)时域波形(2)频域波形(3)连续信号(4)采样(5)波形及频谱图(6)四.仿真实现步骤、代码及结果实验任务1:框(1)、(2)和(3)组成,对音频“F215.wav”分别绘制时域波形、幅频特性曲线代码:[y,fs,bits]=wavread( F215.wav );N=length(y);n=0:N-1;t=n/fs;subplot(2,1,1);plot(t,y);title( audiosignalintimedomain );xlabel( Time(second) );ylabel( Magnitude );subplot(2,1,2);f=fs*n/N;y_freq=abs(fft(y));plot(f,y_freq);title( signalinfrequencydomain );xlabel( Time(second) );ylabel( Amplitude );结果图00.10.20.30.40.50.60.7-0.4-0.200.20.4audiosignalintimedomainTime(second)Magnitude00.511.522.5x1040100200300signalinfrequencydomainTime(second)Amplitude实验任务2:框(4)、(5)和(6)组成,对连续信号用fs=100Hz频率对其采样,分别绘制N=128、1024点波形图以及幅频特性图信号x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)(1)问信号x中有无直流信号,交流分量的幅度、频率、相位分别是多少?答:没有直流信号,交流信号的幅度分别为0.5,2,频率分别为15,40,相位分别为0,0.(2)若用fs=100Hz频率对其采样,请分别绘制N=128、1024点波形图以及幅频特性图代码:n1=0:127;t1=n1/100;x=0.5*sin(2*pi*15*t1)+2*sin(2*pi*40*t1);subplot(2,2,1);plot(t1,x);title( audiosignalintimedomain );xlabel( Time(second) );ylabel( Magnitude );subplot(2,2,2);f=100*n1/128;y_freq1=abs(fft(x));plot(f,y_freq1);title( signalinfrequencydomain );xlabel( Time(second) );ylabel( Amplitude );n2=0:1023;t2=n2/100;x=0.5*sin(2*pi*15*t2)+2*sin(2*pi*40*t2);subplot(2,2,3);plot(t2,x);title( audiosignalintimedomain );xlabel( Time(second) );ylabel( Magnitude );subplot(2,2,4);f=100*n2/1024;y_freq2=abs(fft(x));plot(f,y_freq2);title( signalinfrequencydomain );xlabel( Time(second) );ylabel( Amplitude );结果图:00.20.40.60.811.21.4-3-2-10123audiosignalintimedomainTime(second)Magnitude020406080100020406080100120signalinfrequencydomainfrequency(hz)Amplitude024681012-3-2-10123audiosignalintimedomainTime(second)Magnitude0204060801000200400600800signalinfrequencydomainfrequency(hz))Amplitude(3)根据幅频特性,分析频率组成以及幅度只分析fs/2以下的频谱答:N=128时,f=15hz处幅度为35v,f=40hz处幅度为120v;N=1024时f=15hz处幅度为200v,f=40hz处幅度为800v。(4)分析128点和1024点在频谱分辨率方面的区别答:二者的频谱分辨率相同,与N的大小无关。