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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab求解非线性方程组
162/2返回列表上一页12
查看: 1148  |  回复: 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
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408 344分 求调剂 一志愿 华电计算机技术 +3 solanXXX 2026-03-20 3/150 2026-03-20 16:41 by fxue1114
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 271材料工程求调剂 +7 .6lL 2026-03-18 7/350 2026-03-20 09:10 by xingguangj
[考研] 0703化学调剂 ，六级已过，有科研经历 +12 曦熙兮 2026-03-15 12/600 2026-03-19 19:42 by maocaozhuxi
[考研] 化学求调剂 +3 临泽境llllll 2026-03-17 4/200 2026-03-19 13:59 by houyaoxu
[考研] 346求调剂[0856] +3 WayneLim327 2026-03-16 6/300 2026-03-19 11:21 by WayneLim327
[考研] 本科郑州大学物理学院，一志愿华科070200学硕，346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 化学工程321分求调剂 +15 大米饭！ 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 0854可跨调剂，一作一项核心论文五项专利，省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 326求调剂 +5 上岸的小葡 2026-03-15 6/300 2026-03-17 17:26 by ruiyingmiao
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 311求调剂 +3 冬十三 2026-03-13 3/150 2026-03-13 20:41 by JourneyLucky
信息提示
关闭
请填处理意见
关闭