我有一个正弦函数y = asin(b(x+c)),x为变量;
矩阵A为m*n的矩阵。
取x = 1:m
若y为正,则在A矩阵的第x行前边补round(y)个0,同时在该行的最后边删除round(y)个元素,使该行元素个数保持不变。
若y为负,则在A矩阵的第x行后边补round(y)个0,同时在该行的最前边删除round(y)个元素,使该行元素个数保持不变。
若y为0,则保持不变。
请问除了循环然后对每一行进行修改,有没有更好的方法实现。
可以但没必要,比较麻烦,可以参考如下代码:
a=2;
b=3;
c=2;
m=20;
n=10;
A=repmat(1:n,[m,1]);
x=1:m;
y=round(a.*sin(b.*(x+c)));
nA=[zeros(m,a),A,zeros(m,a)];
Index=((a-y')+(1:n)-1).*m+(1:m)';
B=zeros(m,n);
B(1:m*n)=nA(Index(1:m*n));
disp(A)
disp(y)
disp(B)
A矩阵:
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
y向量(round后):
1 -1 1 -2 2 -2 2 -2 2 -2 2 -2 2 -2 1 -1 1 -1 0 0
B矩阵:
0 1 2 3 4 5 6 7 8 9
2 3 4 5 6 7 8 9 10 0
0 1 2 3 4 5 6 7 8 9
3 4 5 6 7 8 9 10 0 0
0 0 1 2 3 4 5 6 7 8
3 4 5 6 7 8 9 10 0 0
0 0 1 2 3 4 5 6 7 8
3 4 5 6 7 8 9 10 0 0
0 0 1 2 3 4 5 6 7 8
3 4 5 6 7 8 9 10 0 0
0 0 1 2 3 4 5 6 7 8
3 4 5 6 7 8 9 10 0 0
0 0 1 2 3 4 5 6 7 8
3 4 5 6 7 8 9 10 0 0
0 1 2 3 4 5 6 7 8 9
2 3 4 5 6 7 8 9 10 0
0 1 2 3 4 5 6 7 8 9
2 3 4 5 6 7 8 9 10 0
1 2 3 4 5 6 7 8 9 10
1 2 3 4 5 6 7 8 9 10
我也问个问题,首先,你怎么保证,前面有round(y)个补充,后面有round(y)个可以删除?同理,后面补充round(y)个,前面怎么保证有round(y)个可以删除呢?
其二,这种操作,第i步的结果严重依赖上一步即第i-1步的结果,不循环你怎么可能实现呢?