代码搜索:hessenberg
找到约 26 项符合「hessenberg」的源代码
代码结果 26
www.eeworm.com/read/424281/10473411
c hessenberg.c
/* linalg/hessenberg.c
*
* Copyright (C) 2006 Patrick Alken
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as pu
www.eeworm.com/read/478033/6718393
pas hessenberg.pas
www.eeworm.com/read/266227/11235946
m hessenberg.m
function A=hessenberg(A)
%将矩阵A变换为上Heissenberg矩阵
[m,m]=size(A);
for k=1:m-2
%反复应用Householder变换得到上Heissenberg矩阵
[v,b]=householder(A(k+1:m,k));
H1= eye(m-k)-b*v*v';
H2=eye(m);
for i=k+1:m
fo
www.eeworm.com/read/323831/13314181
h mx_hessenberg.h
// Copyright (C) 2003 Zbigniew Leyk (zbigniew.leyk@anu.edu.au)
// and David E. Stewart (david.stewart@anu.edu.au)
// and Ronan Collobert (collober@idiap.ch)
//
www.eeworm.com/read/266227/11235949
m eig_qr.m
function r=eig_qr(A)
%qr方法计算矩阵特征值
A=hessenberg(A);
%将矩阵转化为上Hessenberg矩阵
[q1,r1]=qr_decomposition(A);
%对上Hessenberg矩阵进行qr分解,直到得到的新矩阵的对角线元素趋近于其特征值
A1=r1*q1;
[q2,r2]=qr_decomposition(A1);
A2=r2*q
www.eeworm.com/read/382536/9022667
m hessen.m
function A=hessen(A)
%功能:用Householder变换化矩阵A为Hessenberg形.
%输入: n阶实方阵A.
%输出: A的Hessenberg形.
%调用函数: househ.m
[n,n]=size(A);
for k=1:(n-2)
x=A(k+1:n,k);
H=househ(x);
A(k+1:n,1:n)=H*A
www.eeworm.com/read/424063/10499886
m freqresp.m
function g=freqresp(a,b,c,d,iu,s)
%FREQRESP Low level frequency response function.
%
% G=FREQRESP(A,B,C,D,IU,S)
% G=FREQRESP(NUM,DEN,S)
% Clay M. Thompson 7-10-90
% Copyright (c) 1986-93 by th
www.eeworm.com/read/147096/12584030
m freqresp.m
function g=freqresp(a,b,c,d,iu,s)
%FREQRESP Low level frequency response function.
%
% G=FREQRESP(A,B,C,D,IU,S)
% G=FREQRESP(NUM,DEN,S)
% Clay M. Thompson 7-10-90
% Copyright (c) 1986-93 by th
www.eeworm.com/read/101557/15826532
m freqresp.m
function g=freqresp(a,b,c,d,iu,s)
%FREQRESP Low level frequency response function.
%
% G=FREQRESP(A,B,C,D,IU,S)
% G=FREQRESP(NUM,DEN,S)
% Clay M. Thompson 7-10-90
% Copyright (c) 1986-93 by th
www.eeworm.com/read/382536/9022661
m qrtran.m
function A=qrtran(m,A)
%功能:对矩阵A的左上角的m阶对角块作QR变换:先用Givens变换作QR分解A=QR,
%再作相似变换A:=Q'AQ=RQ.
%输入: n阶HessenbergA,其中A(m+1,m)=0,m>2.
%输出: 变换后的Hessenberg形矩阵A.
Q=diag(ones(1,m));
for i=1:m-1
xi=A(i,i)