实验四 连续时间信号的频域分析.docVIP

实验名称：连续时间信号的频域分析 学生姓名： 班级：08电科一 学号：指导老师：赵丽莉 同组人： 成绩： 实验目的: 分析傅里叶级数以及利用傅氏变换研究信号频谱。 实验要求： 编制Matlab程序，实现以下内容 1 已知门信号，求其傅里叶变换。 2 傅氏级数：将周期为，幅度为1的对称方波分解为多次正弦波之和，用Matlab演示谐波合成情况。（取到9次谐波） 实验方法（步骤）： 实验程序； 1·门信号 R=0.02; t=-2:R:2; f=heaviside(t+1)-heaviside(t-1); w1=2*pi*5; N=500; k=0:N; W=k*w1/N; F=f*exp(-j*t*W)*R; F=real(F); W=[-fliplr(W),W(2:501)]; F=[fliplr(F),F(2:501)]; subplot(2,1,1);plot(t,f); xlabel(t); ylabel(f(t)); title(f(t)=u(t+1)-u(t-1)); subplot(2,1,2);plot(W,F); xlabel(W); ylabel(F(w)); title(f(t)的傅氏变换F(w)); 2· t=0:.01:2*pi;%设定一个时间数组，有101个点 y=sin(t); plot(t,y),figure(gcf),pause %频率为w=1(f=1/2π)基波 y=sin(t)+sin(3*t)/3; plot(t,y),pause %叠加三次谐波 y=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7+sin(9*t)/9;%用1，3，5，7，9次谐波叠加 plot(t,y) %为了绘制三维曲面，要把各次波形数据存为一个三维数组，因此必须重新定义y y=zeros(10,max(size(t))); x=zeros(size(t)); for k=1:2:19 x=x+sin(k*t)/k; y((k+1)/2,:)=x; end %将各波形叠合绘出 pause,figure(1),plot(t,y(1:9,:)),grid line([0,pi+0.5],[pi/4,pi/4])%加上方波幅度线及标注 text(pi+0.5,pi/4,pi/4) %将各半波形绘成三维网格图，看出增加谐波阶次对方波逼近程度的影响 halft=ceil(length(t)/2);pause, figure(2),mesh(t(1:halft),[1:10],y(:,1:halft))%只用正半周波形 学生实验报告 批阅老师 ： 年 月 日 第 1 页 共3 页