应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解非线性方程组
162/2返回列表上一页12
查看: 1031  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

月只蓝

主管区长 (职业作家)
引用回帖:
7楼: Originally posted by 夜雨tju at 2015-08-16 20:11:03
我是化工类的，不太擅长这个，大神可不可以指点细一些，我自己去查资料也可以，大神觉得什么函数可以，我之前也看了智能算法的书。

8楼的dingd是这类问题的大神，求助于他，能让你事半功倍。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
11楼2015-08-19 10:25:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜雨tju

金虫 (小有名气)
引用回帖:
9楼: Originally posted by dingd at 2015-08-16 22:38:40
抱歉，应该是“已知常数”

你好，我刚刚才看到你的回复，一直常数很多是公式算出来的，这里面我用了两个循环
赞一下 回复此楼
12楼2015-08-19 21:16:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜雨tju

金虫 (小有名气)
引用回帖:
9楼: Originally posted by dingd at 2015-08-16 22:38:40
抱歉，应该是“已知常数”

global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531;
yc1=0.92063;
yd1=0.03229;
T=370.72;
HETP1=0.4;
HETP2=0.2;
HETP3=0.42;
HETP4=0.2;
mi=803956.8;            
mig=922.45;            
M1=60;
M2=88;
M3=130;
M4=18;
Ts=40;               
lambda=286.76;        
wi1=0.5;                  
wi2=1;
Le=0.02772;
CA0=465.25;
CB0=475.76;
CC0=5644.11;
CD0=24.46;
k1=0.00000206;
K=3.621;
K1=0.85;
K2=0.9;
K3=0.95;
g=9.8;
for m=1:0.5:2;
    for R=0.01:0.01:0.4;
                S=pi*R^2;  
                a=acos((1+2*m)/(3+2*m));
                S1=a*R^2-((1+2*m)*(R^2)*sin(a))/(3+2*m);
                S5=a*R^2-((1+2*m)*(R^2)*sin(a))/(3+2*m);
                b=acos(1/(3+2*m));
                S2=(b-a)*R^2+((1+2*m)*R^2*sin(a)-(R^2)*sin(b))/(3+2*m);
                S4=(b-a)*R^2+((1+2*m)*R^2*sin(a)-(R^2)*sin(b))/(3+2*m);
                S3=(1-2*b/pi)*pi*R^2+(2*R^2/(3+2*m))*sin(b);
                U=(xa1*M1+xb1*M2+xc1*M3+xd1*M4)*L/(mi*pi*R^2);
                UG=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U1=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U3=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                U5=(ya1*M1+yb1*M2+yc1*M3+yd1*M4)*V/(mig*pi*R^2);
                A2=S2+S4;
                A4=S1+S3+S5;
                U2=S*U1/A2;
                U4=S*U1/A4;
                L1=3600*U;
                F=U2*sqrt(mig/1000);
                Lc1=Ts*U/2;
                dpz=(7.38*10^(4.46*L1*10^(-3)))*(UG*sqrt(mig/1000))^(1.72+L1*3.8*10^(-3));
                hL1=(dpz*Lc1*1000+lambda*(Lc1+Le)*U*mi-Lc1*mi)/(mi-dpz*1000);
                h1=Lc1+hL1;
                Lc2=Lc1;
                h3=h1;
                h2=HETP2*(0.5-h1/HETP1);
                h4=h2;
                lin=L1/S;
赞一下 回复此楼
13楼2015-08-19 21:17:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
13楼: Originally posted by 夜雨tju at 2015-08-19 21:17:31
global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531 ...

代码没贴全吧？
两个for 循环都不封闭？？
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
14楼2015-08-20 13:43:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cooooldog

铁杆木虫 (著名写手)

ส็็็
引用回帖:
13楼: Originally posted by 夜雨tju at 2015-08-19 21:17:31
global K CA0 CB0 CC0 CD0 k1 U L V xa1 xb1 xc1 xd1 ya1 HETP1 HETP2 HETP3 HETP4 mi M1 M2 M3 M4 K1
L=1.168;
V=0.7642;
xa1=0.07039;
xb1=0.07198;
xc1=0.85393;
xd1=0.0037;
ya1=0.02176;
yb1=0.02531 ...

看了你的代码，让人更加困惑了，你既然给常数赋值都是标量，为什么放在嵌套循环里面做无用功？

如果这段代码正确的写法是后面缺少两个end

那么你的代码相当于把 两个 end 去掉之后

for m=1:0.5:2;
    for R=0.01:0.01:0.4;

语句换成
m=2; R=0.4;

效果相同，而且少做了循环嵌套中的很多无意义的计算。

建议你不要用代码，而是用word文档和公式的形式，老老实实详细表达自己的问题，否则别人很难理解你的问题是什么
赞一下 回复此楼
ส็็็็็็็็็็็็็็็็็็็็
15楼2015-08-20 14:13:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

深有同感，自己问问题都描述不清或模棱两可，让他人还得化时间去猜测。求助应该先学会如何简单明了地把自己的问题说清楚，大家不可能像你一样了解该问题及其背景。
赞一下(1人) 回复此楼
16楼2015-08-20 16:00:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 夜雨tju 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭