I obtained different figure after using fsurf and surf
5 次查看(过去 30 天)
显示 更早的评论
Here is the code, I think these two commands will generate the same figure but it did not.
% syms x t;
k = 0.6;
sigma=-0.5;
r =0.5;
s = 1.2;
x=10;t=0;
lambda=k.*sqrt(2.*r.^2./-sigma);
A=lambda.*((k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2))./(jacobiDN(r.*x-s.*t,k.^2)));
AA=lambda.*((k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2))./(jacobiDN(r.*x-s.*t,k.^2)));
B=4.*r.*s.*ellipticE(k.^2).*ellipticK(k.^2).^(-1)+(-2).*r.*s.*jacobiDN(s.*t+(-1).*r.*x,k.^2).^(-2).*(1+(-1).*k.^2+jacobiDN(s.*t+(-1).*r.*x,k.^2).^4);
A;
B;
lambda;
[x,t]=meshgrid(-5:0.2:5);
bb=k.*jacobiSN(r.*x-s.*t,k.^2).*jacobiCN(r.*x-s.*t,k.^2)./(jacobiDN(r.*x-s.*t,k.^2));
figure;
mesh(x,t,abs(bb));
syms c d;
bb=k.*jacobiSN(r.*c-s.*d,d.^2).*jacobiCN(r.*c-s.*d,k.^2).*(jacobiDN(r.*c-s.*d,k.^2)).^(-1);
figure;
fsurf(abs(bb),'MeshDensity',10);
g=400-300+500;
g
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/287990/image.png)
1 个评论
Ameer Hamza
2020-4-27
In your code, the 2nd equation is reciprocal of the first. Shouldn't that change the output.
采纳的回答
Toder
2020-4-28
It appears you may have a typo in your second definition of bb. Your first call to jacobiSN in this line currently reads as
jacobiSN(r.*c-s.*d,d.^2)
Based on your first definition of bb, it looks like it should instead be
jacobiSN(r.*c-s.*d,k.^2)
2 个评论
Toder
2020-4-29
Sorry, I don't know much about Mathematica. My only recommendation is to check that the axes are the same; if they're different than the surf in Matlab the plot could look very different.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Graphics Object Programming 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!