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如何輸入函式,matlab 如何輸入導數
數碼答疑 導數採用diff函式 diff s,v n 其中s為求的函式,v是變數,n為求導次數例如 diff t 6 t 6 720一階導數 diff s 二階導數 diff s,2 diff 函式是matlab的符號函式求導。diff 使用說明 diff s,v 對表示式s中指定符號變數v,計算s...
Matlab如何實現樹結構,Matlab 如何實現樹結構
shine戚七七 由於matlab函式只能傳值傳遞引數的限制,操作起來稍有點不同。建立初始化連結串列空間 l,maintain initlist n,a1 a2 a3 例子 現在需要建立一樹,每個節點包含資料域,左指標域,和右指標域,估計需要的空間是100 l,maintain initlist 1...
matlab如何把矩陣中不同的元素組成新陣列
b為資料來源矩陣,假設為一維 a ka 1 for kb 1 length b if kb滿足某條件或者b kb 元素滿足某條件 a ka b kb ka ka 1 endend 以上 中b也可以為2維的,完全一樣的 但結果a是一維的,刷選出來的有if條件決定 上述 效率可能不佳 人或b如果不是幾個...