shuidixn 发表于 2011-7-11 20:51

求助矩阵求值画图

本帖最后由 shuidixn 于 2011-7-11 20:53 编辑

d1=??
nh=100000;
bc=1.016/nh;
en=(1:nh)*bc;      
k=(2*en).^0.5;
T01=abs((0.5*cosk-0.5./*sink).*cosh(k.*d1);
for q=1:nh-1;
    if    T01(q+1)<=T01(q)&&T01(q+1)<=T01(q+2);
          eN(s1)=en(q+1);
    end
s=eN(1);
如上,d1 我想设为linspace(40.140,100)这100个数,对于d1的每个值,都有一个s,然后画plot(d1,s)
除了for 循环可以,请问还有没有其他的办法,那个循环实在是太慢了,因为后面的程序还很长

ChaChing 发表于 2011-7-12 00:36

本帖最后由 ChaChing 于 2011-7-12 00:43 编辑

个人水平专业有限, 总感觉LZ给的代码, 真的有点乱!?
页: [1]
查看完整版本: 求助矩阵求值画图