[交流] 用MATLAB解非線性方程組，fsolve，各種問題

主函數(shù)
global G ;                %全局變量定義
clear all;
global G ;                %全局變量定義
global L ;
global B ;
global rz ;
global l ;
global f ;
global u;
global K;
global g;
global b;
global cy;
global cx;
global h;
global p;

G=213800;%N
L=4.51;% m
B=2.84;% m
rz=0.2626;%m
l=0.148;%m
h=1.123;%m
f=0.05;
u=0.9;
K=0.015;
g=10;
b=0.29;
cx=0;
cy=0.31;
rou=1:100;
[a,c]=size(rou);
for n=a:1:c;
  p=rou(1,n);
x0=[0.2 0.4 0.2];  % Make a starting guess at the solution
options=optimset('Display','final','TolFun',1e-2,'MaxFunEvals',1000);  % Turn off Displayclc
[x,Fval]=fsolve(@myfunn1,x0,options) ;
z(1,n)=x(1);
z(2,n)=x(2);
z(3,n)=x(3);
n=n+1;
end
plot(rou,z(1,

);

% F=myfun(x0);
%             fun(1)=x(1);
%             fun(2)=x(2);
%             fun(3)=x(3);

子函數(shù)
%子函數(shù)
function F=myfunn1(x)
global G ;
global L ;
global B ;
global l ;
global f ;
global u;
global K;
global g;
global b;
global cy;
global cx;
global h;
global p;
v=0.1;
N1=0.5*G-(h*G*v^2*sqrtm((p*B)^2-x(1).^2)/(B*g*(p*B)*(p*B))-cx*G/B); %R=(p*B)
N2=0.5*G+(h*G*v^2*sqrtm((p*B)^2-x(1).^2)/(B*g*(p*B)*(p*B))-cx*G/B);
  Rf1=N1*f;
  Rf2=N2*f;
a=6*G/(g*b*l^3)*(cy+h*v^2.*x(1)/(g*(p*B)*(p*B)));
Sszn1=-b.*quad(@myfun1,-L/2,L/2);
function z1 = myfun1(y1)
         fai1=(L/2+cy-x(1)-y1)./(sqrt((p.*B)^2-x(1).^2)-B/2+cx-x(2));
jX1=(sqrt((p.*B)^2-x(1).^2)-B/2+cx).*(cos(fai1)-1)-y1.*sin(fai1);
jY1=(sqrt((p.*B)^2-x(1).^2)-B/2+cx).*sin(fai1)-L/2+y1.*cos(fai1);
j1=sqrt(jX1.^2+jY1.^2);
   yx1=-y1./(sqrt(x(2).^2+y1.^2));
z1 = u.*(N1./(l*b)-y1.*a).*(1-exp(-j1/K)).*yx1;
end
%  plot(rou,Sszn1);
%  hold on

Sszw1=quad(@myfu1,-L/2,L/2);  %+quad(@myfun2,1.205,1.501)+quad(@myfun2,0.303,0.599)+quad(@myfun2,-0.599,0.303)+quad(@myfun2,-1.501,-1.205)+quad(@myfun2,-2.255,-1.255));
function z2=myfu1(y2)
fai2=(L/2+cy-x(1)-y2)./(sqrt((p.*B)^2-x(1).^2)+B/2+cx-x(2));%fai2=(L/2+cy-D-y2)/(sqrt(R^2-D^2)+B/2+cx-A1);
jX2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*(cos(fai2)-1)-y2.*sin(fai2);
jY2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*sin(fai2)-L/2+y2.*cos(fai2);
j2=sqrt(jX2.^2+jY2.^2);
yx2=-y2./(sqrt(x(3).^2+y2.^2));
z2=u.*(N2./(l*b)-y2.*a).*(1-exp(-j2/K)).*yx2;
end

%子函數(shù)中對(duì)應(yīng)函數(shù)3  內(nèi)縱向力
%first
Szn1=-b.*quad(@myf1,-L/2,L/2);%+quad(@myfun3,1.205,1.501)+quad(@myfun3,0.303,0.599)+quad(@myfun3,-0.599,0.303)+quad(@myfun3,-1.501,-1.205)+quad(@myfun31,-2.255,-1.255))  ;
function z3=myf1(y3)

   fai1=(L/2+cy-x(1)-y3)./(sqrt((p*B)^2-x(1).^2)-B/2+cx-x(2));%fai2=(L/2+cy-D-y2)/(sqrt(R^2-D^2)+B/2+cx-A1);
jX1=(sqrt((p*B)^2-x(1).^2)-B/2+cx).*(cos(fai1)-1)-y3.*sin(fai1);
jY1=(sqrt((p*B)^2-x(1).^2)-B/2+cx).*sin(fai1)-L/2+y3.*cos(fai1);
j1=sqrt(jX1.^2+jY1.^2);
zx1=x(2)./(sqrt(x(2).^2+y3.^2));
z3=u.*(N1./(l*b)-y3.*a).*(1-exp(-j1/K)).*zx1;
end

Szw1=-b.*quad(@my1,-L/2,L/2);%+quad(@myfun4,1.205,1.501)+quad(@myfun4,0.303,0.599)+quad(@myfun4,-0.599,0.303)+quad(@myfun4,-1.501,-1.205)+quad(@myfun4,-2.255,-1.255));
function z4=my1(y2)
fai2=(L/2+cy-x(1)-y2)./(sqrt((p.*B)^2-x(1).^2)+B/2+cx-x(2));%fai=(L/2+cy-x(1)-y1)./(sqrt((p*B)^2-x(1).^2)-B/2+cx-x(2));
jX2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*(cos(fai2)-1)-y2.*sin(fai2);
jY2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*sin(fai2)-L/2+y2.*cos(fai2);
j2=sqrt(jX2.^2+jY2.^2);
zx2=-x(3)./(sqrt(x(3).^2+y2.^2));
z4=u.*(N2./(l*b)-y2.*a).*(1-exp(-j2/K)).*zx2;
end

%內(nèi)驅(qū)動(dòng)力矩
MDzn=-0.5.*B.*Szn1;
%外驅(qū)動(dòng)力矩
MDzw=-0.5.*B.*Szw1;

%子函數(shù)中對(duì)應(yīng)函數(shù)5  內(nèi)阻力矩
%first
Muzn1=-b.*quad(@myju1,-L/2,L/2);%+quad(@myfun7,1.205,1.501)+quad(@myfun7,0.303,0.599)+quad(@myfun7,-0.599,0.303)+quad(@myfun7,-1.501,-1.205)+quad(@myfun7,-2.255,-1.255));
function z7=myju1(y1)
   fai1=(L/2+cy-x(1)-y1)./(sqrt((p*B)^2-x(1).^2)-B/2+cx-x(2));
jX1=(sqrt((p*B)^2-x(1).^2)-B/2+cx).*(cos(fai1)-1)-y1.*sin(fai1);
jY1=(sqrt((p*B)^2-x(1).^2)-B/2+cx).*sin(fai1)-L/2+y1.*cos(fai1);
j1=sqrt(jX1.^2+jY1.^2);
yx1=-y1./(sqrt(x(2).^2+y1.^2));
z7=y1*u.*(N1./(l*b)-y1.*a).*(1-exp(-j1/K)).*yx1;
end

%子函數(shù)中對(duì)應(yīng)函數(shù)6  外阻力矩
%first
Muzw1=-b.*quad(@mytes1,-L/2,L/2);%+quad(@myfun8,1.205,1.501)+quad(@myfun8,0.303,0.599)+quad(@myfun8,-0.599,0.303)+quad(@myfun8,-1.501,-1.205)+quad(@myfun8,-2.255,-1.255));
  function z8=mytes1(y2)
fai2=(L/2+cy-x(1)-y2)./(sqrt((p.*B)^2-x(1).^2)+B/2+cx-x(2));
jX2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*(cos(fai2)-1)-y2.*sin(fai2);
jY2=(sqrt((p*B)^2-x(1).^2)+B/2+cx).*sin(fai2)-L/2+y2.*cos(fai2);
j2=sqrt(jX2.^2+jY2.^2);
yx2=-y2./(sqrt(x(3).^2+y2.^2));
z8=y2*u.*(N2./(l*b)-y2.*a).*(1-exp(-j2/K)).*yx2;
  end


F(1)=Sszn1+Sszw1-G*v.^2.*sqrt((p*B)^2-x(1).^2)./(g*p*B*p*B);
F(2)=Szn1+Szw1-G*v.^2.*x(1)./(g*p*B*p*B)-Rf1-Rf2;
F(3)=MDzw-MDzn-0.5.*B.*(Rf2-Rf1)+(x(1).*sqrt((p*B).^2-x(1).^2)./(p*B)+cx.*x(1)./(p*B)).*G.*v.^2./(g.*p*B)-Muzw1+Muzn1;
% end
end

回復(fù)此樓

» 猜你喜歡

0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有17人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)
307求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求助matlab高手指點(diǎn)：大型非線性方程組（牛頓拉佛森法）已經(jīng)有7人回復(fù)
matlab求解非線性方程組報(bào)錯(cuò)，請(qǐng)各位大神指點(diǎn)！已經(jīng)有8人回復(fù)
Matlab用fsolve求解非線性方程，卡了三四天了跪求幫助！已經(jīng)有35人回復(fù)
MATLAB解非線性方程組已經(jīng)有9人回復(fù)
matlab解非線性方程組程序糾錯(cuò)與正確運(yùn)行已經(jīng)有6人回復(fù)
親們，求助一下，這個(gè)非線性方程組怎么解？已經(jīng)有11人回復(fù)
matlab解方程組求助攻已經(jīng)有4人回復(fù)
matlab 求解線性方程組Ax=b 已經(jīng)有10人回復(fù)
求教如何用matlab解帶約束的非線性方程組已經(jīng)有4人回復(fù)
matlab 非線性微分方程求解已經(jīng)有3人回復(fù)
用matlab求解一個(gè)非線性方程組的解已經(jīng)有4人回復(fù)
請(qǐng)教一個(gè)matlab求解非線性方程組的問題已經(jīng)有9人回復(fù)
matlab 求解非線性方程組未果已經(jīng)有10人回復(fù)
matlab fsolve得出的解全部為復(fù)數(shù)解已經(jīng)有15人回復(fù)
用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？已經(jīng)有11人回復(fù)
MATLAB求解非線性方程組已經(jīng)有5人回復(fù)
求Matlab解非線性方程代碼已經(jīng)有12人回復(fù)
Matlab：Newton-Raphson 解線非線性性方程組的問題已經(jīng)有12人回復(fù)
使用fsolve求解非線性方程問題已經(jīng)有9人回復(fù)
matlab求解非線性方程組，錯(cuò)誤提示怎么解決已經(jīng)有5人回復(fù)
求助matlab---fsolve解非線性方程組已經(jīng)有6人回復(fù)
求助！matlab用fsolve函數(shù)求解非線性方程組的問題！已經(jīng)有19人回復(fù)
用matlab求解方程出問題，請(qǐng)幫忙看看已經(jīng)有3人回復(fù)
matlab求解非線性方程組已經(jīng)有16人回復(fù)
求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
matlab的fsove 命令求解非線性方程組已經(jīng)有6人回復(fù)
【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(shù)解！已經(jīng)有18人回復(fù)

1樓 2014-10-15 23:36:12

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

will_zhan

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
帖子: 2
在線: 14分鐘
蟲號(hào): 3464364
注冊(cè): 2014-10-10

結(jié)構(gòu)很清晰，就是一步一步表示出表達(dá)式，然后帶入后面的式子中，最后對(duì)求出的表達(dá)式積分得到要帶入最后非線性方程組的式子，麻煩就麻煩再積分的地方，由于新手也可能出現(xiàn)其他一些始料不及的小錯(cuò)誤，多指教

贊一下

回復(fù)此樓

2樓2014-10-15 23:38:57

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +8	薛云鵬 2026-03-13	8/400	2026-03-19 15:36 by haoshis
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +11	yangfz 2026-03-17	11/550	2026-03-19 15:06 by houyaoxu
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7	小材化本科 2026-03-18	7/350	2026-03-19 10:46 by Linda Hu
[考研] 一志愿天大材料與化工（085600）總分338 +5	蔡大美女 2026-03-13	5/250	2026-03-19 10:44 by 是小劉呀～
[考研] 材料專碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 一志愿華中科技大學(xué)，080502，354分求調(diào)劑 +4	守候夕陽CF 2026-03-18	4/200	2026-03-18 22:16 by li123456789.
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3	xbxudjdn 2026-03-18	3/150	2026-03-18 22:14 by zhq0425
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 301求調(diào)劑 +9	yy要上岸呀 2026-03-17	9/450	2026-03-18 08:58 by 無際的草原
[考研] 334求調(diào)劑 +3	志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 293求調(diào)劑 +11	zjl的號(hào) 2026-03-16	16/800	2026-03-18 08:10 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 0856求調(diào)劑 +3	劉夢(mèng)微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 330求調(diào)劑 +3	?醬給調(diào)劑跪了 2026-03-13	3/150	2026-03-14 10:13 by JourneyLucky