MATLAB第九次作业

发布 2022-09-20 20:58:28 阅读 3539

《工程应用数学》第九次作业。

化学工程张少石 622085216041

一、目的。使用dsolve函数求解微分方程的解析解,会使用ode函数求解常微分方程,熟悉ode函数调用的格式,了解各参数的意义。

二、思路。1. 微分方程如下:

dsolve的调用格式为:r=dsolve(‘eq1, eq2,… cond1, cond2, …v’)

2. 微分方程如下:

ode函数分为好几个在建立方程的基础上选择合适的ode方程求解可以快速得到比较满意的结果。

三、计算过程。

1. 调用dsolve求解上述方程的解析解。

> x=sym(‘t');

> y=dsolve('dx=x*(1-x^2)',x(0)=2','t')

y =1/(1-3/4*exp(-2*t))^1/2)

2.使用ode函数对该方程进行计算,将该方程编入m文件。

function dy=myfun01(t,y)

dy=zeros(2,1);

dy(1)=y(2);

dy(2)=(1/t)*y(1)-t*exp(t)*sin(t);

调用ode45函数来对上述方程求解其中积分区间为[-2 2],函数对应的初值为[2 0]。

> [t,y]=ode45(@myfun01,[-2 2],[2 0]);

> plot(t,y(:,1),'o')

对比使用ode45和ode15s解这个方程的结果如下。

t,y]=ode15s(@myfun01,[-2 2],[2 0]);

> hold on

> plot(t,y(:,1),'r-o')

四、结果与讨论。

对于有初值条件的微分方程使用dsolve函数可以快速的求得其解析解,但并不是所有的方程都能顺利接触其解析解,更多的时候只能获得其数值解。

对于一些实际问题往往只能获得数值解解释就需要使用ode函数对其求解,首先更具具体问题建立微分方程,其次要选择合适的ode函数进行求解。

医古文第九次作业及其答案

b型题 1.其病挛痹,其治宜微针。王冰注 微,细小也。王冰的注释属于 解释词义。2.灸之则瘖,石之则狂。王冰注 石,谓以石针开破之。王冰的注释属于 阐述语法。3.我国最早的字典是 说文解字 4.我国最早的词典是 尔雅 5.珠 字在 说文 中的部首是 玉部 6.音 字在 说文 中的部首是 音部 7.六...

医古文第九次作业及其答案

b型题 1.其病挛痹,其治宜微针。王冰注 微,细小也。王冰的注释属于 解释词义。2.灸之则瘖,石之则狂。王冰注 石,谓以石针开破之。王冰的注释属于 阐述语法。3.我国最早的字典是 说文解字 4.我国最早的词典是 尔雅 5.珠 字在 说文 中的部首是 玉部?6.音 字在 说文 中的部首是 音部?7.六...

消费心理学第九次作业

消费心理学第九次作业实战演练。为半明确型消费者建立购买标准。演练目的 为半明确型消费者建立购买标准。演练方法 角色扮演与情景表演。演练要求 1.4人一组,各扮演不同角色,如扮演购买者 扮演商家 扮演影响者等。2.设计某次买卖情景,如买房 买车 买手机 买电脑 买衣服等。3.根据消费者的情况,建立切实...