恰当方程
源代码如下：
% 定义势函数
syms x y;
u = x^2*y + x^3*sin(y) + y^3/3;
% 正确调用 fcontour：使用单向量参数指定范围 [xmin xmax ymin ymax]
fcontour(u, [0 2 0 2]);
% 图形修饰
title('恰当方程的势函数等值线');
xlabel('x');
ylabel('y');
grid on;
一个是符号计算在 Python 插件出来后可以调用 SymPy, 另一个是 fcontour 函数是没有的, 在北太天元中可以使用 contourf 函数配合 meshgrid 来实现:
t = 0:0.1:2;
[X, Y] = meshgrid(t);
Z = X.^2 .* Y + X.^3 .* sin(Y) + Y.^3 / 3;
contourf(Z);
title("恰当方程的势函数等值线");
xlabel("x");
ylabel("y");
grid on;结果如图:
这里北太天元的坐标刻度标的有问题, 其他图像和 MATLAB 的是一样的.
公众号