数值逼近课程设计 080810212 傅强.doc

数值逼近课程设计 姓名： 傅强 学号： 080810212 专业： 信息与计算科学 编程语言： Matlab 南 京 航 空 航 天 大 学 2009年11月 第一题 施密特正交 这是一种施密特正交化的改进，不容易产生在传统施密特正交化中产生的数值不稳定问题。 编写m文件——smit.m function y=smit(A,min,max) s=size(A); n=s(1,2); for m=1:n f=A(1,m); f=f/sqrt(int(f*f,min,max)); B(m,1)=expand(f); A(1,m)=f; i=m+1; while(i=n) f1=A(1,i); A(1,i)=f1-int(f*f1,min,max)*f; i=i+1; end end y=B; 在命令窗口输入命令： clear syms t A=[1,t]; smit(A,0,1) %求在区间[0,1]上关于[1，t]的一组标准正交基 ans = 1 2*3^(1/2)*t-3^(1/2) A=[1,t,t^2]; smit(A,0,1) %求在区间[0,1]上关于[1,t,t^2]的一组标准正交基 ans = 1 2*3^(1/2)*t-3^(1/2) 6*5^(1/2)*t^2+5^(1/2)-6*5^(1/2)*t 第二题最小二乘法 为了制定生产计划，某羊毛衫厂记录了一部一年来羊毛衫的销售量，按月份得到表（1）的数据表，其销售量单位为箱，试建立月份（x）和销售（y）之间的关系。 月份（x） 1 2 3 4 5 6 7 8 9 10 11 12 销量（y） 256 201 159 61 77 40 17 25 103 156 222 345 表（1） 输入a=[1 1 1;1 2 4;1 3 9; 1 4 16;1 5 25;1 6 36;1 7 49;1 8 64;1 9 81;1 10 100;1 11 121;1 12 144]; b=[256;201;159;61;77;40;17;25;103;156;222;345]; x=inv(a*a)*a*b x = 386.0000 -113.4266 9.0420 x1=1:0.01:12; y=9.0420.*x1.^2-113.4266.*x1+386; plot(x1,y) hold on x1=[1,2,3,4,5,6,7,8,9,10,11,12]; y=[256,201,159,61,77,40,17,25,103,156,222,345]; plot(x1,y,.) 第三题 龙格现象 函数定义在[-1,1]上，这是一个很光滑的函数，它的任意阶导数都存在。但对它在[-1,1]上做等距节点插值时，插值多项式就会出现图（1）的情况。 function f=runge(n) syms t z z1; x=-1:2/n:1; y=1./(1+25*x.^2); z=0; for i=1:(n+1) z1=y(1,i); for j=1:(n+1) if(i~=j) z1=z1.*(t-x(1,j))./(x(1,i)-x(1,j)); end end z=z+z1; end f=expand(z); 在命令窗口输入命令： x=-1:0.01:1; y=(1+25*x.^2).^(-1); plot(x,y) hold on y1=runge(4); y1=1250/377*t^4-3225/754*t^2+1; %n=4时的插值多项式 y2=runge(8); y2=1+228601250/3725137*t^4-383000000/3725137*t^6+200000000/3725137*t^87450274*t^2; %n=8时的插值多项式 y3=runge(12); y3=1+367051586875/1847048164*t^4-107641853578125/112669938004*t^6+62017871484375/28167484501*t^8+25628906250000/28167484501*