7Р5.4Р试建立该地全年温度与时间的经验公式。Р程序:Рx=1:1:12;Рy=[3.1,3.8,6.9,12.7,16.8,20.5,24.5,25.9,22,16.1,10.7,5.4];Рplot(x,y,'k.','markersize',25);Рp3=polyfit(x,y,10)Рt=1:1:12;Рs=polyval(p3,t);Рhold onРplot(t,s,'k-','linewidth',2)РgridР结果:Р即y=-0.003*x^9+0.0064*x^8-0.0811*x^7+0.5763*x^6-2.0422*x^5+0.8693*x^4+19.3461*x^3-64.7901*x^2+82.4241*x-33.2083;Р实验十四水塔水流量估计模型实验Р练习一(1)Р题目:Р对函数fx=11+25x2在区间[-10,10]上做等距拉格朗日插值法和分段线性插值,观察插值中出现的龙日现象,要求:在区间[-10,10]上取不同的插值节点数n,做出函数fx的拉格朗日插值多项式和分段多项式;将曲线fx和拉格朗日插值多项式曲线以及分段插值曲线绘在同一坐标系上进行比较。Р程序:(1)Р拉格朗日插值程序:Рfunction p=lagrange(x,y)РL=length(x);РA=ones(L); %产生元素全是1的L阶方阵Рfor j=2:LР A(:,j)=A(:,j-1).*x';РendРX=inv(A)*y';Рfor i=1:LР p(i)=X(L-i+1);РendР拉格朗日插值:Рx=-10:1:10;Рy=(1.0+25.*x.^2).^(-1);Рplot(x,y,'b.','markersize',30)РgridРhold on Рp=lagrange(x,y);Рt=-10:0.01:10;Рu=polyval(p,t);
全文预览
