|
楼主 |
发表于 2010-8-10 14:26
|
显示全部楼层
我刚刚稍微改了一下,用你的方法运行了一下,发现还是不是我想要的结果!
nRow=3;nCol=4;
for n=1:nRow*nCol,s=[nRow,nCol];[i,j]=ind2sub(s,n);
magic_str=['M',int2str(n),'=i+j'];eval(magic_str);
end
========================================
因为我想要的结果是输出一个矩阵A,不是A1,A2....An的,现在我想实现这种功能
++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
j=1;i=1;
while f(j)<F(i)&j<19981
p(j)=4*F(i)/ERBN(i)-(0.0117*4*F(i)/ERBN(i))*(20-51)
if abs(f(j)-F(i))/F(i)<=1
g(j)=abs(f(j)-F(i))/F(i);
else g(j)=1;
end
j=j+1;
end
for k=j:19981
p(k)=4*F(i)/ERBN(i);
if abs(f(k)-F(i))/F(i)<=4
g(k)=abs(f(k)-F(i))/F(i);
else g(k)=4;
end
end
=================================================
上面是我编的一个程序,其中是i=1的时候,而i是要从1到372中情况,而我现在想运行程序后能否得出矩阵p(372x19981)和g(372x19981),现在p和g都是1行19981列,所以现在就碰到这样一个难题了!chaching:你能帮我看看嘛,可能我写的有些乱,不好意思!
===================================================================================
我想表达的意思,下面我用程序改了一下!
for i=1:372
j=1;
while f(j)<F(i)&j<19981
p(i,j)=4*F(i)/ERBN(i)-(0.0117*4*F(i)/ERBN(i))*(20-51)
if abs(f(j)-F(i))/F(i)<=1
g(i,j)=abs(f(j)-F(i))/F(i);
else g(i,j)=1;
end
j=j+1;
end
for k=j:19981
p(i,k)=4*F(i)/ERBN(i);
if abs(f(k)-F(i))/F(i)<=4
g(i,k)=abs(f(k)-F(i))/F(i);
else g(i,k)=4;
end
end
===========================================================
[ 本帖最后由 jinfengshi_2005 于 2010-8-10 15:08 编辑 ] |
|