如何用matlab解方程
一般用matlab求解方程有兩種方法:1、用軟件自帶的solve()函數、fsolve()函數、roots()函數等;2、利用數值分析的方法(如二分法,牛頓法等)自行編程其數值函數。
所以,要根據具體的方程形式來選擇合適的求解函數。請給以把具體的方程貼出來,以便于幫助你。
怎么用matlab解一元二次方程
1、首先打開matlab2017版軟件,新建一個空白的文件:
2、然后在軟件中輸入代碼,首先定義一元二次方程的3個參數,然后定義方程的判別式,接著計算方程的兩個根,最后用disp函數輸出即可:
3、最后運行程序,輸入方程的3個系數,即可得到根的結果:
matlab里怎么解方程組?
可以用vpasolve求解。實現代碼:
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組
matlab解方程
材料/工具
MATLAB R2017b
解一元方程
1
解一元方程:首先打開MATLAB R2017b(如下圖)
2
在“command Window”下直接輸入需要解決的方程的公式如:x=solve('x^2+2*x+1=0','x'),該公式是用于求方程“x^2+2*x+1=0”的解(如下圖)
3
接著點擊回車鍵后即可查看該方程的解(如下圖)
END
解二元方程
1
解二元一次方程:同樣在“command window”下直接輸入需要求解方程的公式,因為二元方程有兩個未知變量,所以解二元方程的公式如下例:[x,y]=solve('5*x+8*y=10','3*x+14*y=12','x','y')(如下圖)
2
方程之間還需要用“逗號”隔開(如下圖)
3
點擊回車鍵即可查看方程的解(如下圖)
4
解二元多次方程:方法與求解二元一次一樣,都是輸入需要求解方程的公式,但是獲取方程的解時,有所不同,如求解方程x^2+y^2=1,x+y=1,輸入公式[x,y]=solve('x^2+y^2=1','x+y=1','x','y')(如下圖)
5
點擊回車鍵查看方程的解,這種非線性方程的解,一般有多個,上面公式的解也可以表示一個圓心在原點,半徑為1的圓,與直線“x+y=1”的交點(如下圖)
END
解多元方程
1
解多元方程:求多元方程的解,方法與求解其他方程類似,只需要添加相應的變量值就行,比如求解一個三元方程:x+y+z=10,x+y-2z=5,3x-y+z=2的解,同樣在“command Window”界面下輸入公式:[x,y,z]=solve('x+y+z=10','x+y-2*z=5','3*x-y+z=2','x','y','z')(如下圖)
2
然后點擊回車鍵即可查看方程的解(如下圖)
如何用matlab解方程組
matlab中解方程組還是很方便的,例如,對于代數方程組Ax=b(A為系數矩陣,非奇異)的求解,MATLAB中有兩種方法:
(1)x=inv(A)*b — 采用求逆運算解方程組;
(2)x=A\B — 采用左除運算解方程組
PS:使用左除的運算效率要比求逆矩陣的效率高很多~
例:
x1+2x2=8
2x1+3x2=13
>>A=[1,2;2,3];b=[8;13];
>>x=inv(A)*b
x =
2.00
3.00
>>x=A\B
x =
2.00
3.00;
即二元一次方程組的解x1和x2分別是2和3。
對于同學問到的用matlab解多次的方程組,有符號解法,方法是:先解出符號解,然后用vpa(F,n)求出n位有效數字的數值解.具體步驟如下:
第一步:定義變量syms x y z ;
第二步:求解[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=0
y^2+4*x+1=0
解法如下:
>>syms x y;
>>[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*i
1.635-3.029*i
-.283
-2.987
y =
1.834-3.301*i
1.834+3.301*i
-.3600
-3.307。
二元二次方程組,共4個實數根;
解答如下:
基本方法是:solve(s1,s2,…,sn,v1,v2,…,vn),即求表達式s1,s2,…,sn組成的方程組,求解變量分別v1,v2,…,vn。
具體例子如下:
x^2 + x*y + y = 3
x^2 - 4*x + 3 = 0
解法:
>> [x,y] = solve('x^2 + x*y + y = 3','x^2 - 4*x + 3 = 0')
運行結果為
x =
1 3
y =
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 =
1 3
y =
1 -3/2
結果一樣,二元二方程都是4個實根。
通過這三個例子可以看出,用matlab解各類方程組都是可以的,方法也有多種,只是用到解方程組的函數,注意正確書寫參數就可以了,非常方便。
2、變參數非線性方程組的求解
對于求解非線性方程組一般用fsolve命令就可以了,但是對于方程組中某一系數是變化的,該怎么求呢?
%定義方程組如下,其中k為變量
function F = 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))];
%求解過程
H=0.32;
Pc0=0.23;W=0.18;
x0 = [2*W; Pc0+2*H]; % 取初值
options = optimt('Display','off');
k=0:0.01:1; % 變量取值范圍[0 1]
for i=1:1:length(k)
kk=k(i);
x = fsolve(@(x) myfun(x,kk), x0, options);%求解非線性方程組
x1(i)=x(1);
x2(i)=x(2);
end
plot(k,x1,'-b',k,x2,'-r');
xlabel('k')
legend('x1','x2')
用matlab解矩陣方程
1、加減法的命令很簡單,直接用加或者減號就可以了。如:c=a+bd=a-b。
2、一般乘法:c=a*b,要求a的列數等于b的行數。如果a,b是一般的向量,如a=[1,2,3] b=[3,4,5]點積:dot(a,b), 叉積:cross(a,b)卷積:conv(a,b)。
3、x=a\b如果ax=b,則x=a\b是矩陣方程的解。x=b/a如果xa=b, 則x=b/a是矩陣方程的解。
4、轉置時,矩陣的第一行變成第一列,第二行變成第二列,。。。x=a。
5、求逆:要求矩陣為方陣。這在矩陣運算中很常用。x=inv(a)。這幾種方式都可以解矩陣方程。