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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 为了模拟，自学的matlab，心力憔悴，救救孩子。
52/1返回列表
查看: 2142  |  回复: 13
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者东翼将赠送您 35 个金币
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

东翼

禁虫 (小有名气)
本帖内容被屏蔽

» 猜你喜欢
1楼 2021-06-30 17:17:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

东翼

禁虫 (小有名气)
本帖内容被屏蔽
14楼2021-07-29 17:45:27
已阅   关注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

至尊木虫 (著名写手)
剩下的我也不会

发自小木虫Android客户端
回复此楼
3楼2021-07-01 15:37:56
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ofwhy

至尊木虫 (著名写手)
N和u有计算关系式吗？

发自小木虫Android客户端
赞一下 回复此楼
4楼2021-07-01 16:34:57
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 14 个回答
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 max7zaruFVS8 2026-10-06 3/150 2026-10-07 16:51 by kUQbiRsIeOpF
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 z1aZaKkVi71s 2026-10-06 3/150 2026-10-07 16:09 by kUQbiRsIeOpF
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 z1aZaKkVi71s 2026-10-06 3/150 2026-10-07 15:41 by kUQbiRsIeOpF
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 z1aZaKkVi71s 2026-10-05 3/150 2026-10-07 15:07 by kUQbiRsIeOpF
[考博] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 z1aZaKkVi71s 2026-10-05 3/150 2026-10-07 14:49 by kUQbiRsIeOpF
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 z1aZaKkVi71s 2026-10-05 3/150 2026-10-07 14:49 by kUQbiRsIeOpF
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 9/450 2026-10-07 11:07 by EXuX9l7TjjIY
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 8/400 2026-10-07 11:01 by EXuX9l7TjjIY
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 FlhJF7cq3uT3 2026-10-04 7/350 2026-10-07 10:51 by EXuX9l7TjjIY
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 FlhJF7cq3uT3 2026-10-04 8/400 2026-10-07 10:21 by EXuX9l7TjjIY
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 FlhJF7cq3uT3 2026-10-04 13/650 2026-10-07 10:09 by EXuX9l7TjjIY
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 FlhJF7cq3uT3 2026-10-04 7/350 2026-10-07 10:07 by EXuX9l7TjjIY
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 FUSxPqiIHRjH 2026-10-04 10/500 2026-10-07 10:01 by EXuX9l7TjjIY
[博后之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 FUSxPqiIHRjH 2026-10-04 11/550 2026-10-07 09:51 by EXuX9l7TjjIY
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 FUSxPqiIHRjH 2026-10-04 12/600 2026-10-07 09:49 by EXuX9l7TjjIY
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 FUSxPqiIHRjH 2026-10-04 12/600 2026-10-07 09:49 by EXuX9l7TjjIY
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 FUSxPqiIHRjH 2026-10-04 9/450 2026-10-06 14:28 by x7yLzftaCUEL
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 FlhJF7cq3uT3 2026-10-05 4/200 2026-10-06 05:08 by z1aZaKkVi71s
[教师之家] 烦躁 +5 水冰月月野兔 2026-10-03 5/250 2026-10-05 21:59 by 火星超人xi
[教师之家] 寻找合作团队 +4 悠悠2003 2026-10-01 4/200 2026-10-04 15:18 by Zh_fish
信息提示
关闭
请填处理意见
关闭