24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1617  |  回复: 13
【悬赏金币】回答本帖问题,作者东翼将赠送您 35 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

东翼

禁虫 (小有名气)

本帖内容被屏蔽

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

东翼

禁虫 (小有名气)

本帖内容被屏蔽

12楼2021-07-29 14:31:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答

ofwhy

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
独孤神宇: 金币+5, 鼓励交流 2021-07-01 13:14:00
东翼: 金币+35, ★★★★★最佳答案, 大佬,和文献图一样啊,拜服。能不能再帮我看下第二问,下一个积分程序,怎么用到这个微分的内容。我的金币不够了,不能在帖子下面回复,太难了。 2021-07-01 14:57:10
clc;
clear all;
close all;

function du=f(xi,u)
%a替换alpha
global R V a;

du1=-u(1)-1./2*V*u(2);   %22 初始值
du2=-V*u(1)-(1+R*exp(-a*xi)*cosh(a*V *xi))*u(2)-(1./2*V-R*exp(-a*xi)*sinh(a*V *xi))*u(3);   %23;a替换alpha
du3=-(1./2*V-R*exp(-a*xi)*sinh(a*V *xi))*u(2)-(1+4*R*exp(-a*xi)*cosh(a*V *xi))*u(3)-(1./2*V-3*R*exp(-a*xi)*sinh(a*V *xi))*u(4);   %24
du4=-(1./2*V-R*exp(-a*xi)*sinh(a*V *xi))*u(3)-(1+9*R*exp(-a*xi)*cosh(a*V *xi))*u(4)';   %25
du=[du1;du2;du3;du4]

end

global R V a;
R=1;
V=1;
a=0.1;
[xi,u]=ode45('f',[0 10],[100 0 0 0]);   %求数值解
plot(xi,u(:,1),'-r',xi,u(:,2),'.k',xi,u(:,3),'*g',xi,u(:,4),':b')
2楼2021-07-01 11:32:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ofwhy

至尊木虫 (著名写手)

3楼2021-07-01 15:37:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ofwhy

至尊木虫 (著名写手)

4楼2021-07-01 16:34:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见