Symbolic Math Toolbox™ обеспечивает аналитическое графическое изображение математических выражений, не генерируя явным образом числовые данные. Эти графики могут быть 2-D или 3-D как линии, кривые, контуры, поверхности или сети.
Эти примеры имеют следующие графические функции, которые принимают символьные функции, выражения и уравнения в качестве входов:
fplot
fimplicit
fcontour
fplot3
fsurf
fmesh
fimplicit3
fplot
Постройте график функции .
syms x
fplot(sin(exp(x)))
Постройте график тригонометрических функций , , и одновременно.
fplot([sin(x),cos(x),tan(x)])
Постройте график функции для , и .
syms x a expr = sin(exp(x/a)); fplot(subs(expr,a,[1,2,4])) legend show
Постройте график функции , его производная , и его интеграл .
syms f(x)
f(x) = x*(1 + x) + 2
f(x) =
f_diff = diff(f(x),x)
f_diff =
f_int = int(f(x),x)
f_int =
fplot([f,f_diff,f_int]) legend({'$f(x)$','$df(x)/dx$','$\int f(x)dx$'},'Interpreter','latex','FontSize',12)
Найти который минимизирует функцию путем решения дифференциального уравнения .
syms g(x,a);
assume(a>0);
g(x,a) = a*x*(a + x) + 2*sqrt(a)
g(x, a) =
x0 = solve(diff(g,x),x)
x0 =
Постройте график минимального значения для от 0 до 5.
fplot(g(x0,a),[0 5]) xlabel('a') title('Minimum Value of $g(x_0,a)$ Depending on $a$','interpreter','latex')
fimplicit
Постройте графики, заданные как с радиусом как целые числа от 1 до 10.
syms x y r = 1:10; fimplicit(x^2 + y^2 == r.^2,[-10 10]) axis square;
fcontour
Графическое изображение контуров функции для уровней контура от -6 до 6.
syms x y f(x,y) f(x,y) = x^3 - 4*x - y^2; fcontour(f,[-3 3 -4 4],'LevelList',-6:6); colorbar title 'Contour of Some Elliptic Curves'
Постройте график аналитической функции .
syms f(x)
f(x) = x*exp(-x)*sin(5*x) -2;
fplot(f,[0,3])
Создайте несколько точек данных из аналитической функции.
xs = 0:1/3:3; ys = double(subs(f,xs));
Постройте график точек данных и сплайн интерполяции, который аппроксимирует аналитическую функцию.
hold on plot(xs,ys,'*k','DisplayName','Data Points') fplot(@(x) spline(xs,ys,x),[0 3],'DisplayName','Spline interpolant') grid on legend show hold off
Найти разложение Тейлора рядом до 5-го и 7-го порядков.
syms x t5 = taylor(cos(x),x,'Order',5)
t5 =
t7 = taylor(cos(x),x,'Order',7)
t7 =
График и его приближения Тейлора.
fplot(cos(x)) hold on; fplot([t5 t7],'--') axis([-4 4 -1.5 1.5]) title('Taylor Series Approximations of cos(x) up to 5th and 7th Order') legend show hold off;
Квадратная волна периода и амплитуды может быть аппроксимировано расширением ряда Фурье
Постройте график квадратной волны с точкой и амплитуды .
syms t y(t) y(t) = piecewise(0 < mod(t,2*pi) <= pi, pi/4, pi < mod(t,2*pi) <= 2*pi, -pi/4); fplot(y)
Постройте график приближения квадратной волны в рядах Фурье.
hold on; n = 6; yFourier = cumsum(sin((1:2:2*n-1)*t)./(1:2:2*n-1)); fplot(yFourier,'LineWidth',1) hold off
Приближение Фурье перерегулирования при разрыве перехода, и «звон» не вымирает, так как в приближение добавляется больше членов. Такое поведение также известно как явление Гиббса.
fplot3
Постройте график спирали, которая определяется для от -10 до 10.
syms t fplot3(sin(t),cos(t),t/4,[-10 10],'LineWidth',2) view([-45 45])
fsurf
Постройте график поверхности, заданный как . Аналитическое графическое изображение с использованием fsurf
(без генерации числовых данных) показывает изогнутые области и асимптотические области около .
syms x y fsurf(log(x) + exp(y),[0 2 -1 3]) xlabel('x')
fsurf
Постройте график многомерной поверхности, заданный как
где .
Установите интервал графика от -5 до 5 и от 0 до 2.
syms f(u) x(u,v) y(u,v) z(u,v) f(u) = sin(u)*exp(-u^2/3)+1.5; x(u,v) = u; y(u,v) = f(u)*sin(v); z(u,v) = f(u)*cos(v); fsurf(x,y,z,[-5 5 0 2*pi])
fmesh
Постройте график многомерной поверхности, заданный как
где . Отобразите нанесенную на график поверхность как сетки при помощи fmesh
. Установите интервал графика от 0 до 2 и от 0 до .
syms s t r = 8 + sin(7*s + 5*t); x = r*cos(s)*sin(t); y = r*sin(s)*sin(t); z = r*cos(t); fmesh(x,y,z,[0 2*pi 0 pi],'Linewidth',2) axis equal
fimplicit3
Постройте график неявной поверхности .
syms x y z f = 1/x^2 - 1/y^2 + 1/z^2; fimplicit3(f)
Постройте график поверхности использование fsurf
. Можно показать контуры на том же графике путем установки 'ShowContours'
на 'on'
.
syms x y f = sin(x)+sin(y)-(x^2+y^2)/20
f =
fsurf(f,'ShowContours','on') view(-19,56)
Далее постройте график контуров на отдельном графике с более мелкими контурными линиями.
fcontour(f,[-5 5 -5 5],'LevelStep',0.1,'Fill','on') colorbar
Найдите градиент поверхности. Создайте 2-D сетки с помощью meshgrid
и замените координаты сетки, чтобы вычислить градиент численно. Отобразите градиент с помощью quiver
.
hold on
Fgrad = gradient(f,[x,y])
Fgrad =
[xgrid,ygrid] = meshgrid(-5:5,-5:5); Fx = subs(Fgrad(1),{x,y},{xgrid,ygrid}); Fy = subs(Fgrad(2),{x,y},{xgrid,ygrid}); quiver(xgrid,ygrid,Fx,Fy,'k') hold off