資源描述：
《matlab求解非線性方程》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、1、解方程最近有多人問(wèn)如何用matlab解方程組的問(wèn)題，其實(shí)在matlab中解方程組還是很方便的，例如，對(duì)于代數(shù)方程組Ax=b(A為系數(shù)矩陣，非奇異)的求解，MATLAB中有兩種方法：(1)x=inv(A)*b—采用求逆運(yùn)算解方程組；(2)x=A—采用左除運(yùn)算解方程組。例:x1+2x2=82x1+3x2=13>>A=[1,2;2,3];b=[8;13];>>x=inv(A)*bx=2.003.00>>x=Ax=2.003.00；即二元一次方程組的解x1和x2分別是2和3。對(duì)于同學(xué)問(wèn)到的用matlab解多次的方程組，有符號(hào)
2、解法，方法是：先解出符號(hào)解,然后用vpa(F,n)求出n位有效數(shù)字的數(shù)值解.具體步驟如下：第一步:定義變量symsxyz...;第二步:求解[x,y,z,...]=solve('eqn1','eqn2',...,'eqnN','var1','var2',...'varN');第三步:求出n位有效數(shù)字的數(shù)值解x=vpa(x,n);y=vpa(y,n);z=vpa(z,n);...。如：解二（多）元二（高）次方程組：x^2+3*y+1=0y^2+4*x+1=0解法如下：>>symsxy;>>[x,y]=solve('x^2+
3、3*y+1=0','y^2+4*x+1=0');>>x=vpa(x,4);>>y=vpa(y,4);結(jié)果是：x=1.635+3.029*i1.635-3.029*i-.283-2.987y=1.834-3.301*i1.834+3.301*i-.3600-3.307。二元二次方程組，共4個(gè)實(shí)數(shù)根；還有的同學(xué)問(wèn)，如何用matlab解高次方程組（非符號(hào)方程組）？舉個(gè)例子好嗎？解答如下：基本方法是：solve(s1,s2,…,sn,v1,v2,…,vn)，即求表達(dá)式s1,s2,…,sn組成的方程組，求解變量分別v1,v2,…,
4、vn。具體例子如下：x^2+x*y+y=3x^2-4*x+3=0解法：>>[x,y]=solve('x^2+x*y+y=3','x^2-4*x+3=0')運(yùn)行結(jié)果為x=13y=1-3/2即x等于1和3；y等于1和-1.5或>>[x,y]=solve('x^2+x*y+y=3','x^2-4*x+3=0','x','y')x=13y=1-3/2結(jié)果一樣，二元二方程都是4個(gè)實(shí)根。通過(guò)這三個(gè)例子可以看出，用matlab解各類(lèi)方程組都是可以的，方法也有多種，只是用到解方程組的函數(shù)，注意正確書(shū)寫(xiě)參數(shù)就可以了，非常方便。citefr
5、om:http://bbs.seu.edu.cn/pc/pccon.php?id=950&nid=14498&tid=02、變參數(shù)非線性方程組的求解對(duì)于求解非線性方程組一般用fsolve命令就可以了，但是對(duì)于方程組中某一系數(shù)是變化的，該怎么求呢？%定義方程組如下，其中k為變量functionF=myfun(x,k)H=0.32;Pc0=0.23;W=0.18;F=[Pc0+H*(1+1.5*(x(1)/W-1)-0.5*(x(1)/W-1)^3)-x(2);x(1)-k*sqrt(x(2))];%求解過(guò)程H=0.32;P
6、c0=0.23;W=0.18;x0=[2*W;Pc0+2*H];%取初值options=optimset('Display','off');k=0:0.01:1;%變量取值范圍[01]fori=1:1:length(k)kk=k(i);x=fsolve(@(x)myfun(x,kk),x0,options);%求解非線性方程組x1(i)=x(1);x2(i)=x(2);endplot(k,x1,'-b',k,x2,'-r');xlabel('k')legend('x1','x2')citefrom:http://foru
7、m.simwe.com/archiver/tid-836299.html3、非線性方程數(shù)值求解matlab里solve如何使用,是否有別的函數(shù)可以代替它.matlab里我解y=9/17*exp(-1/2*t)*17^(1/2)*sin(1/2*17^(1/2)*t)=0這樣的方程為什么只得到0這一個(gè)解,如何可以的到1/2*17^(1/2)*t=n*(pi)這樣一族解??在matlab里面solve命令主要是用來(lái)求解代數(shù)方程（即多項(xiàng)式）的解,但是也不是說(shuō)其它方程一個(gè)也不能解，不過(guò)求解非代數(shù)方程的能力相當(dāng)有限，通常只能給出很
8、特殊的實(shí)數(shù)解。(該問(wèn)題給出的方程就是典型的超越方程，非代數(shù)方程)從計(jì)算機(jī)的編程實(shí)現(xiàn)角度講，如今的任何算法都無(wú)法準(zhǔn)確的給出任意非代數(shù)方程的所有解，但是我們有很多成熟的算法來(lái)實(shí)現(xiàn)求解在某點(diǎn)附近的解。matlab也不例外，它也只能給出任意非代數(shù)方程在某點(diǎn)附近的解，函數(shù)有兩個(gè)：fzero和fsolve,具體用法請(qǐng)用help或
