为什么contourf fcontour给不同的图像相同的matlab代码吗?

16个视图(30天)

显示旧的评论

Athira T Das 约10小时前

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/1974529-why-contourf-and-fcontour-gives-different-images-for-the-same-matlab-code

回答: VBBV 约1小时前

我有代码为fcountour;

信谊x y

m = 1;

sgm = 1;

我们= 0.01;

λ= 532 * 10 ^ 9;

k = 2 *π/λ;

ω= 2;

z = 30;

越南河粉= 0.003;

为u = 1:长度(z)

我= 0;

γ=((1我。* k)。/ (2。* z (u))) +(1. /(越南河粉(u)。^ 2)) +(1. /(卧室。^ 2));

Gamma_star =潜艇(γ1我1我);

α= (Gamma_star) - 1. /((伽马。*越南河粉(u)。^ 4));

C =(1. /(4。*(λ^ 2)。* (z (u) ^ 4)))。*((π^ 2)。/(α。*γ)。*(1. /(2 * 1我。* sqrt(γ))^米);

fx =((ω^ 2)。/(4。*γ))+ (x(1我。*。*。*ω)。/ (2。* z (u)。*γ))(((k * x ^ 2。^ 2)。/ (4 ^ 2 * z (u)。*γ));

财政年度=((ω^ 2)。/(4。*γ))+ (y(1我。*。*。*ω)。/ (2。* z (u)。*γ))- ((k * y ^ 2。^ 2)。/ (4 ^ 2 * z (u)。*γ));

gx =((ω^ 2)。/(4。*γ))——(x(1我。*。*。*ω)。/ (2。* z (u)。*γ))(((k * x ^ 2。^ 2)。/ (4 ^ 2 * z (u)。*γ));

gy =((ω^ 2)。/(4。*γ))——(y(1我。*。*。*ω)。/ (2。* z (u)。*γ))(((k * y ^ 2。^ 2)。/ (4 ^ 2 * z (u)。*γ));

Bx_p =((1我。* k。* x)。/ (2。* z (u))) +((1我。* k。* x)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) +(ω)。/(2。*伽马。*越南河粉(u)。^ 2) + (Omega. / 2);

Bx_n =((1我。* k。* x)。/ (2。* z (u))) +((1我。* k。* x)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) +(ω)。/(2。*伽马。*越南河粉(u)。^ 2)——(Omega. / 2);

By_p =((1我。* k。* y)。/ (2。* z (u))) +((1我。* k。* y)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) +(ω)。/(2。*伽马。*越南河粉(u)。^ 2) + (Omega. / 2);

By_n =((1我。* k。* y)。/ (2。* z (u))) +((1我。* k。* y)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) +(ω)。/(2。*伽马。*越南河粉(u)。^ 2)——(Omega. / 2);

Kx_p =((1我。* k。* x)。/ (2。* z (u))) +((1我。* k。* x)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) -(ω)。/(2。*伽马。*越南河粉(u)。^ 2) + (Omega. / 2);

Kx_n =((1我。* k。* x)。/ (2。* z (u))) +((1我。* k。* x)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) -(ω)。/(2。*伽马。*越南河粉(u)。^ 2)——(Omega. / 2);

Ky_p =((1我。* k。* y)。/ (2。* z (u))) +((1我。* k。* y)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) -(ω)。/(2。*伽马。*越南河粉(u)。^ 2) + (Omega. / 2);

Ky_n =((1我。* k。* y)。/ (2。* z (u))) +((1我。* k。* y)。/ (2。* z (u)。*伽马。*越南河粉(u)。^ 2)) -(ω)。/(2。*伽马。*越南河粉(u)。^ 2)——(Omega. / 2);

X = 0;

为l = 0: m

L =(((我)。* sgm)。^(马丁)。* nchoosek (m, L);

sum_j = 0;

为j = 0: m

J =(((我)。* sgm)。^ (m J)。* nchoosek (m, J);

sum1 = 0.0;

为p = 0:1/2: l / 2

faktor1_P = (((1) ^ p)。* ! (l))。/(γ(p + 1)。* ! (l - (2。* p)))。* (((2。* 1) / (sqrt(γ)))。^ (l - (2。* p)))。* exp (fx)。* exp(年度);

为s1 = 0: l2 * p

faktor1_s1 = nchoosek ((l - (2。* p)), s1)。*((1. /(越南河粉(u)。^ 2)) ^ s1)。*((((1我。* k。* x)。/ (2。* z (u))) + (Omega. / 2))。^ (l - (2。* p) - s1))。*(1. /(((2 * 1我。* sqrt(α)))。^ (j + s1)));

为q = 0:1/2:(马丁)/ 2

faktor1_Q = (((1) ^ q)。* factorial(马丁)。/(伽马(q + 1)。* factorial(马丁- (2。* q)))。* (((2。* 1) / (sqrt(γ)))。^(马丁- (2。* q)));

sum_inner_1 = 0;

为s2 = 0: m-l-2 * q

E1 = (exp (((Bx_p) ^ 2)。/α)。* hermiteH (j + s1,((1我。* Bx_p)。/ (sqrt(α))))。* exp (((By_p) ^ 2)。/α)。* hermiteH (m j + s2,((1我。* By_p)。/ (sqrt(α)))))——(exp (((Bx_n) ^ 2)。/α)。* hermiteH (j + s1,((1我。* Bx_n)。/ (sqrt(α))))。* exp (((By_n) ^ 2)。/α)。* hermiteH (m j + s2,((1我。* By_n)。/ (sqrt(α)))));

sum_inner_1 = sum_inner_1 + nchoosek((马丁- (2。* q)), s2)。*((1. /(越南河粉(u)。^ 2)) ^ s2)。*((((1我。* k。* y)。/ (2。* z (u))) + (Omega. / 2)) ^。(马丁- (2。* q) s2))。*(1. /(((2 * 1我。* sqrt(α)))。^ (m j + s2))) * E1;

结束

sum1 = sum1 + faktor1_P。* * faktor1_Q faktor1_s1。* sum_inner_1;

结束

sum_j = sum_j + j . * (sum1);

结束

X = X + (l * sum_j);

结束

X = X (c *);

我= + abs (vpa (X));

图;

fcontour(真正的(我),(1 1 1)* 1 e 1,“填满”,“上”,“MeshDensity”,150);

colormap (“热”);

轴(“平等”);colorbar (“绿色”);

结束

然后我contourf阴谋;

图;

m = 1;

sgm = 1;

我们= 0.1;

λ= 532 * 10 ^ 9;

k = 2 *π/λ;

ω= 2;

z = 30

越南河粉= 0.003

越南河粉= 0.0030

(x, y) = meshgrid (1:0.06:1);