应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 将一个微分方程组的解应用到另一个方程
11/1返回列表
查看: 633  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

m76acxb

铜虫 (小有名气)

[求助] 将一个微分方程组的解应用到另一个方程

第一个方程组M文件：
CODE:
function dy=eq3(t,y)
dy=zeros(2,1);
dy(1)=[1.6+0.2*(1-y(1))^2+0.8*sin(200*pi*t)]/[1.6+0.4*(1-y(1))^2];
dy(2)=0.2*y(2)*[0.8+0.8*sin(200*pi*t)]*(y(1)-1)/[0.8+0.2*(1-y(1))^2]^2;

第二个方程M文件，x即dy(1)，y即dy(2)，将x，y代入到p的方程，最后画plot（x，p）
CODE:
function f=fx(x,y);
f=(0.8+0.2*x.^2)^1.5*x.^2/y;
p=log(0.8+0.2*x.^2)-0.2857*log(0.0351+0.01348*quad(@(x)fx(x,y),x,1));%第二个方程，里面含有积分，积分范围【x 1】
H0=1;
[t,y]=ode45('eq3',[0 1],[0 H0]);
plot(x,p)

1.程序运行有错（方程输入没错），将第一个M文件的解代入到第二个M文件的P的方程，
2.第二个M文件中的x即dy(1)，y即dy(2)，然后画p和x的关系，希望有人帮忙修改一下
回复此楼

» 猜你喜欢
1楼 2015-06-05 15:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 m76acxb 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭