版块导航
正在加载中...
客户端APP下载
论文辅导
调剂小程序
登录
注册
帖子
帖子
用户
本版
应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！
24小时热门版块排行榜
>
论坛更新日志
(3209)
>
考研
(451)
>
虫友互识
(361)
>
导师招生
(332)
>
硕博家园
(97)
>
文献求助
(71)
>
论文投稿
(56)
>
休闲灌水
(46)
>
考博
(29)
>
博后之家
(22)
>
招聘信息布告栏
(17)
>
基金申请
(15)
>
教师之家
(14)
>
数学
(12)
>
公派出国
(11)
>
找工作
(9)
小木虫论坛-学术科研互动平台
»
计算模拟区
»
程序语言
»
其它
»
运用
26
3/3
返回列表
上一页
1
2
3
查看: 1088 | 回复: 25
只看楼主
@他人
存档
新回复提醒
(忽略)
收藏
在APP中查看
lxsdkjdx
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁
引用回帖:
19楼
:
Originally posted by
月只蓝
at 2015-08-18 10:56:56
两处 global a0 a1 a2 a3 a4 b0 b1 b2 b3 b4 中
把新加入的a5 b5等补上：
global a0 a1 a2 a3 a4 b0 b1 b2 b3 b4 a5 a6 b5 b6
两处以global引出的语句都要补...
我现在都不知道对不对
赞
一下
回复此楼
21楼
2015-08-18 11:02:38
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
月只蓝
主管区长
(职业作家)
专家经验: +1059
应助: 1712
(讲师)
贵宾: 8.888
金币: 68129.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖:
计算模拟区
★ ★ ★
小木虫: 金币+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的回帖
lxsdkjdx
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁
引用回帖:
22楼
:
Originally posted by
月只蓝
at 2015-08-18 11:13:26
要验证程序，实际运行一下就可以了。
a0~a7 b0~b7是随意赋值的，代码如下：
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 ...
谢谢你！！
回复此楼
23楼
2015-08-18 11:16:24
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
lxsdkjdx
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁
引用回帖:
6楼
:
Originally posted by
月只蓝
at 2015-08-18 09:34:16
那Φ1c Φ4c上面有一点具体表达式是什么？
要求x1,x2,x3,x4,x5,x6，6个方程6个未知因变量，可以计算，前提是方程中其他变量诸如带点的Φ1c Φ4c能写出用x1~x6表达的式子。...
就是用matlab编个方波
就是这个图形
赞
一下
回复此楼
24楼
2015-09-07 10:37:13
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
lxsdkjdx
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁
送红花一朵
引用回帖:
2楼
:
Originally posted by
月只蓝
at 2015-08-18 09:22:02
Φ1c Φ4c的表达式是什么？
亲，在不，
我想问你一个问题，就是我以前问过你的一个问题
还是用matlab解决龙哥库塔法求解四阶微分方程的
前一段时间找工作没做课题
最近几天我一直在调试，可是想得到一个近似的曲线
这个是以前你给我编写的程序
这个是结果
这个是那个要求解的方程
这个是我想得到的结果
这个是原方程。
赞
一下
回复此楼
25楼
2015-12-07 19:35:35
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
lxsdkjdx
新虫
(初入文坛)
应助: 0
(幼儿园)
金币: 10.6
散金: 8
帖子: 39
在线: 28.4小时
虫号: 3352479
注册: 2014-08-04
专业: 无机非金属类高温超导与磁
引用回帖:
2楼
:
Originally posted by
月只蓝
at 2015-08-18 09:22:02
Φ1c Φ4c的表达式是什么？
就是这一期的
回复此楼
26楼
2015-12-08 09:58:45
已阅
回复此楼
关注TA
给TA发消息
送TA红花
TA的回帖
相关版块跳转
第一性原理
量子化学
计算模拟
分子模拟
仿真模拟
程序语言
我要订阅楼主
lxsdkjdx
的主题更新
26
3/3
返回列表
上一页
1
2
3
如果回帖内容含有宣传信息，请如实选中。否则帐号将被全论坛禁言
普通表情
龙
兔
虎
猫
高级回复
(可上传附件)
百度网盘
|
360云盘
|
千易网盘
|
华为网盘
在新窗口页面中打开自己喜欢的网盘网站，将文件上传后，然后将下载链接复制到帖子内容中就可以了。
最具人气热帖推荐
[查看全部]
作者
回/看
最后发表
[
考研
]
一志愿085502，267分求调剂
+17
再忙也要吃饭啊
2026-04-08
18/900
2026-04-13 09:57
by
darren_yym
[
考研
]
302求调剂
+6
易！?
2026-04-13
6/300
2026-04-13 09:39
by
fqwang
[
考研
]
生物学调剂，一志愿西南大学348，Top期刊一区二作、二区三作，三等奖学金三次
+5
candyyyi
2026-04-09
5/250
2026-04-13 09:02
by
可淡不可忘
[
考研
]
314求调剂
+24
wakeluofu
2026-04-09
25/1250
2026-04-13 08:58
by
lhj2009
[
考研
]
291 求调剂
+30
化工2026届毕业�
2026-04-09
30/1500
2026-04-13 08:52
by
lhj2009
[
考研
]
071000生物学，一志愿深圳大学296分，求调剂
+13
TIckLw
2026-04-06
14/700
2026-04-13 08:49
by
可淡不可忘
[
考研
]
调剂求收留
+29
果然有我
2026-04-10
30/1500
2026-04-12 22:14
by
zxcwyt
[
考研
]
339求调剂
+3
hanwudada
2026-04-12
3/150
2026-04-12 21:43
by
lbsjt
[
考研
]
生物学308求调剂
+5
相信必会光芒万�
2026-04-11
5/250
2026-04-12 18:14
by
zhouxiaoyu
[
考研
]
299求调剂
+8
ZVVZ13
2026-04-08
8/400
2026-04-12 00:40
by
蓝云思雨
[
考研
]
288求调剂，一志愿华南理工大学071005
+18
ioodiiij
2026-04-08
18/900
2026-04-11 20:25
by
liyun12321
[
考研
]
调剂
+4
电气300求调剂不
2026-04-08
7/350
2026-04-11 10:44
by
紫曦紫棋
[
考研
]
284求调剂
+12
archer..
2026-04-10
13/650
2026-04-11 08:44
by
zhq0425
[
考研
]
309求调剂
+14
wdhw
2026-04-10
15/750
2026-04-10 21:06
by
zhouxiaoyu
[
考研
]
266求调剂
+29
阳阳哇塞
2026-04-07
29/1450
2026-04-10 16:20
by
高维春
[
论文投稿
]
mdpi小修rvr时间四五天了
20
+3
哈哈high
2026-04-08
5/250
2026-04-10 16:02
by
北京莱茵润色
[
考研
]
材料专硕初试分332一志愿西北工业大学，
+12
故人??
2026-04-09
12/600
2026-04-09 18:34
by
Ccclqqq
[
考研
]
一志愿华东理工085601材料工程303分求调剂
+15
a1708
2026-04-06
15/750
2026-04-08 16:23
by
luoyongfeng
[
考研
]
277、学硕，求调剂 数一104，
+11
瓶子PZ
2026-04-07
12/600
2026-04-07 23:30
by
一只好果子?
[
考研
]
327考研调剂推荐
+6
呜呜呜呜呢
2026-04-06
6/300
2026-04-06 21:39
by
啵啵啵0119
信息提示
关闭
请填处理意见
关闭
确定
回复此楼
送红花一朵 
20+3