Solving systems of equations graphically and finding where they cross.
4 views (last 30 days)
Show older comments
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)')
0 Comments
Accepted Answer
More Answers (2)
Fabio Freschi
on 21 Oct 2023
Edited: Fabio Freschi
on 21 Oct 2023
@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)))
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
See Also
Categories
Find more on Particle & Nuclear Physics in Help Center and File Exchange
Products
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!