恰當(dāng)方程
源代碼如下:
% 定義勢(shì)函數(shù)
syms x y;
u = x^2*y + x^3*sin(y) + y^3/3;
% 正確調(diào)用 fcontour:使用單向量參數(shù)指定范圍 [xmin xmax ymin ymax]
fcontour(u, [0 2 0 2]);
% 圖形修飾
title('恰當(dāng)方程的勢(shì)函數(shù)等值線');
xlabel('x');
ylabel('y');
grid on;
一個(gè)是符號(hào)計(jì)算在 Python 插件出來(lái)后可以調(diào)用 SymPy, 另一個(gè)是 fcontour 函數(shù)是沒(méi)有的, 在北太天元中可以使用 contourf 函數(shù)配合 meshgrid 來(lái)實(shí)現(xiàn):
t = 0:0.1:2;
[X, Y] = meshgrid(t);
Z = X.^2 .* Y + X.^3 .* sin(Y) + Y.^3 / 3;
contourf(Z);
title("恰當(dāng)方程的勢(shì)函數(shù)等值線");
xlabel("x");
ylabel("y");
grid on;結(jié)果如圖:

這里北太天元的坐標(biāo)刻度標(biāo)的有問(wèn)題, 其他圖像和 MATLAB 的是一樣的.
