matlab如何算矩陣函式?如A2 0 0 1 1 1 1 1 3怎麼求exp A sin A

時間 2021-08-30 09:37:20

1樓:匿名使用者

用funm函式計算矩陣函式。

示例:>> a=[ 2 0 0; 1 1 1; 1 -1 3 ];

>> funm(a, @exp)

ans =

7.3891   -0.0000   -0.

00007.3891    0.0000    7.

38917.3891   -7.3891   14.

7781>> funm(a, @sin)

ans =

0.9093   -0.0000    0.

0000-0.4161    1.3254   -0.

4161-0.4161    0.4161    0.

4932

2樓:

直接輸入即可:

a=[ 2 0 0; 1 1 1; 1 -1 3 ]

a =2 0 0

1 1 1

1 -1 3

>> exp(a)

ans =

2431/329 1 1

1457/536 1457/536 1457/536

1457/536 536/1457 3053/152

>> sin(a)

ans =

401/441 0 0

1327/1577 1327/1577 1327/1577

1327/1577 -1327/1577 441/3125

你那是泰勒了,問題都不會問。

matlab程式設計 用乘冪法近似計算矩陣a的主特徵值和特徵向量,其中a=[2 -1 0;-1 2 0

3樓:匿名使用者

雅克賓方法

clc;

clear all;

%矩陣a

a=[2 ,-1,0;-1,2,-1;0,-1,2]%取矩陣a的維數

n=max(size(a));

%迭代誤差

eps=1e-5;

r=1;

%最大迭代次數為100

m=100;

k=1;

%小於迭代次數或迭代誤差進入計算

while r>=eps & k<=m

p=1;

q=1;

amax=0;

for i=1:n

for j=1:n

if i~=j & abs(a(i,j))>amaxamax=abs(a(i,j));

p=i;

q=j;

endend

endr=amax;%計算當前迭代誤差

%以下為構造正交矩陣u

l=-a(p,q);

u=(a(p,p)-a(q,q))/2;

if u==0

w=1;

else

w=sign(u)*l/sqrt(l*l+u*u);

ends=-w/sqrt(2*(1+sqrt(1-w*w)));

c=sqrt(1-s*s);

u=eye(n);

u(p,p)=c;

u(q,q)=c;

u(p,q)=-s;

u(q,p)=s;

%旋轉計算

a=u'*a*u%顯示每步計算a的計算結果k=k+1;

endif k>m

disp('a矩陣不收斂');

else

for i=1:n

d(i)=a(i,i);

enddisp('a特徵值為:');

dend

qr方法以a=[1 -1 2;-2 0 5;6 -3 6]為例不改了,自己改

構造矩陣

>>a=[1 -1 2;-2 0 5;6 -3 6]a =1 -1 2

-2 0 5

6 -3 6

將矩陣a變換為相似的擬上三角矩陣(即為上hessenberg矩陣)>>h=hess(a)

h =1.0000 2.2136 -0.

31626.3246 4.8000 -1.

40000 6.6000 1.2000對h矩陣作qr分解:

>>[q,r]=qr(h)

q =-0.1562 0.2101 -0.

9651-0.9877 -0.0332 0.

15260 0.9771 0.2127r =-6.

4031 -5.0868 1.43220 6.

7546 1.15260 0 0.3468作50次迭代計算(具體迭代次數可依具體實驗矩陣進行)>>for i=1:

50b=r*q;

[q,r]=qr(b);

end>>r*q

ans =

5.0000 7.4864 0.

5929-0.0000 3.0000 4.

96000 0.0000 -1.0000由以上結果可得到迭代計算的特徵值為,可見基本qr法的迭代精度還是很高的.

matlab計算c=d\c d=[0;1] a=[1,1;1,1]

4樓:我行我素

你上面說的算式有問題,

>> d=[0;1]; a=[1,1;1,1];

c=a\d

警告: 矩陣為奇異工作專精度

屬。 c =

-inf

inf>> d=[5;2]; a=[10,18;21,10];c=a\d

c =-0.0504

0.3058

matlab中compan函式怎麼具體計算??a=[1,2,3];compan(a)怎麼具體計算的呢?

5樓:匿名使用者

這個是求解伴隨矩陣或者叫友矩陣。

具體的方法是

a = compan(u) %u為多項式系統向量,a為友矩陣,a的第1行元素為 -u (2:n)/u(1),其中u (2:n)為u的第2到第n個元素,a為特徵值就是多項式的特徵根。

6樓:匿名使用者

你輸入:

doc compan

這條命令。裡面有說,怎麼生成的。

7樓:匿名使用者

可以在命令視窗裡面輸入help compan,matlab就給你英文的解釋和定義

matlab中已知矩陣為a=[1 2;3 4] 怎麼求某一行的和,b(1,1)=sum (a(1,))

8樓:匿名使用者

試試 b(1,1)=sum (a(1,:))

9樓:匿名使用者

b(1,1)=sum (a(1,:))

可以實現,a(1,:)表示一行所有列的元素相加a=[1 2;3 4];

>> b(1,1)=sum (a(1,:));

>> b

b =3

10樓:

a=[1 2;3 4];

b(1,1)=sum (a(1,:)) % 少冒號

《matlab》怎麼讓二值影象中0變成1,1變成0?

11樓:匿名使用者

很簡單,matlab專門提供復了一種取反操作,制就是符bai號“~”,就是數字鍵1左邊du的那個鍵的上檔鍵,zhi

比如daoa=[1,1,1;0,0,0];

那麼b=~a的結果就是b=[0,0,0;1,1,1];

如果是二值影象,也這樣做就行了。

matlab求和函式用什麼函式如 n 2 7n 3 n

柔情西瓜啊 求和是可以用sum 函式,得先用迴圈把每個元素計算出來 s 0 for i 1 100 s s i 2 7 i 3 或者 syms k symsum k.2 7 k 3,1,100 就可以了。擴充套件資料 函式定義 1 傳統定義 一般的,在乙個變化過程中,假設有兩個變數x y,如果對於任...

如何用matlab求逆矩陣,如何用MATLAB求逆矩陣

白鯨 第一步 啟動matlab。第二步 輸入 clear 和 clc 清屏 第三步 根據你的需求設定一個矩陣。圖中示例設定為矩陣a 1 2 3 4 a 可以定義為你需要的任何字母 第四步 用 b inv a b 可以定義為你需要的其他字母,inv 裡的字母為你需要求逆的矩陣。第五步 驗證自己求解的逆...

MATLAB如何求相關係數,MATLAB如何求相關係數

小心corrcoef的結果和您預期的不一樣啊。不過按你的要求是這樣的 x abcd y efgh cc corrcoef x,y result cc 2 result是乙個 1,1 區間內的數,表示x y兩個向量的相關係數。 簡單的說就是用corrcoef函式來計算。這是求相關度的結果,對於一般的矩...