How do I solve complex Integro-Differential equations in MATLAB?

2018 8 6
1 个回答
回答已采纳
更新时间：2018 8 10
3 次查看（30 天）

0 个投票

I am somewhat new to solving these kinds of equations but I would like to obtain some graphs/solutions using Analytical method preferebly.

2 个评论
显示 无 隐藏 无

David Goodmanson
David Goodmanson 2018-8-9
Hi Aditya,
what variable are you trying to solve this for? Kappa, perhaps?
Aditya Sreekumar
Aditya Sreekumar 2018-8-9
Kappa is a constant and the del is varrying btw. But I just need a code example to plot any little complex intergo-diff equations to clear my basics. You can give any example. :)

请先登录，再进行评论。

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

 采纳的回答

Dimitris Kalogiros
编辑：Dimitris Kalogiros 2018-8-10

1 个投票

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 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Calculus 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by