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函式來計算。這是求相關度的結果,對於一般的矩...