function question_43
clear;clc
vf(2,2,1,1,1)
function x=vf(Q,y,z,wc,hc)
x=(((1-cosh(pi*y/hc)/cosh(pi*wc/(2*hc)))*cos(pi*z/(2*hc)))-1/27*((1-cosh(3*pi*y/hc)/cosh(3*pi*wc/(2*hc)))*cos(3*pi*z/(2*hc))))*((((1-2*hc/(pi*wc)*tanh(pi*wc/(2*hc))))+1/81*((1-2*hc/(3*pi*wc)*tanh(3*pi*wc/(2*hc)))))^(-1))*pi*Q/(2*100*20)
end
end

2楼: Originally posted by pdl9527 at 2017-04-18 15:49:38
我调用了一下并没有出错。


function question_43
clear;clc
vf(2,2,1,1,1)
function x=vf(Q,y,z,wc,hc)
x=(((1-cosh(pi*y/hc)/cosh(pi*wc/(2*hc)))*cos(pi*z/(2*hc)))-1/27*((1-cosh(3*pi*y/hc)/cosh(3*pi ...

3楼: Originally posted by lla1016 at 2017-04-20 14:03:31
>> vf(,1,1,100,20)

vf =

  Columns 1 through 13

    0.0014    0.0028    0.0043    0.0057    0.0071    0.0085    0.0099    0.0114    0.0128    0.0142    0.0284    0.0426    0.0568

   ...

function question_43
clear;clc
x=[0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20];
y=vf(x,1,1,100,20);
plot(x,y,'ro-')
function x=vf(Q,y,z,wc,hc)
x=(((1-cosh(pi*y/hc)/cosh(pi*wc/(2*hc)))*cos(pi*z/(2*hc)))-1/27*((1-cosh(3*pi*y/hc)/cosh(3*pi*wc/(2*hc)))*cos(3*pi*z/(2*hc))))*((((1-2*hc/(pi*wc)*tanh(pi*wc/(2*hc))))+1/81*((1-2*hc/(3*pi*wc)*tanh(3*pi*wc/(2*hc)))))^(-1))*pi*Q/(2*100*20);
end
end