求助高手看下程序怎么修改(对隐函数画图)
本帖最后由 cypro 于 2011-1-15 08:04 编辑a0=2,deltae=0.01,ro0/a0=0.8;
c0=0.01,epson0=5.8e7,epsone=10*epson0;
mu0=pi*4e7,mue=mu0,mus=mu0;
deltas=1e5,eta0=120*pi;
for f=0:1e7:1e9 omiga=2*pi*f;
gama=alfa+j*beta;
gamae=j*mue*omiga*(deltae+j*epsone*omiga);
gama0=sqrt(-omiga^2*mu0*epson0);
zs=(sqrt(mus*omiga/(2*deltas)))*(1+j)/(2*pi*c0);
v=sqrt(gama0^2-gama^2);
u=sqrt(gamae^2-gama^2);
eps=0;
for m=-20:20 %这个小的for循环是对eps进行累加
zm=-j*mue*omiga*diff('besselk(m,u*a0)',u)/(u*besselk(m,u*a0));
ym=j*gamae^2*diff('besselk(m,u*a0)',u)/(u*besselk(m,u*a0));
delta0=(j*m*gama/a0)^2*(1/v^2-1/u^2)^2/(((gama0/v)*diff('besseli(m,v*a0)',u)/besseli(m,v*a0))+(zm/eta0));
rm=(((gama0/v)*diff('besselk(m,u*a0)',u)/besselk(m,u*a0))+ym+delta0)/(((gama0/v)*diff('besseli(m,v*a0)',u)/besseli(m,v*a0))+ym+delta0);
eps=eps+rm*besselk(m,v*a0)/besseli(m,v*a0)*besseli(m,v*ro0)*besseli(m,v*(c0+ro0));
end
-j*mu0*omiga*v^2/(2*pi*gama0)*(besselk(0,v*c0)-eps)=zs;
plot(f,alfa);
end
如上，倒数第三行是一个关于gama和omiga的方程，其他参数都由前述表示了。中间的小循环是对eps的求和。gama和omiga分别是f和alfa的函数，其他的函数关系如上书。现在要画f和alfa的函数关系图，感觉这样编程可能有点问题，本人才疏学浅，有同学能帮帮忙么。谢谢了。。。
页:
[1]