system of equations (solve exponentiel equations
2 次查看(过去 30 天)
显示 更早的评论
Hi,
I want to solve this system of equations
P= exp(-r.^2.*P+R.^2.*Q).*sum((r.^2.*P +R.^2.*Q).^m./factorial(m)) // with m=(0:1:k-1); ;k=5
Q=exp(-r.^2.*P + R.^2.*Q.*sum((r.^2*P + R.^2.*Q).^n./factorial(n)) // n=(k:1:K_avg); K_avg=10
I tried to solve this system with solve. the result of execution : error eqns.
do you have an idea ?
th = -80:-3:-106;
list1 = zeros(size(th));
list2 = zeros(size(th));
r=[ 15.8 20.8 27.5 36.3 47.6 63.9 83.1 109.6 115.2];
R=[ 8.12 10.9 14.4 19.2 25.6 33.1 43.5 75.6 80.2];
for i = 1:numel(th)
K_avg=17;
k=9;
m=(0:1:k-1);
n=(k:1:K_avg);
syms P Q
eqns= [exp(-r.^2.*P+R.^2.*Q).*sum((r.^2.*P +R.^2.*Q).^m./factorial(m))-P==0,exp(-r.^2.*P + R.^2.*Q.*sum((r.^2*P + R.^2.*Q).^n./factorial(n))-Q==0)];
vars = [P Q];
[sol1,sol2]=solve(eqns,vars);
list1(i) = sol1;
list2(i) = sol2;
end
0 个评论
采纳的回答
madhan ravi
2018-9-10
@Marwen you get error/s because the number of elements m,n and R,r are not the same in order to make the loop perform calculations all the arrays should have the same number of elements.
9 个评论
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Calculus 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!