I obtained different figure after using fsurf and surf
5 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
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 Kommentar
Ameer Hamza
am 27 Apr. 2020
In your code, the 2nd equation is reciprocal of the first. Shouldn't that change the output.
Akzeptierte Antwort
Toder
am 28 Apr. 2020
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 Kommentare
Toder
am 29 Apr. 2020
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.
Weitere Antworten (0)
Siehe auch
Kategorien
Mehr zu Graphics Object Programming finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!