| 查看: 1902 | 回复: 11 | |||
[求助]
牛顿迭代解方程画曲线 已有2人参与
|
|||
|
最近需根据一个方程计算用matlab来画一个图，但是总是不出结果，在此向大家求助。 方程如下图，方程中有3个符号变量，其中a和C0是通过前面的数值计算计算出来的，然后在a为某一值时，做出C0-Cw曲线。这里希望就是对a一定时，C0=0.1:0.1:10，分别求出对应的Cw的各个数据，而不是利用隐函数ezplot()来画图，如下图所示。 我现在编的若已知a和Co求Cw的话，可以做出来这单个的一组，但是一旦来做这条曲线时就怎么弄都弄不出来。希望大家帮忙看一下，指点一二（后面赋有一个matlab文件）。  公式.png  曲线.jpg |
回复此楼» 猜你喜欢
0703化学/290求调剂/本科经历丰富/工科也可
已经有11人回复
-
南京大学化学调剂
已经有15人回复
-
292求调剂
已经有4人回复
-
材料专硕调剂
已经有8人回复
-
085601材料工程找调剂
已经有17人回复
-
面上5B能上会吗？
已经有6人回复
-
化学308分调剂
已经有7人回复
-
化学工程085602 305分求调剂
已经有25人回复
-
求化学调剂
已经有11人回复
-
085602 化学工程专硕 340分求调剂
已经有3人回复
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9850.4
- 散金: 89
- 红花: 19
- 帖子: 3623
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
【答案】应助回帖
感谢参与，应助指数 +1
|
把你的程序稍微改动了一下，不多说，看code: clear; clc; a=0.6; syms Cw Cw0; C=0.1; i=1; while C<=10 f=Cw-1/(1-a^4)*(4/pi*(1/3*(2.5-1/(Cw+C)^2)*(1-1/(Cw+C)^2)^0.5+(Cw+C)/2*asin(1/(Cw+C)))-a^4*(Cw+C));%%牛顿迭代形式 df=diff(f,Cw); eps=1e-7; Cw0=10.0; cnt=0; MAXCNT=200; %最大循环次数 while cnt<MAXCNT %防止无限循环 Cw1=Cw0-subs(f,Cw,Cw0)/subs(df,Cw,Cw0); %去掉这个分号,可以看到迭代过程。 if (abs(Cw1-Cw0)<eps) break; end Cw0=Cw1; cnt=cnt+1; end error=abs(Cw1-Cw0)/2; %%误差 if cnt==MAXCNT || error > eps disp(['不收敛@C=',num2str(C)]); continue; end Y(i) =Cw0; i= i+1; C = C+0.1; end X = 0.1:0.1:10; plot(X,Y,'r','LineWidth',2); xlabel('{\itC}_0'); ylabel('{\itC}_w'); disp('done...') |
3楼2016-12-13 15:21:16
2楼2016-12-12 10:22:56
4楼2016-12-13 21:09:49
5楼2016-12-13 21:20:42
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9850.4
- 散金: 89
- 红花: 19
- 帖子: 3623
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
6楼2016-12-14 09:21:57
7楼2016-12-14 09:45:51
bluesine
铁杆木虫 (职业作家)
科苑小木虫
- 数学EPI: 5
- 应助: 132 (高中生)
- 贵宾: 1.991
- 金币: 9850.4
- 散金: 89
- 红花: 19
- 帖子: 3623
- 在线: 377.5小时
- 虫号: 869544
- 注册: 2009-10-12
- 性别: GG
- 专业: 数学物理
8楼2016-12-14 11:36:08
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
9楼2016-12-14 17:50:48
10楼2016-12-14 21:16:35
