【答案】应助回帖

» 猜你喜欢

导师想让我从独立一作变成了共一第一已经有9人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

1楼2016-12-12 10:18:36

回帖支持 ( 显示支持度最高的前 50 名 )

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9501.3
散金: 89
红花: 19
帖子: 3534
在线: 364.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...')

赞一下(1人)

板凳要做十年冷文章不发一个字

3楼2016-12-13 15:21:16

普通回帖

klding

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 96
散金: 3
帖子: 91
在线: 26小时
虫号: 4264065
注册: 2015-12-03
专业: 机构学与机器人

matlab程序。

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : SHIYAN.m

2016-12-12 10:19:40, 1.03 K

2楼2016-12-12 10:22:56

klding

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 96
散金: 3
帖子: 91
在线: 26小时
虫号: 4264065
注册: 2015-12-03
专业: 机构学与机器人

引用回帖:

3楼: Originally posted by bluesine at 2016-12-13 15:21:16
把你的程序稍微改动了一下，不多说，看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 ...

你好，很感谢你。但是我把你的程序导进去后发现进入死循环了，停不下来啊

4楼2016-12-13 21:09:49

Sky青云

银虫 (小有名气)

应助: 1 (幼儿园)
金币: 145.6
散金: 180
帖子: 154
在线: 45.6小时
虫号: 5082834
注册: 2016-10-03
专业: 金属结构材料

在学

青云志

5楼2016-12-13 21:20:42

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9501.3
散金: 89
红花: 19
帖子: 3534
在线: 364.5小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

【答案】应助回帖

引用回帖:

4楼: Originally posted by klding at 2016-12-13 21:09:49
你好，很感谢你。但是我把你的程序导进去后发现进入死循环了，停不下来啊...

并不是死循环，具体看电脑配置。。因为没有输出中间结果，需要等一会，我的大概在15秒左右出结果。看附件

test.jpg

板凳要做十年冷文章不发一个字

6楼2016-12-14 09:21:57

klding

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 96
散金: 3
帖子: 91
在线: 26小时
虫号: 4264065
注册: 2015-12-03
专业: 机构学与机器人

引用回帖:

6楼: Originally posted by bluesine at 2016-12-14 09:21:57
并不是死循环，具体看电脑配置。。因为没有输出中间结果，需要等一会，我的大概在15秒左右出结果。看附件

test.jpg
...

你好，还是不行啊，我的电脑配置不低啊，i7-4790了，8G内存啊，现在都运行看10分钟了，还是没有出结果，我其他的程序都好着，就这个不行，现在好像matlab也卡死了。
应该还是程序问题吧。

7楼2016-12-14 09:45:51

bluesine

铁杆木虫 (职业作家)

科苑小木虫

数学EPI: 5
应助: 132 (高中生)
贵宾: 1.991
金币: 9501.3
散金: 89
红花: 19
帖子: 3534
在线: 364.5小时
虫号: 869544
注册: 2009-10-12
性别: GG
专业: 数学物理

引用回帖:

7楼: Originally posted by klding at 2016-12-14 09:45:51
你好，还是不行啊，我的电脑配置不低啊，i7-4790了，8G内存啊，现在都运行看10分钟了，还是没有出结果，我其他的程序都好着，就这个不行，现在好像matlab也卡死了。
应该还是程序问题吧。...

你把循环里面C=C+0.1的分号去掉，还有中间的Cw1后面的分号去掉，看看有没有结果

板凳要做十年冷文章不发一个字

8楼2016-12-14 11:36:08