MAtlab傅里叶变换实验报告

1、matlab傅里叶变换实验报告 班级 信工 142 学号 22 姓名 何岩 实验组别 实验日期 室温 报告日期 成绩 报告内容:( 目得与要求, 原理, 步骤, 数据, 计算, 小结等) 1、求信号得离散时间傅立叶变换并分析其周期性与对称性; 给定正弦信号 x(t)=2s(2p*1),s=100z,求其tf。 ()代码: =10;=1f;w=1:0、2:10; t1=0:0、000:;0:0。0:1; n1=-2;2=8;=;nn1:0。0:n; 5=0.01; x=2cos(2*pi1); x2*co(2f*pit); x3=(exp().(t2*w); x=*x; subplot(2,2,

2、1);plo(t1,x1); axs( 1、1*min(x2) 1。1ma(2); xlae(x(n);ylabe(x(n)"); tit(原信号 x1"); lbe("t");ylabe("x1); sbpot(2,2,3);sem(t2,2); axis( 1 1、1min(x2) 1。max(x2); title(原信号采样结果 x2); xlabel();lael(x2"); ubplot(2,2,);ste(n,x5); axi(0 1 1、1*mi(x5) 1.1*mx(x5); xab(n);abel(x2")

3、; titl(采样函数); ublt(2,2,4);tm(t2,4); axis(0 1 0、21。1*min(x) 、1*x(x4); xlabel(t");ylabel(x4"); title("f结果 x4); (b)结果: 、用以下两个有限长序列来验证 dtft 得线性、卷积与共轭特性; x1(n)= 3 5 6 7 8 9 10 11 ;x2()=r 10 (n) (1) 线性:(a)代码: w=linae(-8,8,10000); 1:11; nx20:9; 1=1 2 4 8 9 10 1 ; x2=1 1 1 1 1 1 ; x3=x2,es(,(

4、legth(x1)lngth(x2); x4=2x1+3x3; 1=x1*ex(-jx1w);频率特性 x3=x*exp(-j*1);%频率特性 x44*exp(*1*w);频率特性 suplt(5,),em(nx1,x),xi(-1,3,0,15);tile(x), label("(n); suot(5,3,2),stem(nx,x),xis(,3,0,5);ite("x2); subplot(5,3),ste(nx1,x4),axis(-1,13,0,2);titl(x4=21+3"); sbplt(,3,4),(w,abs(x1); ylabl(幅度) su

5、bplot(5,7),lot(w,agle(1);abel(相位) subpt(5,10),plo(w,ral();yabel(实部) sblt(,3,13),pot(w,imag(x); ybl("虚部) uplo(5,3,5),pot(,as(x); ubpo(5,3,8),plot(w,ange(3); sublot(,3,11),plo(w,rel(x3); subplt(5,3,4),plot(w,imag(x3); subpl(,3,6),lot(,abs(x4); sbplot(5,3,9),plot(w,ane(); subpot(5,3,2),plo(w,ral(x

6、); sbplo(5,15),plot(w,mag(4); (b)结果: (2)卷积:(a)代码: nx=0:11; 2:9; x3=:20; w=nspace(-,40); w=8,8分 10000 份 x11 2 3 6 9 10 11 12; x1 1 1 1 1 1 1 1 1 1; x3=cnv(x,); x1 卷积 x2 4=x1exp(-nx1"*w); x频率特性 x5=2*exp(2w);% x2 频率特性 x6=x3*exp(-jn3"w);% x1 卷积 x频率特性 x7x、*x5; spt(2,2,),stm(x1,x1),ax(1,15,0,15)

7、,te(x1"); su plo (2,2,2), em(nx2, ),ax s(1, 5,0,5),title(x2); uplot(2,1,),stem(x3,3),axi(1,25,0,80);itle(x卷积2 结果 x3); figue,sublt(2,1),st(x,"filled),tite("得dtft 结果4); subplo(2,2,2),tem(x,"filled),titl(x得 dtft结果 x5); plot(2,3),st(6,fild),it(3得 dtft 结果 x6); sublot(2,2,4),te(x7,&quo

8、t;flled),itle(x4 得tft 结果7); figure,sulot(3,1),stem(,ab(6), ylabel("幅度),le(x1 卷积 x2 得 dtf); upl(4,2,3),ste(w,ale(6),ybe("相位") sbpo(4,2,5),stem(w,el(6),ylabl("实部) subpl(4,2,7),sem(w,mag(x6),ylabe(虚部) ublot(,2,2),stem(w,abs(x), ttl(x1 与 x2 得 dtf得乘积); sbplot(,2,4),tem(w,ngle(x); subp

9、lot(4,2,6),st(w,real(7); subpot(4,2,),s(w,ag(x7); ()结果: (3)共轭:(a)代码: x1n=1 2 3 4 5 6 7 9 1 1 12; w=10:10; n=ngth(1n);n=0:n11; x1=real(x1n); x2=ia(n); x2n=1jx2; 1x2n*(ex(j)、(n1); x2=xn*(exp(j)、(n1w); x3=ral(x2); ximg(2); =xj*x4; fgue,subpt(211);stem(w,1,".);ti("x1共轭得 dft); suplot(22);sem(,2

10、,"、);tle("1n 得 dft 取共轭且反折"); (b)结果: 3。 求 lti 系统得频率响应 给定系统 h()=b()/a(z),a0。8777 0。318 0、256 b=0.98997 0.989 。9897,求系统得幅频响应与相频响应、(要求使用ler(b,(n)求解。 (a)结果: a=0、98777 -。313 0、6; =0。899 0、89 0、899; c=1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 y=filtr(b,a,c); supot(2,2,1);stem(,、);tit(原始序列&qu

11、ot;); mag=abs(); =nle(); ph=ph*10i; sbplot(,2,2);m(m,"、);title(幅频特性); xlbel(时间信号");ylabel(信号幅度); suplt(2,2,);stm(ph,"、);itl("相频特性"); xla("时间信号 n);yabl("信号相位"); (b)结果: 4. 采样与频谱混叠 给定信号(t)=100*ep(-100*)co(2*i500t),求该信号得频谱;当采样频率分别为 fs=200z,f2=1000h;fs3=00hz; f=200

12、hz,时输出序列得 df。 (a)代码: x=100xp(-100*t)、cos(2p500t); =2:、1:2;w=10:0。1:10; y*(ex(-j)、(*); subpot(,1,1),plot(t,); ubpot(2,1,2),lt(w,);ttle(原始信号得频谱); figre,12021;ts1s1;n1=-2:t1:2; fs210;s2=1/fs;n2=2:ts2:2; fs3=500;ts3=1/f;n3-2:ts3:2; fs4=200;ts4=1/;n42:t4:2; x100。ep(10n1)。*cos(*pi*500*);=1(ep(-j)。(n1"

13、;w); sbplo(221);plt(w,);title("经 2021hz 采样后信号得 dtft"); x2=100。*ex(-0*2)、*os(*500n);2x2(e(-j)、(n2*w); subpl(222);ot(w,y);tle(经 1000h采样后信号得 dtt); x3=00、xp(00*n3)、co(pi50n3); 3x3*(ex(j)、("w); supt(22);lt(,y3);titl(经00hz 采样后信号得 dtf"); x4=100.exp(100*n4)。cs(*pi*5*4);y4=x*(exp()、(n4);

14、suplt(24);plo(w,y4);tile(经 2h采样后信号得 dtt); (b)结果: 收获及感想: df针对得就是有限长数字信号得傅立叶变换或傅立叶时频分析问题。但 以前得傅立叶变换就是定义在整个时间轴上得,而且一般针对得就是连续信号 ,获得得就是一个连续得频谱。 离散傅里叶变换(dft),就是傅里叶变换在时域与频域上都呈现离散得形式,将时域信号得采样变换为在离散时间傅里叶变换(dtt)频域得采样。在形式上,变换两端(时域与频域上)得序列就是有限长得,而实际上这两组序 列都应当被认为就是离散周期信号得主值序列。即使对有限长得离散信号作dt,也应当将其瞧作经过周期延拓成为周期信号再作变换。在实际应用中通常采用快速傅里叶