Solving systems of equations graphically and finding where they cross.

6 次查看（过去 30 天）
Sumer
Sumer 2023-10-21
评论： Sumer 2023-10-21
Ran in:
I've got the following non-linear equation that I need to solve: F(x) = 2x - 3sin(x) +5. I know how to graph it and my code is below. However, I don't know how to find the intersection points. Am I supposed to use fplot or smth?
How can I find the exact coordinates for where y and F(x) intersect?
Could I have found the x-intercepts of F(x) without considering y?
%define x, y and F(x)
x = linspace(-20,20,2000);
F = 2*x-3*sin(x)+5;
y = zeros(size(x));
% plot the graphs on the same graph
hold on
plot(x,F)
plot(x,y)
% make the graph lookd neater
grid on
xlabel('x')
ylabel('F(x)')
title(' plot of F(x)')

请先登录，再回答此问题。

采纳的回答

Matt J
Matt J 2023-10-21
编辑：Matt J 2023-10-21
Ran in:
x_intersect=fzero(@(x) 2*x-3*sin(x)+5 ,[-20,+20])
x_intersect = -2.8832

更多回答（2 个）

Fabio Freschi
Fabio Freschi 2023-10-21
编辑：Fabio Freschi 2023-10-21
Ran in:
@Matt J solution is the cleanest. I have understood you want to find the graphical intersection. If you want to do it yourself you can calculate the line between the two points before and after the crossing and find the zero of the line.
clear variables, close all
%define x, y and F(x)
x = linspace(-20,20,2000);
F = @(x)2*x-3*sin(x)+5; % <- anonymous function
y = zeros(size(x));
% plot the graphs on the same graph
hold on
plot(x,F(x))
plot(x,y)
% make the graph lookd neater
grid on
xlabel('x')
ylabel('F(x)')
title(' plot of F(x)')
% find where you are crossing the 0
idx = find(diff(sign(F(x))));
% find crossing of linear interpolation of these two points
x0 = x(idx)+F(x(idx))*(x(idx+1)-x(idx))/(F(x(idx+1))-F(x(idx)))
x0 = -2.8997
plot(x0,F(x0),'o');
Additional notes:
  • You should also check if by chance, one of your x is a solution, that is F(x) = 0.
  • If you have more intersections, idx is a vector and you can loop over the vector length to find all intersections
Walter Roberson
Walter Roberson 2023-10-21
How can I find the exact coordinates for where y and F(x) intersect?
You cannot. There is no known closed-form expression for the coordinates. The best you can do is get a numeric approximation. You could use fzero or fsolve for that, or vpasolve

类别

在 Help Center 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by