24小时热门版块排行榜

返回列表

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有16人回复
求个博导看看已经有19人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
matlab如何求解一个非线性微分方程组已经有8人回复
MATLAB求解非线性方程组已经有5人回复
求Matlab解非线性方程代码已经有12人回复
matlab解非线性方程组解不开已经有3人回复
二阶非线性微分方程求解已经有14人回复
matlab求解非线性方程组，错误提示怎么解决已经有5人回复
急切求助！一个mathematica求解一个简单一元非线性方程的问题已经有4人回复
matlab求解方程中的参数已经有21人回复
求助matlab---fsolve解非线性方程组已经有6人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
matlab求解非线性方程组，求助！已经有6人回复
matlab求解非线性方程组已经有16人回复
求matlab求解一个方程组的问题已经有6人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】MATLAB 有限差分法（FDM）求解偏微分方程已经有22人回复
【求助】积分微分方程matlab求解已经有6人回复
【求助】求MAtlab求解程序！！！已经有20人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2013-07-10 15:43:16

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

tgvictor99

木虫 (知名作家)

迷糊领域荣誉领主

应助: 40 (小学生)
金币: 2667.4
散金: 1907
红花: 49
沙发: 1
帖子: 5238
在线: 923.8小时
虫号: 2524003
注册: 2013-06-27
专业: 流体力学

【答案】应助回帖

感谢参与，应助指数 +1

妹纸，程序里的k2你写错了，应该是除以p0.,你怎么写成了k2=sqrt(r1*Ec*E0/x);应该是：
k2=sqrt(r1*Ec*E0/p0);

赞一下

回复此楼

每个人的人生都是时代的产物。每个人的血统、命运和梦想组成了现在这个时代。

2楼2013-07-10 17:11:53

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

月只蓝

主管区长 (职业作家)

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

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zyxzyh921: 金币+20, ★★★很有帮助, 灰常感谢~ 问题解决了~~~ 2013-07-11 10:03:59

1.解x有虚数问题
解答：把x的初值取一点，取0.01，就不会出现虚数了：

function zyxzy
clear all;clc
format long
x0=0.01;
options=optimset('Display','iter'); %创建options结构变量，其中display的值为iter，
[x,fval]=fsolve(@myfun,x0,options)

function f=myfun(x)
r1=0.001;
rn=0.02;
Ec=4653000;
E0=8.85e-12;
U=25000;
k1=sqrt(r1*Ec*x/E0);
k2=sqrt(r1*Ec*E0/x);
k3=sqrt(k2^2-r1^2);
f1=sqrt(rn^2+k2^2-r1^2);
f=U-k1*(f1-k2+k3*(log(rn/r1)+log(k2+k3)-log(k3+f1)));

结果：

                                       Norm of    First-order Trust-region
Iteration  Func-count    f(x)       step       optimality radius
   0       2 4.59128e+008                   4.26e+010             1
   1       3 4.59128e+008    0.0107764    4.26e+010             1
   2       5    2.462e+008    0.0026941    3.6e+010       0.00269
   3       7    4.3381e+007    0.00673525    4.09e+010       0.00674
   4       9 1.46051e+006    0.00106182    5.2e+009       0.00674
   5       11       1394.16 0.000280992    1.51e+008       0.00674
   6       13    0.00124421 9.23407e-006    1.42e+005       0.00674
   7       15 4.93763e-016 8.73987e-009       0.0897       0.00674
   8       17             0 5.50577e-015             0       0.00674

Equation solved.

fsolve completed because the vector of function values is near zero
as measured by the default value of the function tolerance, and
the problem appears regular as measured by the gradient.

<stopping criteria details>

x =

0.001922707332546

fval =

   0

2.把p0=0.0019,代回原方程，phi=91.8586的问题
解答：请注意你构建的myfun函数中出现了1e-12级别的极小数E0，如果p0取0.0019，仅有0.0001的精度，这样的精度与E0想乘除以后，误差会被剧烈放大，所以phi=91.8586
如果你代入x =

0.001922707332546
就会发现结果为：1.469743438065052e-009，已经非常接近于0了。

赞一下(2人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2013-07-10 17:26:43

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