武漢科技大學《信號與系統(tǒng)》實驗四實驗報告參考模板

1/7專業(yè)班級電信(DB)1102班姓名XX學號 3505X驗時間2013年6月18日指導教師陳華麗成績1.掌握離散信號譜分析的方法：序列的傅里葉變換、離散傅里葉級數(shù)、離散傅里葉變換、快速傅里理解這些變換之間的關系；葉級數(shù)、離散傅里葉變換、快速傅里葉變換的Matlab實現(xiàn)；4.學習用FFT對連續(xù)信號和離散信號進行譜分析的方法，了解可能出現(xiàn)的分析誤差及其原因，以便1.對連續(xù)信號xa(t)=Ae-atsin(業(yè)0t)u(t)1 (A=444.128，a=502"，業(yè)0=502")進行理0.6想采樣，可得采樣序列a0.4圖1給出了xa(t)的幅頻特性曲線，由此圖可以確00100200300400500定對xa(t)采用的采樣頻率。分別取采樣頻率為f/Hz(4)觀察上述三種情況下，x(n)的幅度曲線是否一致？為什么？3.(1)編制信號產(chǎn)生子程序，產(chǎn)生以下典型信號供譜分析用。"x(n)=cosn24"x(n)=sinn384(2)對信號x(n)，x(n)，x(n)進行兩次譜分析，F(xiàn)FT的變換區(qū)間N分別取8和16，觀察兩23 (3)連續(xù)信號x(n)的采樣頻率f=64Hz，N=16,32,64。觀察三次變換的結果是否一致？為什4s武漢科技大學《信號與系統(tǒng)》實驗四實驗報告functionx=ff(T)%定義函數(shù)文件n=0:50;A=444.128;a=50*sqrt(2)*pi;w0=50*sqrt(2)*pi;x=A*exp(-a*n*T).*sin(w0*n*T);k=-250:250;w=pi/125*k;X=x*(exp(-j*pi/125)).^(n'*k);實驗記錄及個人小結(包實驗源程注結果分析與討論Y=abs(y);plot(w/pi,Y)ndy=ff(1/1000);title('f=1000Hz')200HZ：y=ff(1/300);title('f=300Hz')7000600050004000300020000002000800004002000-1.5-1.5-0.5-0.5f=1000Hz0f=300Hz00.50.5111.5.522300HZ：y=ff(1/200);title('f=200Hz')00512武漢科技大學《信號與系統(tǒng)》實驗四實驗報告functiony=x(n)%定義函數(shù)文件y=cos(0.48*pi*n)+cos(0.52*pi*n);end(1)n=0:10;f1=x(n);f2=(fft(f1));stem(n,f2)title('x(n)(n=0:10)')axis([0,10,-2,8])8765432100123x(n)(n=0:10)456789(2)f3=zeros(1,20);forn=0:10a=2*n+1;f3(a)=n;endf4=x(f3);f5=fft(f4)stem(f5)8642005x(n)(n=0:20)2025結(包括：(3)n=0:100;f1=x(n);f2=(fft(f1));stem(n,f2)002030405060708090x(n)的幅度曲線是否一致(1)N=8;n=0:N-1;x1=(n+1).*(u(n-4)-u(n))+(8-n).*(u(n-8)-u(n-4));x2=cos(pi*n/4);x3=sin(pi*n/8);%x4=cos(8*pi.*t)+cos(16*pi.*t)+cos(20*pi.*t);y1=fft(x1,N);subplot(2,2,1);stem(n,abs(y1));y2=fft(x2,N);subplot(2,2,2);武漢科技大學《信號與系統(tǒng)》實驗四實驗報告y2=fft(x2,N);subplot(2,2,2);stem(n,abs(y2));y3=fft(x3,N);subplot(2,2,3);stem(n,abs(y3));50064200224466884321002468(2)行周期延拓后的傅里葉級數(shù)變換的主值系列,里葉級數(shù)也不相同,所以變換后的結果也不相同.(3)n=0:15;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);x4k16=fft(x4,16);結(包括：subplot(3,1,1)stem(n,abs(x4k16))title('x4FFT變換N=16')axis([0,14,0,15])n=0:31;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);xsubplot(3,1,2)5stem(n,abs(x4k5title('x4FFT變換N=32')axis([0,31,0,20])0204681020551015510152025000武漢科技大學《信號與系統(tǒng)》實驗四實驗報告n=0:63;f=64;x4=cos(8*pi*n/f)+cos(16*pi*n/f)+cos(20*pi*n/f);x4FFT變換N=x4FFT變換N=16subplot(3,1,3)stem(n,abs(x4ktitle('x4FFT變換N=64')axis([0,63,0,40])502406824068x4FFT變換N=3255列也不相同,所以傅里葉級數(shù)也不相同,所以變換后的結果也不相同.0051505x4FFT變