for lambda=1:0.1:2
syms x
qr=1.449*lambda.*(1-0.1416*lambda.^2).^3.0303;
lambda1=vpasolve(1.57744*x*(1-0.1667*x.^2).^2.5==qr)
end
運行結果
matlab解方程組lnx表示成log(x)而lgx表示成log10(x)1-exp(((log(y))/x^0.5)/(x-1))1、解方程最近有多人問如何用matlab解方程組的問題，其實在matlab中解方程組還是很方便的，例如，對於代數方程組Ax=b(A為系數矩陣，非奇異)的求解，MATLAB中有兩種方法：(1)x=inv(A)*b—采用求逆運算解方程組； (2)x=A\B—采用左除運算解方程組PS：使用左除的運算效率要比求逆矩陣的效率高很多~例:x1+2x2=82x1+3x2=13>>A=[1,2;2,3];b=[8;13];>>x=inv(A)*bx=2.003.00 >>x=A\Bx=2.003.00；即二元壹次方程組的解x1和x2分別是2和3。對於同學問到的用matlab解多次的方程組，有符號解法，方法是：先解出符號解,然後用vpa(F,n)求出n位有效數字的數值解.具體步驟如下：第壹步:定義變量symsxyz...;第二步:求解[x,y,z,...]=solve('eqn1','eqn2',...,'eqnN','var1','var2',...'varN');第三步:求出n位有效數字的數值解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*y+1=0','y^2+4*x+1=0');>>x=vpa(x,4);>>y=vpa(y,4);結果是：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個實數根；還有的同學問，如何用matlab解高次方程組（非符號方程組）？舉個例子好嗎？解答如下：基本方法是：solve(s1,s2,…,sn,v1,v2,…,vn)，即求表達式s1,s2,…,sn組