- 1、本文档共14页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
专业资料
word完美格式
目 录
TOC \o 1-2 \h \z 目录 1
1. 傅里叶级数 2
2.数学模型 2
2.1 理论分析 2
2.2 实例分析 2
3.程序探究 3
4.算法设计 4
5.程序设计 4
6.图形 界面设计 6
6.1 界面设计 6
6.2 主程序代码 6
7测试数据及结果 11
8. 改进之处 12
9.总结 12
参考文献 14
关于傅里叶级数的课程设计
一.傅里叶级数
1804年,傅里叶提出“在有限区间上由任意图形定义的任意函数可以表示为单纯的正弦和余弦函数之和”。傅氏级数的展开被称为最辉煌大胆的猜想。从分析的角度来看,一些复杂的周期现象用无限多个正弦函数余弦函数叠加来表示。从物理意义上讲,信号可以分解成为一系列的简谐波的复合,并可借由此来分析信号波的一些基本特征。此次课程设计的主要目的是用动画演示来说明当项数逐渐增大时,函数傅里叶级数的图像逐渐逼近原函数。
二. 数学模型
= 1 \* GB3 ①理论分析
数学上定义:若函数f(x)在区间[-π,π]可积,则称 (n=0,1,2,...), (n=0,1,2,3,...)是函数f(x)的傅里叶系数。以函数的傅里叶为级数的三角级数
当讨论到级数,不容忽视的是其部分和表示的是其前2n+1项的和。
若x是函数f(x)的第一类间断点,则函数f(x)的傅里叶级数收敛于函数f(x)在点x的左、右极限平均值,即,若x为连续点,则,则函数f(x)的傅里叶级数收敛于f(x)。
= 2 \* GB3 ②实例分析
将函数 展成傅里叶级数。
可求得
于是有
三.程序探究
对于这个特例函数展开成傅里叶级数,我首先想到的是作图表现其趋近程度。然而通过数学计算,可以得到这个函数的傅里叶级数展开形式,所以,没有通过计算机来算其傅里叶展开系数,直接用已知的函数g(x)进行绘图。程序如下:
function FULIYE(hedit,hlist)
n=str2num(get(hedit,String));%获取编辑框字符并转化为数字
n1=get(hlist,Value);%获取列表框选项序号
colmat=[g,r,y,k];%创建列表框取值的字符向量
axis([-4,4,-2,2])
x=-pi:0.005:pi;
g=0;
plot([0,4],[1,1],colmat(n1))
hold on
plot([0,-4],[-1,-1],colmat(n1))
for i=1:1:n
y=sin((2*i-1)*x)/(2*i-1);
g=g+(4/pi)*y;
plot(x,g,colmat(1+mod(i,3)))
pause(1)
end
hold off
四.算法设计
现在由特例的函数向一般化转换,当然问题只能一步一步来,先在[-pi,0]和[0,pi]分别输入函数(向量形式),然后使用积分函数,算出a0,即初始化S,然后进入循环,不断重复积分算出a(i),b(i),可由其构造出关于S的表达式,再将x从[-pi,pi]的函数用图画出。
第1步 获取由用户输入的向量形式的函数,转化成在线函数。同时获取用户输入的项数n
第2步 将x符号化,并用int函数积分算出s,即对s的初始化。
第3步 对于i=1,2,3,...n,反复做以下操作
= 1 \* GB2 ⑴ 将函数积分算得a(i)和b(i)
= 2 \* GB2 ⑵ 令x从-pi到pi,计算s的值,并将s转化为数值形式
= 3 \* GB2 ⑶ 作图画出s的图像,并将图形保持
= 4 \* GB2 ⑷ 再将x符号化,以便再次积分
= 5 \* GB2 ⑸ 每次画图需暂停来观察效果
第4步 循环结束
五.程序设计
这里将原函数和级数展开的函数分开,避免其中的语句相互影响,也使程序更加简明清晰。
= 1 \* GB4 ㈠原函数作图的函数
function ORIGIN(hlist,hedit1,hedit2)
n1=get(hlist,Value);%获取列表框选项序号
colmat=[g,r,y,k];%创建列表框取值的字符向量
axis([-4,4,-2,2])
x=-pi:0.005:0;
y1=inline(get(hedit1,String));
g1=y1(x);
plot(x,g1,colmat(n1))
x=0:0.005:pi;
y2=inli
文档评论(0)