24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lxsdkjdx

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁

[交流] 运用

方程式

[ Last edited by lxsdkjdx on 2015-8-19 at 09:06 ]

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有8人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有5人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

1楼2015-08-16 19:52:29

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

lxsdkjdx

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁

引用回帖:

17楼: Originally posted by 月只蓝 at 2015-08-18 10:41:12
只要忽略掉积分式，程序代码就和之前帖子中的类似，常数是没有影响的。...

function solve_odes
clear all;clc
global a0 a1 a2 a3 b0 b1 b2 b3

x0=[0 0 0 0 0 0];
tspan=linspace(0,2,100);

[t x]=ode45(@odes,tspan,x0);
[t x]
figure(1)
subplot(2,2,1);
plot(t,x(:,1));
subplot(2,2,2);
plot(t,x(:,2));
subplot(2,2,3);
plot(t,x(:,3));
subplot(2,2,4);
plot(t,x(:,4));

subplot(2,2,5);
plot(t,x(:,5));
subplot(2,2,6);
plot(t,x(:,6));

function f=odes(t,x)
global a0 a1 a2 a3 a4 b0 b1 b2 b3 b4

x1=x(1)
x2=x(2)
x3=x(3)
x4=x(4)
x5=x(5)
x6=x(6)

f(1)=x3;
f(2)=x4;
f(3)=a0+a1*x3^2+a2*x3*x4+a3*x4^2;
f(4)=b0+b1*x3^2+b2*x3*x4+b3*x4^2;
f(5)=a7+a4*x3^2+a5*x1+a6*x5;
f(6)=b7+b4*x3^2+b5*x1+b6*x5;

f=f'

这样行不行