电子科技大学数学实验4：数值计算实验

最新推荐文章于 2024-08-08 08:45:06 发布

夏夏夏至

最新推荐文章于 2024-08-08 08:45:06 发布

阅读量978

点赞数 5

文章标签：数学建模 matlab 开发语言

本文链接：https://blog.csdn.net/qq_45945307/article/details/131205366

版权

数值计算实验

基础训练

1.方程求根

编程调用fzero求解方程 $2x^3-3x^2+4x-5=0$ ，并将所求根赋给变量xp，编写一个函数调用fzero，并返回xp。

解：

function xp=myfun

f=@(x)(2*x.^3-3*x.^2+4*x-5);

xp=fzero(f,0);

ans =

  1.3711

2. 求解二阶微分方程

$\frac{d^2x}{dt^2}=20(1-x^2)dx/dt+0.5x\\ x(0)=2; x'(0)=0$

编写函数调用ode工具箱函数返回x在点0:0.1:5处的函数值，用列向量存储这些函数值。此列向量为double型数组.

解：

function z=test

y0=[2;0];

tn=[0:0.1:5];

[t,y]=ode23(@myfun1,tn,y0);

z=y(:,1);

function f=myfun1(t,x)

f=[x(2);20*(1-x(1).^2)*x(2)+0.5*x(1)];

3. 二次多项式拟合

某种产品在生产过程中的性能指标y与它所含的某种材料的含量x有关，现将试验所得16组数据记录列于下表。

x	20.05	22.09	24.13	26.24	28.11	30.29	32.09	34.23
y	26.5	10.46	2.75	3.53	11.67	29.98	52.26	87.19
x	36.23	38.2	40.27	42.27	44.07	46.05	48.47	50.08
y	128.11	176.24	235.17	300.25	365.66	445.1	552.84	631

要求拟合y与x的函数关系。用多项式拟合函数polyfit进行二次多项式拟合。编写函数文件返回2个参数：

第1个返回参数为二次多项式系数组成的行向量p（元素由高次到低次排列）；

第2个返回参数为拟合函数在点x=25:0.4:60处的函数值（用1个行向量表示）。

程序文件第1行参考格式如下：

function [p,v]= myfun

解：

function [p,v]=myfun

x=[20.05  22.09  24.13  26.24  28.11  30.29  32.09  34.23 36.23 38.2  40.27  42.27  44.07  46.05  48.47  50.08];

y=[26.5 10.46  2.75  3.53  11.67  29.98  52.26  87.19 128.11  176.24 235.17 300.25 365.66 445.1  552.84 631];

p=polyfit(x,y,2);

v=polyval(p,[25:0.4:60]);

P=1.0000 -49.9995 626.9857

V=1.0e+03 *[0.0020  0.0022  0.0026  0.0034  0.0046...1.1580  1.1854  1.2130]

综合训练

一．实验任务

请用Euler法和Matlab函数ode23求解下列微分方程：
$\frac{dy}{dt} = 0.02(1-0.001y)y\\ y(0)=10$

并将Euler求解结果与Matlab的ode23函数求解结果対比（精确解），并绘图。(t:[0,300])

二. 实验目的

认识Euler法。熟悉Matlab解微分方程数值解的函数.

三.实验过程

function test

%此处令h=1

h=1;

y=10;

Y1=[10];

for i=0:h:299

  y=y+h*myfun1(i,y);

  Y1=[Y1;y];

end

%下面是ode求解

[a,Y2]=ode23(@myfun1,[0:h:300],10);

%合并两种运算的结果求误差并输出

wc=abs(Y1-Y2);

all=[a,Y1,Y2,wc];

fprintf('    t   y_Euler y_exact  error\n')

disp(all)

%绘图

plot(a,Y1,'r',a,Y2,'b')

function f=myfun1(t,y)

f=0.02*(1-0.001*y).*y;