信号与系统matlab实验傅里叶分析及应用报告答案
更新时间:2024-04-30 06:25:01 阅读量: 综合文库 文档下载
实验二 傅里叶分析及应用
姓名学号班级
一、实验目的
(一)掌握使用Matlab进行周期信号傅里叶级数展开和频谱分析
1、学会使用Matlab分析傅里叶级数展开,深入理解傅里叶级数的物理含义 2、学会使用Matlab分析周期信号的频谱特性
(二)掌握使用Matlab求解信号的傅里叶变换并分析傅里叶变换的性质
1、学会运用Matlab求连续时间信号的傅里叶变换 2、学会运用Matlab求连续时间信号的频谱图
3、学会运用Matlab分析连续时间信号的傅里叶变换的性质 (三) 掌握使用Matlab完成信号抽样并验证抽样定理
1、学会运用MATLAB完成信号抽样以及对抽样信号的频谱进行分析 2、学会运用MATLAB改变抽样时间间隔,观察抽样后信号的频谱变化 3、学会运用MATLAB对抽样后的信号进行重建
二、实验条件
需要一台PC机和一定的matlab编程能力
三、实验内容
2、分别利用Matlab符号运算求解法和数值计算法求下图所示信号的FT,并画出其频谱图(包括幅度谱和相位谱)[注:图中时间单位为:毫秒(ms)]。
1
符号运算法: Ft=
sym('t*(Heaviside(t+2)-Heaviside(t+1))+Heaviside(t+1)-Heaviside(t-1)+(-t)*(Heaviside(t-1)-Heaviside(t-2))'); Fw = fourier(Ft); ezplot(abs(Fw)),grid on; phase = atan(imag(Fw)/real(Fw)); ezplot(phase);grid on; title('|F|'); title('phase');
3、试用Matlab命令求F(j?)?
104-的傅里叶反变换,并绘出其时域信号图。
3?j?5?j?2
[注意:(1)写代码时j?i]
syms t
Fw = sym('10/(3+iw)-4/(5+iw)'); ft = ifourier(Fw,t); F = abs(ft);
ezplot(F,[-3,3]),grid on;
4、已知门函数自身卷积为三角波信号,试用Matlab命令验证FT的时域卷积定理。 [注:即验证门函数FT的平方与相应三角波信号的FT后结果是否一致,可结合频谱图观察分析]
3
f = sym('Heaviside(t+1)-Heaviside(t-1)'); Fw = simplify(fourier(f)); F = Fw.*Fw;
ezplot(abs(F)),grid on title('Fw*Fw');
triangle = sym('(t+2)*Heaviside(t+2)-2*t*Heaviside(t)+(t-2)*Heaviside(t-2)'); Fttriangle = fourier (triangle); F = simplify (Fttriangle); ezplot(abs(F),[-6,6]),grid on; title('triangle FT');
Hz;现在使用5、设有两个不同频率的余弦信号,频率分别为f1?100Hz,f2?3800Hz对这三个信号进行抽样,使用MATLAB命令画出各抽样信号的波形抽样频率fs?4000和频谱,并分析其频率混叠现象[建议:抽样信号的频谱图横坐标范围不小于-10000Hz~10000Hz或-20000*pi~20000*pi rad/s]。
F1=100Hz时,
实验代码如下:
Ts=1/4000;
4
dt=0.0001;
t1=-0.006:dt:0.006; ft=cos(2*pi*100*t1); subplot(221); plot(t1,ft),grid on;
axis([-0.006 0.006 -1.6 1.6]); xlabel('Time(sec)'),ylabel('f(t)') title('余弦信号波形');N=5000; k=-N:N;
W=2*pi*k/((2*N+1)*dt); Fw=dt*ft*exp(-j*t1'*W); subplot(222); plot(W,abs(Fw)); grid on;
axis([-20000 20000 0 0.006]); xlabel('\\omega'),ylabel('F(w)'); title('余弦信号的频谱'); t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2); subplot(223)
plot(t1,ft,':'),hold on stem(t2,fst),grid on
axis([-0.006 0.006 -1.5 1.5])
xlabel('Time (sec)'),ylabel('fs(t)') title('抽样后的信号'),hold off Fsw=Ts*fst*exp(-j*t2'*W); subplot(224)
plot(W,abs(Fsw)),grid on axis([-20000 20000 0 0.006]) xlabel('\\omega'),ylabel('Fsw') title('抽样信号的频谱')
5
F2=3800Hz时,
实验代码如下: Ts=1/4000; dt=0.0001;
t1=-0.006:dt:0.006; ft=cos(2*pi*3800*t1); subplot(221); plot(t1,ft),grid on;
axis([-0.006 0.006 -1.6 1.6]); xlabel('Time(sec)'),ylabel('f(t)') title('余弦信号波形');N=5000; k=-N:N;
W=2*pi*k/((2*N+1)*dt); Fw=dt*ft*exp(-j*t1'*W); subplot(222); plot(W,abs(Fw)); grid on;
axis([-20000 20000 0 0.006]); xlabel('\\omega'),ylabel('F(w)'); title('余弦信号的频谱'); t2=-0.006:Ts:0.006; fst=cos(2.*pi.*100*t2);
6
subplot(223)
plot(t1,ft,':'),hold on stem(t2,fst),grid on
axis([-0.006 0.006 -1.5 1.5])
xlabel('Time (sec)'),ylabel('fs(t)') title('抽样后的信号'),hold off Fsw=Ts*fst*exp(-j*t2'*W); subplot(224)
plot(W,abs(Fsw)),grid on axis([-20000 20000 0 0.006]) xlabel('\\omega'),ylabel('Fsw') title('抽样信号的频谱')
(六) 结合抽样定理,利用MATLAB编程实现Sa(t)信号经过冲激脉冲抽样后得到的抽样信号fs?t?及其频谱[建议:冲激脉冲的周期分别取4*pi/3 s、pi s、2*pi/3 s三种情况对比],并利用fs?t?构建Sa(t)信号。
周期取4*pi/3 s:
7
周期取pi s:
wm=2;
wc=1.2*wm; Ts=pi; dt=0.1; t1=-4:dt:4;
ft=sinc(t1/pi).*(uCT(t1+10)-uCT(t1-10)); N=500; k=-N:N;
W=2*pi*k/((2*N+1)*dt); n=-100:100; nTs=n*Ts;
fst=sinc(nTs/pi).*(uCT(nTs+10)-uCT(nTs-10)); subplot(221);
plot(t1,ft,':'),hold on stem(nTs,fst),grid on axis([-4 4 -1 1]);
xlabel('Time (sec)'),ylabel('fs(t)') title('Sa(t)抽样后的信号'),hold off; Fsw=Ts*fst*exp(-j*nTs'*W); subplot(222)
8
plot(W,abs(Fsw),'c'),grid on axis([-10 10 -3 10])
xlabel('\\omega'),ylabel('Fsw') title('Sa(t)抽样信号的频谱') t=-10:dt:10;
f=fst*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t)))); subplot(223);
plot(t,f,'m:'),grid on; axis([-10 10 -2 9]);
xlabel('t'),ylabel('f(t)');
title('由f(nTs)信号重建得到Sa(t)信号');
周期取2*pi/3 s: 实验代码如下: wm=2;
wc=1.2*wm; Ts=2*pi/3; dt=0.1; t1=-4:dt:4;
ft=sinc(t1/pi).*(uCT(t1+10)-uCT(t1-10)); N=500; k=-N:N;
W=2*pi*k/((2*N+1)*dt); n=-100:100;
9
nTs=n*Ts;
fst=sinc(nTs/pi).*(uCT(nTs+10)-uCT(nTs-10)); subplot(221);
plot(t1,ft,':'),hold on stem(nTs,fst),grid on axis([-4 4 -1 1]);
xlabel('Time (sec)'),ylabel('fs(t)') title('Sa(t)抽样后的信号'),hold off; Fsw=Ts*fst*exp(-j*nTs'*W); subplot(222)
plot(W,abs(Fsw),'c'),grid on axis([-10 10 -3 10])
xlabel('\\omega'),ylabel('Fsw') title('Sa(t)抽样信号的频谱') t=-10:dt:10;
f=fst*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t)))); subplot(223);
plot(t,f,'m:'),grid on; axis([-10 10 -2 9]);
xlabel('t'),ylabel('f(t)');
title('由f(nTs)信号重建得到Sa(t)信号');
10
四、实验结论和讨论
信号的时域与频域呈离散与周期的对应关系。最低的抽样频率应该为2f,只要抽样间
隔不大于1/2t,时域中的波形就不会发生混叠。
五、实验思考
要时刻考虑使用的函数的参数列表的问题,这次试验一开始出现的问题就是参数不匹
配,本来应该输入字符串,我却输入了表达式,导致后来没办法画出函数图像。还有就是要选取适当的观察范围,有时候坐标长度选得不合适看不出结果。
11
四、实验结论和讨论
信号的时域与频域呈离散与周期的对应关系。最低的抽样频率应该为2f,只要抽样间
隔不大于1/2t,时域中的波形就不会发生混叠。
五、实验思考
要时刻考虑使用的函数的参数列表的问题,这次试验一开始出现的问题就是参数不匹
配,本来应该输入字符串,我却输入了表达式,导致后来没办法画出函数图像。还有就是要选取适当的观察范围,有时候坐标长度选得不合适看不出结果。
11
正在阅读:
信号与系统matlab实验傅里叶分析及应用报告答案04-30
小学教师性别结构问题研究04-03
通航安全技术论证的教学案例汇编01-31
人教新目标英语八年级下册Unit 5 单元教案 最新人教版初二英语下06-13
16春天大《钢结构设计原理》在线作业二05-05
格拉斯哥大学预科项目06-16
“遥观若画”前三岛04-06
2016-2022年中国脱硫设备产业运营格局态势及十三五发展策略研究06-10
第一章 医学影像设备学概论12-31
- 小学生造句大全
- 增压泵投资项目可行性研究报告(模板)
- 高中语文人教版粤教版必修1-5全部文言文知识点归纳
- 两学一做专题民主生活会组织生活会批评与自我批评环节个人发言提
- 管理处环境保洁工作操作标准作业指导书
- 2012六一儿童节活动议程 - 图文
- 移树申请报告
- 《贵州省市政工程计价定额》2016定额说明及计算规则
- 计算机长期没有向WSUS报告状态
- 汉语拼音教学策略研究
- 发展西部领先的航空货运枢纽
- 司法所上半年工作总结4篇
- 如何提高银行服务水平
- 发电厂各级人员岗位职责
- 丰田汽车的外部环境分析
- 2017—2018年最新冀教版四年级数学下册《混合运算》教案精品优质
- 中建八局样板策划 - 图文
- 戚安邦《项目管理学》电子书
- 2015年高级项目经理笔记
- 弯桥的设计要点
- 信号
- 答案
- 实验
- 报告
- 应用
- 分析
- matlab
- 系统
- 傅里
- 五四红旗团支部创建工作总结(彩色完整版)
- 站直了做人
- 2014科普知识竞赛300题
- 怎样编写校志
- 抗菌药物临床医师处方权限管理措施
- LSI 9364-8i配置Raid 5及单盘、多盘掉线处理 - 邱禹迪 - 图文
- 关于进一步明确在党员发展工作中实行公示制和票决制的通知
- 2012年基本能力测试模拟试卷三
- 产业及产业结构的概念
- 机械设计基础 - 第2版 - 朱龙英主编课后习题部分答案(有些题号
- 2011年最新公开选拔副科级领导干部试题及答案
- 2014年中石油思想政治理论考试单项选择题库答案
- 图书管理系统
- 安底镇关于开展用人单位遵守劳动用工和社会保险法律法规情况专项
- 2011年专业学位研究生入学统一考试-347-心理学专业综(精)
- 男装标准尺码对照表
- 信息技术会考-选择题(基础、多媒体)
- 黄石市换届政府工作报告
- 国民经济行业分类与代码(2011年11月采用)
- 优课《生活中的经济学》作业答案