您现在的位置:网站首页答辩论文论文专题

[数值分析] 数值分析上机实验报告

  • 简介:上机实验题——土木研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 端点条件为 ,用......
    • 请与管理员联系购买资料 QQ:5739126
  • 论文简介
  • 相关论文
  • 论文下载
[页数]: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位网友发表了看法
  • 验证码: