24小时热门版块排行榜

返回列表

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

后会有期wq

新虫 (初入文坛)

应助: 2 (幼儿园)
金币: 640.5
帖子: 27
在线: 16.9小时
虫号: 3830458
注册: 2015-04-24
专业: 光学

[求助] matlab求解超越方程已有1人参与

程序如下：
p=0:0.01:100;ni=4e15*p;
% 超越方程 nf/ni=1+6.94e25/ni*log(nf/ni)
for i=1:length(ni);
nf=fsolve(@(x)x/ni-1-6.94e25/ni*log(x/ni),1.63e27);
end;

我想要得出nf和ni之间的关系，画出他们之间的曲线
用fsolve求解的话，总报错说我初值写的不对。
Objective function is returning undefined values at initial point. FSOLVE cannot continue.
不知道该怎么改了。求大神帮忙！

回复此楼

» 猜你喜欢

1楼 2016-05-27 20:44:30

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68161.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3292.2小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

CODE:

clear all;clc

format long

p=0.01:1:100;nni=4e15*p;

% 超越方程 nf/ni=1+6.94e25/ni*log(nf/ni)

for i=1:length(nni);

     ni=nni(i);

     [nf(i) fval(i)]=fsolve(@(x)x/ni-1-6.94e25/ni*log(x/ni),ni);

end;

[4e15*p ; nf ]'

fval'

figure(1)

plot(4e15*p,nf)

两个问题：
其一，如果 p=0:0.01:100;ni=4e15*p; 可知，第一个ni=0，在方程中有两处ni作分母，会出现 x/0 的情况，所以报错“Objective function is returning undefined values at initial point. FSOLVE cannot continue.”

其二，我给出代码可以运行，但实际上通过简单的分析，就会发现 nf=ni 就是原方程的解，即nf和ni的关系是y=x的直线关系，写代码已经没有意义。