我做了#1
,但在#2
我需要在代码中插入Aij = cos(...)
吗?
我有代码,但不知道我应该在哪里插入给定的条目Aij
。任何帮助将不胜感激!谢谢!
Code
n=8;
for i=1:n
for j=1:n
a(i,j)=1/(i+j-1);
end
end
c=ones(n,1);
b=a*c;
aorig=a;borig=b;
%Elimination
for i=1:n-1
for j=i+1:n
m=a(j,i)/a(i,i);
for k=i:n
a(j,k)=a(j,k)-m*a(i,k);
end
b(j)=b(j)-m*b(i);
end
end
%Back substitution
for i=n:-1:1
for j=i+1:n
b(i)=b(i)-a(i,j)*x(j);
end
x(i)=b(i)/a(i,i);
end
x'
答案 0 :(得分:1)
用#2中的余弦表达式替换希尔伯特矩阵定义a(i,j)=1/(i+j-1)
。