関数ハンドルのハンドルの渡し方について

Question

kazuma kaneda 2021-11-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1590714-

回答： Atsushi Ueno 2021-11-22

采纳的回答： Atsushi Ueno

下のように、ωについての積分をするコードを書きました。

・・・・１

と

は以下の通りです。

・・・・２

以下がコードです。

A=1;m=1;s=1;wA=1;const=10^-10; %定数を決定

%２の計算　まず、eで積分

J=ffun(@(e)A.*exp(-(e-m).^2./(2*pi.*s.^2))./(w-e),-inf,w-const,w+const,inf)-@(w)1i*pi*A.*exp(-(w-m).^2./(2*pi.*s.^2));

%１の計算　次に、wで積分

K=@(t)gfun(@(w)exp(-1i*w.*t)./(w-wA-J),-inf,inf);

function f=ffun(fun,a,b,c,d)

f=integral(fun,a,b)+integral(fun,c,d);

end

function g=gfun(fun,a,b)

g=integral(fun,a,b);

end

こなままだと、まずωが変数として認識されません。また、ハンドルの渡し方についても違うと思うのですがわかる方がいましたらご教授お願いします。

2 个评论
显示无隐藏无

Atsushi Ueno 2021-11-21

2行目の”@(w)”を先頭に移動すれば、ωに関するエラー無く動く事は間違いないです。

ただ、続くKの演算におけるgfunの中身(積分する式)が「ωを含む数式-(ωを引数とする関数ハンドル)」と、数式の内部で関数ハンドルを演算しており、この形が成立するのかどうかが疑問です

Atsushi Ueno 2021-11-22

無名関数 - MATLAB & Simulink - MathWorks 日本：複数の無名関数

すなわち「無名関数の入れ子」の考え方が必要ですね

请先登录，再进行评论。

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

Answer 1

Atsushi Ueno 2021-11-22

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1590714-#answer_837474

ηによる積分範囲を分割して下記のようにした理由が良く分かりませんが、とにかくそのように致します。

JとKを離せない問題があって長く見難いですが、これで引数t,ω,ηが入れ子になった無名関数として成立します。

A=1;m=1;s=1;wA=1;const=10^-10; %定数を決定
K=@(t)gfun(@(w)exp(-1i*w.*t)./(w-wA-ffun(@(e)A.*exp(-(e-m).^2./(2*pi.*s.^2))./(w-e),-inf,w-const,w+const,inf)-1i*pi*A.*exp(-(w-m).^2./(2*pi.*s.^2))),-inf,inf)
K = function_handle with value:
    @(t)gfun(@(w)exp(-1i*w.*t)./(w-wA-ffun(@(e)A.*exp(-(e-m).^2./(2*pi.*s.^2))./(w-e),-inf,w-const,w+const,inf)-1i*pi*A.*exp(-(w-m).^2./(2*pi.*s.^2))),-inf,inf)
function f=ffun(fun,a,b,c,d)
f=integral(fun,a,b)+integral(fun,c,d);
end
function g=gfun(fun,a,b)
g=integral(fun,a,b);
end