查看: 203|回復: 0|關注: 0
打印 上一主題 下一主題

辽宁快乐12怎么玩法: [未答復] 如何對一個含有矩陣的函數進行積分上下限為負無窮到正無窮的積分

[復制鏈接]

新手

6 麥片

財富積分


辽宁快乐12预测一号码推荐号码 www.exedlj.com.cn 050


2

主題

10

帖子

0

最佳答案
跳轉到指定樓層
1#
發表于 2020-2-18 17:13:45 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
我嘗試了三種積分方式,但是都不能得到滿意的答案,請求大佬幫助:handshake
第一種方式:(用for循環)

clc
clear all
M=1024*1;  %橫向分辨率與縱向分辨率
L=1*10^3; %入射場直徑  1*10^4個波長=5.32mm
[x,y]=meshgrid(-L/2:L/(M-1):L/2,L/2:-L/(M-1):-L/2);%空間坐標(-5 : 10/1023 : 5 , 5 : -10/1023 : -5)
r=sqrt(x.^2+y.^2);%為什么./(L/4)
p=0.1;
c=0;
chi=0;

for s=-inf:inf
    Pe=exp(i.*(s.^4+(r./p).*s.^2+(chi.*r./p+c).*s));
end

警告: FOR 循環索引太大。將截斷為 9223372036854775807。
> In Untitled (line 12)


第二種方式:用(int)

clc
clear all
M=1024*1;  %橫向分辨率與縱向分辨率
L=1*10^3; %入射場直徑  1*10^4個波長=5.32mm
[x,y]=meshgrid(-L/2:L/(M-1):L/2,L/2:-L/(M-1):-L/2);%空間坐標(-5 : 10/1023 : 5 , 5 : -10/1023 : -5)
r=sqrt(x.^2+y.^2);%為什么./(L/4)
p=0.1;
c=0;
chi=0;

syms s r
Pe=int(exp(i.*(s.^4+(r./p).*s.^2+(chi.*r./p+c).*s)),s,-inf,inf)


結果:Pe =

int(exp(s^4*1i + r*s^2*10i), s, -Inf, Inf)
這根本沒有運算積分呀,最后想得到double型的結果,這樣的結果也不能用double函數:'(
第三種方法:用(integral)

clc
clear all
M=1024*1;  %橫向分辨率與縱向分辨率
L=1*10^3; %入射場直徑  1*10^4個波長=5.32mm
[x,y]=meshgrid(-L/2:L/(M-1):L/2,L/2:-L/(M-1):-L/2);%空間坐標(-5 : 10/1023 : 5 , 5 : -10/1023 : -5)
r=sqrt(x.^2+y.^2);%為什么./(L/4)
p=0.1;
c=0;
chi=0;

fun = @(s) exp(i.*(s.^4+(r./p).*s.^2+(chi.*r./p+c).*s));
q = integral(fun,0,1)


錯誤使用  .*
矩陣維度必須一致。


出錯 Untitled1>@(s)exp(i.*(s.^4+(r./p).*s.^2+(chi.*r./p+c).*s))


出錯 integralCalc/iterateScalarValued (line 314)
                fx = FUN(t);


出錯 integralCalc/vadapt (line 132)
            [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);


出錯 integralCalc (line 75)
        [q,errbnd] = vadapt(@AtoBInvTransform,interval);


出錯 integral (line 88)
Q = integralCalc(fun,a,b,opstruct);


出錯 Untitled1 (line 13)
q = integral(fun,0,1)

請問怎樣對一個含有矩陣的函數進行積分上下限為負無窮到正無窮的積分?
您需要登錄后才可以回帖 登錄 | 注冊

本版積分規則

關閉

站長推薦上一條 /4 下一條

快速回復 辽宁快乐12预测一号码推荐号码 返回列表