function[output]=tridiag(d,l,r)
A=zeros(3);
for i=1:5
A(i,i)=d(i);
for j=2:5
A(j,(i-1))=l(j);
for k=1:4
A(k,(i+5))=r(k);
end
end
end
A
end
在这部分上面写着,
for j=2:5
A(j,(i-1))=l(j);
我希望变量"i“从1~ 4执行,而不是从5执行。这就是为什么我将i-1放在那里,但很明显,它不能工作。
如何使嵌套的j-循环重复i=1:4而不是i=1:5?同样的问题也适用于k循环。
+update:所以我希望外部I-循环从1:5开始执行,但是对于内部j-循环,我想要i=1:4。
发布于 2017-03-05 21:47:46
对于你最初的问题:在内环中使用循环控制。
for i=1:5
A(i,i)=d(i);
for j=2:5
if i==5
continue
end
A(j,(i-1))=l(j);
for k=1:4
A(k,(i+5))=r(k);
end
end
end
我猜您希望创建一个三对角矩阵,使用d
作为主要对角线,l
作为左/下对角线,r
作为右/上对角线。由于索引错误,您的代码似乎无法工作。
d = 1:5;
l = 10:10:50;
r = 100:100:500;
N = length(d);
A=zeros(N);
for ii=1:N
A(ii,ii)= d (ii);
if ii~=1
A(ii,(ii-1)) = l(ii);
end
if ii~=N
A(ii,(ii+1)) = r(ii);
end
end
A
A =
1 100 0 0 0
20 2 200 0 0
0 30 3 300 0
0 0 40 4 400
0 0 0 50 5
发布于 2017-03-05 21:26:50
A(j,(i-1))
我从1开始。所以在第一次迭代中,下面给出如下内容:
A(j,(0))
由于MATLAB中数组的索引从1开始,这是绝对错误的。因此,将您的代码更改为类似这样的代码将有效:
for j=1:4
A(j,(i))=l(j);
发布于 2017-03-06 02:27:52
function[output]=tridiag(d,l,r)
A=zeros(3);
for i=1:5
A(i,i)=d(i);
j=2;
while(j<6 & i~=5)
A(j,i)=l(j);
j=j+1;
for k=1:4
A(k,(i+5))=r(k);
end
end
end
A
end
以下三行确保除i=5外,“i”的每个值都从'2‘到'5’之间。
j=2;
while(j<6 & i~=5)
j=j+1;
https://stackoverflow.com/questions/42618065
复制