请问这个程序出错在什么地方？
k_step = 1e-5;k = 0.0001 : k_step:0.1;
c = 1/2/(tanh(350*k).*(7320000*k.^4+49/5)./k).^(1/2)*((350-350*tanh(350*k).^2).*(7320000*k.^4+49/5)./k+29280000*tanh(350*k).*k.^2-tanh(350*k).*(7320000*k.^4+49/5)./k.^2);
semilogx(k,c) k_step = 1e-5; k = 0.0001 : k_step:0.1;
c1 = (tanh(350*k).*(7320000*k.^4+49/5)./k).^(1/2);
c2 = (350-350*tanh(350*k).^2).*(7320000*k.^4+49/5)./k;
c3 = 29280000*tanh(350*k).*k.^2-tanh(350*k).*(7320000*k.^4+49/5)./k.^2;
c = 1/2./c1.*(c2+c3);
semilogx(k,c)
程序出错在什么地方？
g = 9.8; H = 350; D = 7.32e6; k_step = 1e-5; k = 0.0001 : k_step:0.1;c=1/2./(k*tanh(k*H).*(D*k.^4+g)).^(1/2).*(tanh(k*H).*(D*k.^4+g)+k*(1-tanh(k*H).^2).*H*(D*k.^4+g).+4*k.^4.*tanh(k*H).*D).
semilogx(k,c)
[ 本帖最后由 ChaChing 于 2010-7-7 20:54 编辑 ] 个人打字慢, 所以2F就懒得说明! LZ可否发现我使用了c1/c2/c3等过渡参数?
个人以为式子写成那么长, 不说别人不想看, 自己也不好检查!
自己先试试看吧! 相信LZ可找出问题 谢谢提醒 ！
页:
[1]