帮我看看matlab程序错误function [h,q]=multioderarnoldi(A)n=1000;p=10;%端口数syms q;phi0=ones(n);phik=zeros(n);B=ones(n,p);R=inv(phi0)*B;%由上式求R,n*pK=zeros(n,p);%n*p[q0,K]=qr(R); M=4850;%多项式系统阶数q=(n-1)*p;for k=1:n-1Ak=-i
来源：学生作业帮助网 编辑：六六作业网 时间：2024/12/24 11:48:53
帮我看看matlab程序错误function [h,q]=multioderarnoldi(A)n=1000;p=10;%端口数syms q;phi0=ones(n);phik=zeros(n);B=ones(n,p);R=inv(phi0)*B;%由上式求R,n*pK=zeros(n,p);%n*p[q0,K]=qr(R); M=4850;%多项式系统阶数q=(n-1)*p;for k=1:n-1Ak=-i
帮我看看matlab程序错误
function [h,q]=multioderarnoldi(A)
n=1000;
p=10;%端口数
syms q;
phi0=ones(n);
phik=zeros(n);
B=ones(n,p);
R=inv(phi0)*B;%由上式求R,n*p
K=zeros(n,p);%n*p
[q0,K]=qr(R);
M=4850;%多项式系统阶数
q=(n-1)*p;
for k=1:n-1
Ak=-inv(phi0)*phik;%算Ak,n*n
m=min(k,M);
for j=1:m
while k-j>=0
qk(:,:,0)=Aj*q(k-j);
end
end
for j=1:k
H=q(k-j)'*qk(:,:,j-1);
qk(:,:,j)=qk(:,:,j-1)-q(k-j)*H;
[qk,K]=qr(qk(:,:,k));%qk^(k)的qr分解
end
end
Q=q(1:k,1:k)
Warning:Matrix is singular to working precision.
> In multioderarnoldi at 9
Warning:Matrix is singular to working precision.
> In multioderarnoldi at 15
Undefined function or variable 'Aj'.
Error in ==> multioderarnoldi at 19
qk(:,:,0)=Aj*q(k-j);
帮我看看matlab程序错误function [h,q]=multioderarnoldi(A)n=1000;p=10;%端口数syms q;phi0=ones(n);phik=zeros(n);B=ones(n,p);R=inv(phi0)*B;%由上式求R,n*pK=zeros(n,p);%n*p[q0,K]=qr(R); M=4850;%多项式系统阶数q=(n-1)*p;for k=1:n-1Ak=-i
可以把错误提示贴出来吗?这样不好解决.