恰當(dāng)方程
源代碼如下:
% 定義勢函數(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ù)等值線');
xlabel('x');
ylabel('y');
grid on;
一個是符號計(jì)算在 Python 插件出來后可以調(diào)用 SymPy, 另一個是 fcontour 函數(shù)是沒有的, 在北太天元中可以使用 contourf 函數(shù)配合 meshgrid 來實(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ù)等值線"); xlabel("x"); ylabel("y"); grid on;
結(jié)果如圖:
這里北太天元的坐標(biāo)刻度標(biāo)的有問題, 其他圖像和 MATLAB 的是一樣的.