当前位置:首页 >> 理学 >>

数值分析第二次上机答案


第二次计算实习题 姓名郑延学号 1140509051
1、P95-1:对于给函数 f(x)=1/(1+25x2)在区间[-1,1]上取 xi=-1+0.2i(i=0,1,……10),试求 3 次曲线拟合,试画出拟合曲线并打印出方程,与第 2 章计算实习题 2 的结果比较。 程序: xi = -1 + (0:10) * 0.2; f = 1 ./ (1+25*xi.^2); p = polyfit(xi,f,3); x = -1 : 0.01 : 1; y1 = 1 ./ (1+25*x.^2); y2 = polyval(p,x); a=inv(g)*d p=a(1,1)+a(2,1)*x.^1+a(3,1)*x.^2+a(4,1)*x.^3 plot(xi,f,'bo',x,y1,'b--',x,y2,'r'); legend('节点','函数 f(x)','三次曲线拟合') 结果:

三次曲线方程: p=-4.4862e-17x3 - 0.57518x2 + 2.0447e-17x + 0.48412 结论: 可见三次曲线虽然不出现龙格现象,但拟合的效果不好,取样点的值与原函数误差 较大,不适合作为 f(x)的逼近。 2、补充:取切比雪夫多项式的零点作为取样点, =cos2( +1) ,k=1,2,…,n+1,
2 +1

n=20 程序:

k=0:1:20 m=2*k+1 x=cos(m*pi/42) y=1./(1+25*x.^2) D=vander(x); A=fliplr(D) B=[y].' C=inv(A)*B p=fliplr(C.') y2=polyval(p,x) plot(x,y,'-b',x,y2,'*k') grid on legend('切比雪夫正交多项式') title('y=1/(1+25*x2)') xlabel('x') ylabel('y')

——切比雪夫正交多项式

结论: 取切比雪夫正交多项式零点拟合的效果较好,且不出现龙格现象,可以采用。


赞助商链接
相关文章:
数值分析第二次上机作业
数值分析第二次上机作业一、切比雪夫零点插值多项式 1、Matlab 程序 clc clear all N=20;%等分段数 %插值点 j=1:N+1; x0(j)=cos((2*j-1)*pi/(2*...
数值分析第二次程序题
数值分析第二次程序题 - 数值分析第二次程序作业 数值分析第二次程序题——插值法 1.对 Runge 函数 R ( x ) ? 并对结果进行分析。 (1) 以 xi ? -1 ...
北航数值分析B第二次上机作业算法实现_图文
北航数值分析B第二次上机作业算法实现 - 数值分析 B 第二次上机作业 课后题第 8 题 关于 x , y , t , u , v , w 的方程组(A.3) ?0.5cost ? ...
数值分析上机作业第二次
数值分析上机作业第二次 隐藏>> 第二次上机作业数分五班 程序分配 1-2 题: 学生姓名:马羽龙 班学 3-4 题: 学生姓名:马纪翔 班学 5-6 题: 学生姓名:...
数值分析第二次作业
数值分析第二次上机作业算法设计 1. 矩阵 A 的拟上三角化 对矩阵 A 进行拟...数值分析第二次作业 36页 免费 数值分析第二次作业解答... 3页 免费 ...
清华大学贾仲孝老师高等数值分析第二次实验
清华大学贾仲孝老师高等数值分析第二次实验 - 高等数值分析第二次实验作业 高等数值分析第二次实验作业 高等数值分析第二次实验作业 T1.构造例子特征值全部在右半...
电子科大matlab与数值分析第二次上机实践报告
电子科大matlab与数值分析第二次上机实践报告 - Matlab 上机实践报告 实践内容:Simulink 的建模与仿真。 学院: 姓名: 学号: 指导老师: 实践日期: 题目(一) ...
清华大学高等数值分析_第二次实验作业
清华大学高等数值分析_第二次实验作业_数学_自然科学_专业资料。贾仲孝老师的第二次实验作业的答案 供大家参考!第二次实验作业 清华大学 高等数值分析 贾仲孝老师...
北航数值分析第二次大作业
北航数值分析第二次大作业_天文/地理_自然科学_专业资料。数值分析第二次大作业姓名:李潇 学号:SY1303314 题目:使用带双步位移的 QR 分解法求矩阵 A ? [aij ...
清华大学数值分析A第二次作业
清华大学数值分析A第二次作业_理学_高等教育_教育专区。1、1)证明: ||x ||...2013年注会设计统考真题及答案67份文档 九妖笑话 2014年笑话大全之让你笑个够 ...
更多相关标签: