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

4 次查看（过去 30 天）
kazuma kaneda
kazuma kaneda 2021-11-19
回答： Atsushi Ueno 2021-11-22
下のように、ωについての積分をするコードを書きました。
・・・・１
とは以下の通りです。
・・・・２
以下がコードです。
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
Atsushi Ueno 2021-11-21
2行目の”@(w)”を先頭に移動すれば、ωに関するエラー無く動く事は間違いないです。
ただ、続くKの演算におけるgfunの中身(積分する式)が「ωを含む数式-(ωを引数とする関数ハンドル)」と、数式の内部で関数ハンドルを演算しており、この形が成立するのかどうかが疑問です
Atsushi Ueno
Atsushi Ueno 2021-11-22
無名関数 - MATLAB & Simulink - MathWorks 日本：複数の無名関数
すなわち「無名関数の入れ子」の考え方が必要ですね

请先登录，再进行评论。

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

采纳的回答

Atsushi Ueno
Atsushi Ueno 2021-11-22
Ran in:
ηによる積分範囲を分割して下記のようにした理由が良く分かりませんが、とにかくそのように致します。
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

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 関数の基礎 的更多信息

标签

产品


版本

R2021b

Community Treasure Hunt

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

Start Hunting!