?? givens.m
字號:
function [K,L]=givens(H,n,eig1)
J=eye(n,n);
c=zeros(n,n);
s=zeros(n,n);
for i=2:n-1
for j=1:n
if H(i,j)*H(i,i)+H(i+1,j)*H(i+1,i)~=0
c(i,i)=H(i,j)*eig1(i,1)/(H(i,j)*H(i,i)+H(i+1,j)*H(i+1,i));
s(i,j)=H(i+1,j)*eig1(i,1)/(H(i,j)*H(i,i)+H(i+1,j)*H(i+1,i));
J(i,i)=c(i,i);
J(i,j)=s(i,j);
J(i+1,i)=s(i,j);
J(i+1,j)=c(i,i);
end
end
end
K=J*H
L=eig(J)
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -