[页数]:5 [字数]:537 [正文] 1.(1) 编制求第一类3次样条插值函数的通用程序。 (2) 已知汽车门曲线型值点的数据如下: ... 端点条件为 ,用编制的程序求车门的3次样条插值函数s(x),打印出S(i+0.5)(i=0,1,2,……,9)的值。 function f=spl() x=[0:1:10]; y=[2.51 3.30 4.04 4.70 5.22 5.54 5.78 5.40 5.57 5.70 5.80]; xx=[0.5:1:9.5]; yy=csapi(x,y,xx); f=yy; 2.给定初值 及容许误差ε,编制Newton发解方程 根的通用程序。 %其中fc为原函数,df为导数,x0为初值,q为允许误差,x1为根,n为迭代次数 function y=newton(x0,q) x1=x0-fc(x0)/df(x0); n=1; while (abs(x1-x0)>=q)&(n<=10000000) x0=x1; x1=x0-fc(x0)/df(x0);n=n+1; end x1 n 3. 编制RK4方法的通用程序,用编制的程序求解一个微分方程的初值问题。 %其中fun为微分方程,x0、x1为x的边界值,y0为初值,xn、yn为计算结果 %建立fun.m文件,若微分方程为y'=-2y+2x^2+2x,则fun.m文件如下: %fun.m %function f=fun(x,y) %f=-2*y+2*x.^2+2*x; function f=RK() [xn,yn]=ode45('fun',[x0,x1],y0); ...... [原文截取] 上机实验题 ——土木研0407 刘伟 (1) 编制求第一类3次样条插值函数的通用程序。 (2) 已知汽车门曲线型值点的数据如下: X 0 1 2 3 4 5 6 7 8 9 10 Y 2.51 3.30 4.04 4.70 5.22 5.54 5.78 5.40 5.57 5.70 5.80 端点条件为 ,用编制的程序求车门的3次样条插值函数s(x),打印出S(i+0.5)(i=0,1,2,……,9)的值。 function f=spl() x=[0:1:10]; y=[2.51 3.30 4.04 4.70 5.22 5.54 5.78 5.40 5.57 5.70 5.80]; xx=[0.5:1:9.5]; yy=csapi(x,y,xx); f=yy; 给定初值 及容许误差ε,编制Newton发解方程 根的通用程序。 %其中fc为原函数,df为导数,x0为初值,q为允许误差,x1为根,n为迭代次数 function y=newton(x0,q) x1=x0-fc(x0)/df(x0); n=1; while (abs(x1-x0)>=q)&(n<=10000000) x0=x1..... |
[数值分析] 数值分析上机实验报告
查看评论
已有0位网友发表了看法