- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
信号与系统 实验三 连续时间周期信号的傅里叶级数
实验三 连续时间周期信号的傅里叶级数
一、实验目的
掌握连续时间周期信号的傅里叶级数的展开和合成,理解吉布斯现象,掌握周期矩形脉冲信号的频谱及脉冲宽度、周期对周期信号频谱的影响。
二、实验内容
1、周期信号的傅里叶级数的展开和合成
画出如下图对称方波(取E=1、T=1),并采用有限项傅里叶级数对原函数进行逼近,画出对称方波的1、3、5、7、9、11次谐波的傅里叶级数合成波形,观察吉布斯现象。
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:1
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:3
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:5
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:7
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:9
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
sum=0;
t=-3:0.01:3;
E=1;T=1;ta=T/2;w=2*3.14159/T;
for n=1:11
fn=(2*E*ta/T)*sin(w*ta*n/2)/(w*ta*n/2);
f=(E*ta/T)+cos(n*w*t)*fn-E/2;
sum=sum+f;
end
plot(t,sum)
2、周期矩形脉冲信号的频谱
a. 取E=1,(=1, 画出周期矩形脉冲(教材P83图3-6)的傅里叶级数的频谱(教材P83图3-7);
n=-12:12;
E=1;t=1;T=5*t;w=2/T;
fn=(E*t/T)*sinc(w*t*n/2);
stem(n,fn,filled);
hold on
k=-12:0.01:12;
f=abs(E*t/T)*sinc(w*t*k/2);
plot(k,f,--);
b. 取E=1,(=1, 画出教材P85图3-8(a);
t=-12:0.01:12;
y=u(t+1/4)-u(t-1/4)+u(t-19/4)-u(t-21/4)-u(t+19/4)+u(t+21/4)+u(t-39/4)-u(t-41/4)-u(t+39/4)+u(t+41/4);
subplot(2,1,1);plot(t,y);axis([-12 12 -0.1 1.1]);xlabel(t);ylabel(f(t));
n=-12:12;
E=1;t=1;T=10*t;w=2/T;
fn=(E*t/T)*sinc(w*t*n/2);
subplot(2,1,2);stem(n,fn,filled);hold on;
k=-12:0.01:12;
f=abs(E*t/T)*sinc(w*t*k/2);
plot(k,f,--);
xlabel(w);ylabel(Fn);
c. 取E=1,(=1, 画出教材P85图3-8(c)。
t=-12:0.01:12;
y=u(t+1/4)-u(t-1/4)+u(t-39/4)-u(t-41/4)-u(t+39/4)+u(t+41/4);
subplot(2,1,1);plot(t,y);axis([-12 12 -0.1 1.1]);xlabel(t);ylabel(f(t));
n=-12:12;
E=1;t=1;T=5*t;w=2/T;
fn=(E*t/T)*sinc(w*t*n/2);
subplot(2,1,2);s
您可能关注的文档
最近下载
- SOAP(2023年)剖宫产术后快速康复共识.docx
- 网店运营推广 子任务3.3.1:商品标题制作与优化 3-3-1 商品标题优化.pptx
- 2024年北京市东城区八年级初二(下)期末道德与法治试卷(含答案).pdf
- 网店运营推广 子任务3.1.1:认识SEO工作原理 3-1-1 SEO工作原理.pptx
- 《第十二夜(莎士比亚剧本插图珍藏本)》—[英]莎士比亚著;朱生豪译.pdf
- 学习中国工会第十八次代表大会精神知识竞答活动试卷含答案.doc VIP
- WW-T 0016-2008馆藏文物保存环境质量检测技术规范.pdf VIP
- 2014年永固电力金具样本.pdf
- 食品制造公司安全风险分级管控清单.docx VIP
- 宏观经济学二十五讲中国视角.pptx VIP
文档评论(0)