How to calculation this integral ?
1 次查看(过去 30 天)
显示 更早的评论
L0 = 0:20
M0 = -20:20
N0 = 0:20
C = 10 ,
b = 0 ,
V = 2/3 ,
V1 = 2.0001/3 ,
K = 2*pi ,
"EpsilonN : E(n=0) = 1 and for E(n>=1)=2 ",
EpsilonN = [1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2]
%-------------------------------------------------------------------------------------------------------
for L1 = 0:length(L0)
for M1 = 0:length(M0)
for N1 = 0:length(N0)
if (L1/V1)+M1 ~= (N1/V0) %SA1 ~= abs(SA2)
Clmn = (EpsilonN/2*V0*pi)*besselj(M1,K*B)*((sin(((L1/V1)+M1+(N1/V0))*V0*pi)/(L1/V1)+M1+(N1/V0))+...
(sin(((L1/V1)+M1-(N1/V0))*V0*pi)/(L1/V1)+M1-(N1/V0))); %<part 1>
else N1 ~= 0 %n ~= 0 ,
Clmn = (EpsilonN/2)*besselj(M1,K*B) ; %<part 3>
if (L1/V1)+M1 == (N1/V0) %SA1 == SA2 ,
Clmn = EpsilonN*besselj(M1,K*B); %<part 2>
end
end
end
end
end
Clmn
1 个评论
Walter Roberson
2021-3-3
What is the difference between this and https://www.mathworks.com/matlabcentral/answers/725422-find-an-and-bl-unknown-coefficient?s_tid=srchtitle ?
回答(1 个)
Walter Roberson
2021-2-3
format long g
L0 = 0:20;
M0 = -20:20;
N0 = 0:20;
C = 10;
b = 0;
V = 2/3;
V1 = 2.0001/3;
K = 2*pi;
%"EpsilonN : E(n=0) = 1 and for E(n>=1)=2 ",
EpsilonN = 2 * ones(size(N0));
EpsilonN(N0 == 0) = 1;
%-------------------------------------------------------------------------------------------------------
for L1 = 1:length(L0)
L = L0(L1);
for M1 = 1:length(M0)
M = M0(M1);
for N1 = 1:length(N0)
N = N0(N1);
EpsN = EpsilonN(N1);
if abs((L/V1)+M) ~= abs(N/V)
C = (EpsN/2*V*pi)*besselj(M,K*b)*((sin(((L/V1)+M+(N/V))*V*pi)/(L/V1)+M+(N/V))+...
(sin(((L/V1)+M-(N/V))*V*pi)/(L/V1)+M-(N/V))); %<part 1>
elseif N ~= 0
C = (EpsN/2)*besselj(M,K*b) ; %<part 3>
else
C = EpsN*besselj(M,K*b); %<part 2>
end
Clmn(L1, M1, N1) = C;
end
end
end
Clmn
Clmn =
Clmn(:,:,1) =
NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN 1 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324541345203 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324538639666 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324534130437 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324527817516 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324519700904 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.0002193245097806 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324498056251 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.00021932448452892 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324469197819 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324452062476 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324433124846 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324412381271 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324389834943 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324365484953 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324339331957 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324311374587 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324281613548 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324250049389 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000219324216680976 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000219324181508891 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Clmn(:,:,2) =
NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649082691512 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.0004386490772802 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649068260572 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438649055634953 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649039401948 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438649019561557 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.00043864899611285 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648969058151 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648938395602 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648904126131 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648866249274 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648824761776 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648779670612 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648730970667 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648678662871 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648622750013 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648563227445 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648500098885 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648433362474 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648363017747 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Clmn(:,:,3) =
NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN 0 NaN NaN NaN NaN NaN 0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438649082690117 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649077279735 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438649068261502 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649055636348 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438649039402413 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438649019561557 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648996111919 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648969057221 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648938394672 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648904125201 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648866249739 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 -0.000438648824761776 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.000438648779670612 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
3 个评论
Arash Shahpasand
2021-2-4
Dear Walter Roberson,
I hope you are doing well and healthy,
thank you very very very much for the quick response.
when get calculation of summation , i get error "Error using symengine Singularity."
syms l m
AS1 = symsum(symsum(Clmn.*besselj((l/V1)+m,K*C),m,-20,20),l,0,20)
%--------------------------------------------------------------------------------------------------------------
syms l m B
AS2 = symsum(B*symsum(Clmn.*besselj((l/V1)+m,K*C),m,-20,20),l,0,20)
and for this calculation to find An and Bl, need your help, Plz .
"linear algebra equations of the unknown coefficient A(n) and B(l)"
an = (2/V).*EpsilonN.*exp((-1i.*N0.*pi)/(2.*V)).*cos(N0.*gama/V)
AX=B ,
X = [Bl ; An]
A = [A11 A12 ; A21 A22]
B = [R11; R22]
X = linsolve(A,B)
or use "X = inv(A)*B "
Regards
Walter Roberson
2021-2-5
Do not use symsum() for that purpose. Calculate an array of besselj((l/V1)+m,K*C) values that is length(L0) by length(M0), and use .* against Clmn, and sum() that across the second dimension. Construct a vector of B values that is a column length(L0) tall, and .* that by the result of the sum, and sum() the result of the multiplication across the first dimension. The result should be 1 x 1 x length(N0) .
Walter Roberson
2021-2-5
I do not want to get involved in solving infinite numbers of linear equations. Any such a proposal is numeric nonsense, and needs to be approached through theoretical techniques, which might include:
- "renormalization" of infinities -- the sort of mathematics that "proves" that 1-1+1+1+1... infinity "equals" 1/12
- limit processes... which would not give you the individual values anyhow, since there are an infinite number of individual values
- switching from summation to integration (which still will not give you all infinite number of results.)
The determinant of an infinite matrix is going to be either 0 or infinite or exactly 1. In the first two cases, you cannot find the solution for the equations; in the third case, the solution is the matrix times the marginal vector and no simultaneous equation work needs to be done.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Bessel functions 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!发生错误
由于页面发生更改,无法完成操作。请重新加载页面以查看其更新后的状态。
您也可以从以下列表中选择网站:
如何获得最佳网站性能
选择中国网站(中文或英文)以获得最佳网站性能。其他 MathWorks 国家/地区网站并未针对您所在位置的访问进行优化。
美洲
- América Latina (Español)
- Canada (English)
- United States (English)
欧洲
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
亚太
- Australia (English)
- India (English)
- New Zealand (English)
- 中国
- 日本Japanese (日本語)
- 한국Korean (한국어)