24小时热门版块排行榜

返回列表

【悬赏金币】回答本帖问题，作者东翼将赠送您 35 个金币

» 猜你喜欢

第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有239人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有17人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复

1楼2021-06-30 17:17:02

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

ofwhy

至尊木虫 (著名写手)

应助: 3 (幼儿园)
金币: 13429.6
散金: 93
红花: 7
帖子: 2309
在线: 177.9小时
虫号: 2777064
注册: 2013-11-04
性别: GG
专业: 工程热力学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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')