24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知(长期有效)
查看: 980  |  回复: 25
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

lxsdkjdx

新虫 (初入文坛)

[交流] 运用

方程式

[ Last edited by lxsdkjdx on 2015-8-19 at 09:06 ]
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
lxsdkjdx: 金币+2 2015-08-18 11:16:03
引用回帖:
21楼: Originally posted by lxsdkjdx at 2015-08-18 11:02:38
我现在都不知道对不对...

要验证程序,实际运行一下就可以了。
a0~a7 b0~b7是随意赋值的,代码如下:
CODE:
function solve_odes
clear all;clc
global a0 a1 a2 a3 a4 b0 b1 b2 b3 b4 a5 a6 b5 b6 a7 b7

a0=0.1;
a1=0.2;
a2=0.3;
a3=0.4;
a4=0.5;
a5=0.6;
a6=0.7;
a7=0.8;

b0=0.7;
b1=0.6;
b2=0.5;
b3=0.4;
b4=0.3;
b5=0.2;
b6=0.1;
b7=0.1;



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

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

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

function f=odes(t,x)
global a0 a1 a2 a3 a4 b0 b1 b2 b3 b4 a5 a6 b5 b6 a7 b7

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';

计算结果:
CODE:
  
    0         0         0         0         0         0         0
    0.0202    0.0000    0.0001    0.0020    0.0141    0.0163    0.0020
    0.0404    0.0001    0.0006    0.0040    0.0283    0.0328    0.0041
    0.0606    0.0002    0.0013    0.0061    0.0424    0.0495    0.0062
    0.0808    0.0003    0.0023    0.0081    0.0566    0.0665    0.0083
    0.1010    0.0005    0.0036    0.0102    0.0708    0.0837    0.0105
    0.1212    0.0007    0.0051    0.0123    0.0850    0.1012    0.0127
    0.1414    0.0010    0.0070    0.0143    0.0992    0.1190    0.0150
    0.1616    0.0013    0.0092    0.0165    0.1135    0.1369    0.0173
    0.1818    0.0017    0.0116    0.0186    0.1278    0.1552    0.0196
    0.2020    0.0021    0.0143    0.0208    0.1421    0.1737    0.0220
    0.2222    0.0025    0.0173    0.0230    0.1564    0.1925    0.0244
    0.2424    0.0030    0.0206    0.0253    0.1708    0.2115    0.0268
    0.2626    0.0035    0.0242    0.0276    0.1853    0.2309    0.0293
    0.2828    0.0041    0.0281    0.0300    0.1998    0.2505    0.0318
    0.3030    0.0047    0.0323    0.0324    0.2144    0.2704    0.0344
    0.3232    0.0054    0.0368    0.0348    0.2290    0.2906    0.0370
    0.3434    0.0062    0.0416    0.0374    0.2437    0.3111    0.0397
    0.3636    0.0069    0.0466    0.0400    0.2585    0.3319    0.0424
    0.3838    0.0078    0.0520    0.0426    0.2733    0.3530    0.0451
    0.4040    0.0087    0.0577    0.0454    0.2882    0.3744    0.0479
    0.4242    0.0096    0.0636    0.0482    0.3032    0.3962    0.0508
    0.4444    0.0106    0.0699    0.0511    0.3183    0.4183    0.0537
    0.4646    0.0117    0.0765    0.0541    0.3336    0.4407    0.0566
    0.4848    0.0128    0.0834    0.0572    0.3489    0.4634    0.0596
    0.5051    0.0140    0.0906    0.0604    0.3643    0.4865    0.0627
    0.5253    0.0152    0.0981    0.0637    0.3798    0.5099    0.0658
    0.5455    0.0166    0.1060    0.0671    0.3955    0.5337    0.0689
    0.5657    0.0179    0.1141    0.0706    0.4113    0.5578    0.0722
    0.5859    0.0194    0.1226    0.0742    0.4272    0.5823    0.0754
    0.6061    0.0209    0.1314    0.0780    0.4433    0.6072    0.0788
    0.6263    0.0226    0.1405    0.0819    0.4595    0.6324    0.0822
    0.6465    0.0243    0.1499    0.0860    0.4759    0.6581    0.0856
    0.6667    0.0260    0.1597    0.0902    0.4925    0.6841    0.0892
    0.6869    0.0279    0.1698    0.0945    0.5092    0.7105    0.0927
    0.7071    0.0299    0.1803    0.0991    0.5261    0.7374    0.0964
    0.7273    0.0319    0.1911    0.1038    0.5433    0.7646    0.1001
    0.7475    0.0340    0.2022    0.1086    0.5606    0.7923    0.1039
    0.7677    0.0363    0.2137    0.1137    0.5781    0.8204    0.1078
    0.7879    0.0386    0.2256    0.1190    0.5959    0.8490    0.1117
    0.8081    0.0411    0.2378    0.1245    0.6139    0.8780    0.1157
    0.8283    0.0437    0.2504    0.1302    0.6322    0.9075    0.1198
    0.8485    0.0464    0.2634    0.1361    0.6508    0.9374    0.1240
    0.8687    0.0492    0.2767    0.1423    0.6696    0.9678    0.1283
    0.8889    0.0521    0.2904    0.1487    0.6887    0.9987    0.1326
    0.9091    0.0552    0.3045    0.1554    0.7082    1.0301    0.1370
    0.9293    0.0584    0.3190    0.1624    0.7279    1.0620    0.1415
    0.9495    0.0617    0.3340    0.1697    0.7480    1.0944    0.1462
    0.9697    0.0653    0.3493    0.1773    0.7685    1.1273    0.1509
    0.9899    0.0689    0.3650    0.1852    0.7894    1.1608    0.1557
    1.0101    0.0727    0.3812    0.1934    0.8107    1.1949    0.1606
    1.0303    0.0767    0.3978    0.2020    0.8324    1.2295    0.1656
    1.0505    0.0809    0.4148    0.2110    0.8545    1.2646    0.1707
    1.0707    0.0853    0.4323    0.2204    0.8772    1.3004    0.1759
    1.0909    0.0898    0.4502    0.2303    0.9004    1.3368    0.1813
    1.1111    0.0946    0.4687    0.2405    0.9241    1.3738    0.1867
    1.1313    0.0995    0.4876    0.2513    0.9483    1.4114    0.1923
    1.1515    0.1047    0.5070    0.2625    0.9732    1.4497    0.1980
    1.1717    0.1101    0.5269    0.2743    0.9988    1.4887    0.2039
    1.1919    0.1158    0.5474    0.2866    1.0250    1.5284    0.2099
    1.2121    0.1217    0.5683    0.2995    1.0520    1.5687    0.2160
    1.2323    0.1279    0.5899    0.3131    1.0797    1.6098    0.2224
    1.2525    0.1344    0.6120    0.3273    1.1083    1.6517    0.2288
    1.2727    0.1411    0.6347    0.3423    1.1378    1.6943    0.2355
    1.2929    0.1482    0.6579    0.3580    1.1683    1.7377    0.2423
    1.3131    0.1556    0.6819    0.3745    1.1997    1.7819    0.2493
    1.3333    0.1634    0.7064    0.3919    1.2323    1.8270    0.2565
    1.3535    0.1715    0.7317    0.4102    1.2661    1.8730    0.2639
    1.3737    0.1799    0.7576    0.4295    1.3011    1.9199    0.2715
    1.3939    0.1888    0.7842    0.4499    1.3375    1.9677    0.2794
    1.4141    0.1981    0.8116    0.4714    1.3754    2.0166    0.2875
    1.4343    0.2079    0.8398    0.4942    1.4149    2.0664    0.2959
    1.4545    0.2181    0.8688    0.5183    1.4562    2.1173    0.3045
    1.4747    0.2288    0.8987    0.5439    1.4993    2.1693    0.3135
    1.4949    0.2401    0.9294    0.5710    1.5445    2.2225    0.3228
    1.5152    0.2519    0.9611    0.5998    1.5919    2.2769    0.3324
    1.5354    0.2643    0.9937    0.6305    1.6417    2.3326    0.3424
    1.5556    0.2774    1.0274    0.6632    1.6943    2.3897    0.3528
    1.5758    0.2911    1.0622    0.6982    1.7499    2.4482    0.3637
    1.5960    0.3056    1.0982    0.7356    1.8087    2.5082    0.3750
    1.6162    0.3209    1.1353    0.7757    1.8713    2.5699    0.3869
    1.6364    0.3370    1.1738    0.8189    1.9379    2.6332    0.3994
    1.6566    0.3540    1.2137    0.8654    2.0089    2.6985    0.4125
    1.6768    0.3720    1.2550    0.9154    2.0848    2.7656    0.4263
    1.6970    0.3910    1.2979    0.9696    2.1663    2.8350    0.4409
    1.7172    0.4112    1.3425    1.0284    2.2541    2.9066    0.4563
    1.7374    0.4326    1.3890    1.0924    2.3492    2.9808    0.4728
    1.7576    0.4553    1.4375    1.1625    2.4527    3.0579    0.4904
    1.7778    0.4796    1.4882    1.2395    2.5656    3.1381    0.5093
    1.7980    0.5055    1.5413    1.3244    2.6895    3.2218    0.5297
    1.8182    0.5333    1.5971    1.4199    2.8283    3.3100    0.5521
    1.8384    0.5630    1.6558    1.5266    2.9828    3.4027    0.5767
    1.8586    0.5950    1.7177    1.6448    3.1532    3.5003    0.6033
    1.8788    0.6294    1.7830    1.7760    3.3421    3.6031    0.6325
    1.8990    0.6665    1.8524    1.9237    3.5541    3.7125    0.6647
    1.9192    0.7069    1.9264    2.0929    3.7963    3.8300    0.7009
    1.9394    0.7511    2.0058    2.2899    4.0782    3.9579    0.7425
    1.9596    0.7998    2.0916    2.5229    4.4115    4.0990    0.7909
    1.9798    0.8537    2.1850    2.8018    4.8102    4.2565    0.8480
    2.0000    0.9137    2.2870    3.1376    5.2907    4.4343    0.9161

运用
附图1.png

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
22楼2015-08-18 11:13:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 26 个回答

月只蓝

主管区长 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
Φ1c Φ4c的表达式是什么?

» 本帖已获得的红花(最新10朵)

MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
2楼2015-08-18 09:22:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lxsdkjdx

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by 月只蓝 at 2015-08-18 09:22:02
Φ1c Φ4c的表达式是什么?

这两个都是系数,
要求的就是运用四阶龙哥库塔求解x1,x2,x3,x4,x5,x6
3楼2015-08-18 09:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
引用回帖:
3楼: Originally posted by lxsdkjdx at 2015-08-18 09:27:00
这两个都是系数,
要求的就是运用四阶龙哥库塔求解x1,x2,x3,x4,x5,x6...

Φ1c Φ4c上面有一点的 也是系数?
MATLAB、MS小问题、普通问题请发帖求助!时间精力有限,恕不接受无偿私信求助。
4楼2015-08-18 09:28:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见