




下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、实验一:数值定积分求面积200820401015 代云雅一问题叙述用数值积分法求由y=-x2+115,y=0,x=0与x=10围成的图形面积,并讨论步长和积分方法对精度的影响。二问题分析用矩形法和梯形法分别求数值积分并作比较,步长的变化用循环语句实现。设x向量的长度为n,即将积分区间分为n-1段,各段长度为(i=1,2n-1)。算出各点的(i=1,2n+1),则矩形法数值积分公式为梯形法的公式为比较两个公式,它们之间的差别只是0.5()。三实验程序及注释clf,for dx=2,1,0.5,0.1 %设不同步长x=0:.1:10;y=-x.*x+115; %取较密的函数样本plot(x,y),
2、hold on %画出被积曲线,并保持x1=0:dx:10;y1=-x1.*x1+115; %求取样点上的y1n=length(x1);s=sum(y1(1:n-1)*dx; %用欧拉法求积分,末尾要去掉一点q=trapz(y1)*dx; %用梯形法求积分stairs(x1,y1),plot(x1,y1) %画出欧拉法及梯形法的积分区域dx,s,q,pause(10),hold off %显示步长及两种积分方法所得的面积rsums('115-x.2',0,10) %矩形法积分end四实验数据结果及分析ans =2 910 810ans =1 865 815ans =0.5000
3、 841.2500 816.2500ans =0.1000 821.6500 816.6500 步长dx矩形法解s梯形法解q29108101865815.5841.25816.25.1821.65816.65图1矩形法和梯形法的积分面积图2 矩形法积分的几何演示五实验结论 用解析法求得的精确解为2450/3=816.6666。dx=2时矩形法和梯形法的积分面积如图1。在曲线的切线斜率为负的情况下,矩形法的积分结果一定偏大,梯形法是由各采样点的连线包围的面积,在曲线曲率为负(上凸)时,其积分结果一定偏小,因此精确解在这两者之间。矩阵法积分的演示结果如图2。由此结果也能看出,步长相同时,梯形法的精
4、度比矩形法高。实验二:傅里叶级数200820401015 代云雅一问题叙述编写计算以x=-,为周期的任意函数的傅里叶系数的程序。二问题分析任何周期为2的满足狄利克雷条件的函数f(x),都可以用欧冠傅里叶级数表示为(-)其中:三实验程序及注释x=linspace(-pi,pi,1001);dx=2*pi/1000; %-pi,pi内长为1001的x数组及步长f=input('输入f=(长度为1001点的数组)'); %用户输入长为1001的f数组n=input('傅里叶系数的阶数n='); %用户给出所需要的阶数a0=trapz(f)/pi*dx %计算傅里叶系数
5、a0for k=1:n a(k)=trapz(f.*cos(k*x)/pi*dx; %计算傅里叶系数a(k) b(k)=trapz(f.*sin(k*x)/pi*dx; %计算傅里叶系数b(k) disp(k,a(k),b(k) %显示系数endpause,f1=a0/2*ones(size(x); %以a0为基础,构造傅里叶级数for k=1:n f1=f1+a(k)*cos(k*x)+b(k)*sin(k*x); %累加各项傅里叶级数endsubplot(1,2,1),plot(x,f),subplot(1,2,2),plot(x,f1) %在两个分图上画出四实验数据结果及分析 在输入f=
6、 x(1:501),zeros(1,500), 阶数n=9后,结果如下:a0 = -1.5708 1.0000 0.6366 1.0000 2.0000 -0.0000 -0.5000 3.0000 0.0707 0.3333 4.0000 -0.0000 -0.2500 5.0000 0.0255 0.2000 6.0000 0.0000 -0.1666 7.0000 0.0130 0.1428 8.0000 0.0000 -0.12509.0000 0.0079 0.1111f= x(1:501),zeros(1,500)ka(k)b(k)0-1.570810.63661.00002-0.0000-0.500030.07070.33334-0.0000-0.250050.02550.200060.0000-0.166670.01300.142880.0000-0.125090.00790.1111表1 函数所对应的傅里叶系数图1 信号的傅里叶变换前
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 通信工程基础知识
- 机械设备合作经营协议书范本
- 货仓库储存合同范本
- 联合开发合作协议书范本
- 个人委托公司代缴社保协议书范本
- 的工人工伤赔偿协议书范本
- 2025年度智能家居公寓租赁安全协议书
- 2025年度矿山工人安全赔偿保障协议
- 二零二五年度教育培训机构教育投资合作协议
- 2025年度股东经营协议书:智能化产业园区股权合作框架
- GB/T 3317-2006电力机车通用技术条件
- GB/T 30133-2013卫生巾用面层通用技术规范
- 二年级科学 《磁极与方向》优教
- 沥青路面病害课件
- 安全周知卡-酒精
- 《中学语文课程标准与教材研究》教学大纲
- 我国钢铁企业环境会计信息披露问题研究以宝钢为例13.26
- 测绘工程产品价格-国测财字20023号-测量费
- 罗氏试剂盒说明书 T3 11810456001V18
- 风机盘管机组巡检记录表
- 果汁加工工艺
评论
0/150
提交评论