How do I solve complex Integro-Differential equations in MATLAB?
1 次查看(过去 30 天)
显示 更早的评论
I am somewhat new to solving these kinds of equations but I would like to obtain some graphs/solutions using Analytical method preferebly.
2 个评论
采纳的回答
Dimitris Kalogiros
2018-8-9
编辑:Dimitris Kalogiros
2018-8-10
Hi Aditya
In order to solve this kind of equations symbolically, you should have Symbolic Math Toolbox.
Provided that, I'm giving you a tutorial piece of code:
clear all;
clc; close all;
% define symbolic variables
syms r n E w k d
syms mySum leftTerm
% left term of equation
leftTerm=(r/E)*(w/k)^2
% define term of summation
f=((-1)^(n+1))*(factorial(3)/factorial(2*n+1))*(d^(2*n-2))*(k^(2*n-2))
%calculate sum
mySum=symsum(f , n, 1, inf);
mySum=rewrite(mySum, 'sin')
% substitute variables other than d
rightHand=subs(mySum , k , 1)
leftHand=subs(leftTerm, [r E w k], [1 2 1 1]);
%plot left and right part of equation
fplot(rightHand, [-10 10]); hold on;
fplot(leftHand, [-10 10]); zoom on; grid on;
% numerical solving
equ= leftHand== rightHand
solNeg=vpasolve(equ, d, [-5 0])
solPos=vpasolve(equ, d, [0 5])
I suggest to run this script within a "Live Script"
At the following picture, you can see the results:
0 个评论
更多回答(0 个)
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!