解方程命令

2021-04-27 13:44:41 字數 3677 閱讀 3952

求極限

limit:

例:limit(f,x,a,left);對表示式f求極限,變數為x,從左邊趨近a。

inf:正無窮;

-inf:負無窮。

求解線形方程

solve,linsolve

例:a=[5 0 4 2;1 -1 2 1;4 1 2 0;1 1 1 1];

%矩陣的行之間用分號隔開,元素之間用逗號或空格

b=[3;1;1;0]

x=zeros(4,1);%建立乙個4元列向量

x=linsolve(a,b)

diff(fun,var,n):對表示式fun中的變數var求n階導數。

例如:f=sym('u(x,y)*v(x,y)');   %sym()用來定義乙個符號表示式

diff(f);  %matlab區分大小寫

pretty(ans) %pretty():用習慣書寫方式顯示變數;ans是答案表示式

非線性方程求解

fsolve(fun,x0,options)

其中fun為待解方程或方程組的檔名;

x0位求解方程的初始向量或矩陣;

option為設定命令引數

建立檔案fun.m:

function y=fun(x)

y=[x(1)-0.5*sin(x(1))-0.3*cos(x(2)), ...

x(2) - 0.5*cos(x(1))+0.3*sin(x(2))];

>>clear;x0=[0.1,0.1];fsolve(@fun,x0,optimset('fsolve'))

注:...為續行符

m檔案必須以function為檔案頭,呼叫符為@;檔名必須與定義的函式名相同;fsolve()主要求解複雜非線性方程和方程組,求解過程是乙個逼近過程

不定積分與定積分

不定積分:int(fun,var)

例:求∫sinaxsinbxsincxdx

syms a b c x

y=sin(a*x)*sin(b*x)*sin(c*x);

int(y,x);

pretty(ans)

定積分:int(fun,var,a,b)

其中a,b分別為上下限

foll

應用到的函式:

i=imread('影象檔名') :讀取影象資料,儲存在矩陣i中;

imshow(i) :顯示灰度影象i,其他用法見matlab幫助;

i2=im2double(i1) :把影象陣列i1轉換成double精度型別;

d=dctmtx(n) :二維離散余弦變換函式,返回n*n離散余弦變換矩陣。

乙個n*n的變換矩陣t被定義成:

tpq=1/sqrt(n)                                   

,當p=0,0<=q<=m-1;

tpq=sqrt(2/n)*cos[pi*(2q+1)*p/2n]    

,當1<=p<=m-1,0<=q<=m-1。

b=blkproc(a,[m n],fun,p1,p2...) :塊操作函式。

對影象a的每個不同的m*n塊應用fun函式,p1,p2等為fun函式引數。在影象邊緣用0來擴充套件;只有當fun的返回矩陣是m*n矩陣時,b和a的大小才相同。

figure:強制生成乙個新的個繪圖視窗;

空間曲線

例:求方程組的空間曲線

1。化為引數方程組

x=t

y=sqrt[t(1-t)]

z=sqrt[1-x^2-y^2]

2。 程式

t=0:0.001:1;

x=t;

y=sqrt(t.*(1-t));

z=sqrt(1-x.^2-y.^2);

plot3(x,y,z,z,-y,z)

其中.*為陣列相乘,.^為陣列乘方

空間曲面

mesh()函式語法:

mesh(z):

mesh(x,y,z,c):其中c是用來定義相應點顏色等屬性的陣列

例:求x^2+y^2=z的空間曲面

x=-4:4;

y=x;

[x,y]=meshgrid(x,y);%生成x,y座標

z=x.^2+y.^2;

mesh(x,y,z)

級數的求和與收斂

symsum(fun,var,a,b):其中fun是通項表示式,var為求和變數,a為求和起點,b為求和終點

例如:i為1/[n*(2n+1)]從1到正無窮的和,求i

syms n;

f1=1/(n*(2*n+1));

i=symsum(f1,n,1,inf)

計算結果為:

i =2-2*log(2)

矩陣的基本運算

行列式求值:det(a)

矩陣加減:+、-

矩陣相:*

矩陣左除:a/b     %相當於inv(a)*b

矩陣右除:a/b     %相當於a*inv(b)

矩陣的:^

矩陣轉置:'

矩陣求共軛(實部相同,虛部相反):conj(x)

矩陣求:inv(x)

matlab之基本繪圖函式

clear:清空記憶體中的變數;

figure:強制生成乙個新的個繪圖視窗;

symsx y t :宣告變數;

fplot(函式表示式,繪圖區間);

plot(橫座標向量,縱座標向量,顏色/線形等引數)

ezplot(函式表示式):簡單的fplot,easy fplot

axis([xmin xmax ymin ymax ...]):設定座標軸顯示範圍

基本解方程

ncl是一家專門從事計算器改良與公升級的實驗室,最近該實驗室收到了某公司所委託的乙個任務 需要在該公司某型號的計算器上加上解一元一次方程的功能。實驗室將這個任務交給了乙個剛進入的新手zl先生。為了很好的完成這個任務,zl先生首先研究了一些一元一次方程的例項 4 3x 8 6a 5 1 2 2a 5 ...

數學 解方程

簡而言之,本題任務就是解方程。共有兩個子任務。作為小學生,我們只會解一元一次方程,一元一次方程最終都可以化為 ax n 的形式。現在問 對於給定的 n 要使得 x 有正整數解,總共可以取多少個不同的 a 呢?作為中學生,我們只會解二元一次不定方程,二元一次不定方程最終都可以化為 ax by n 的形...

Mysql解方程 解方程(codevs 3732)

題目描述 已知多項式方程 a0 a1x a2x 2 anx n 0 求這個方程在 1,m 內的整數解 n 和m 均為正整數 輸入輸出格式 輸入格式 輸入檔名為equation in。輸入共n 2 行。第一行包含2 個整數n m 每兩個整數之間用乙個空格隔開。接下來的n 1 行每行包含乙個整數,依次為...