1

このトピックで説明されているのと同じ種類の問題があります: Using fzero: タイプ 'sym' の入力引数に対して未定義の関数またはメソッド 'isfinite'

彼らの答えは本当に助けになりましたが、私はまだ立ち往生しています。

w の関数の零点も見つけなければなりません。この関数はいくつかのステップで定義されます。

したがって、唯一の未知数は w であり、次のような他のオブジェクトを定義しました。

lambda= @(w) ((16*rho(i)*A(i)*w^2*Lprime(i)^2)/(E(j)*I(i)))^0.25;

beta=@(w) lambda*b(i)^0.5;

gamma=@(w) lambda*Lprime(i)^0.5;

次に、4*4 行列 M2 を定義します。

M2=@(w) [besselj(4,beta) bessely(4,beta) besseli(4,beta) besselk(4,beta);
               besselj(3,beta) bessely(3,beta) besseli(3,beta) -besselk(3,beta);
               besselj(2,gamma) bessely(2,gamma) besseli(2,gamma) besselk(2,gamma);
               besselj(4,gamma) bessely(4,gamma) besseli(4,gamma) besselk(4,gamma)];

この場合、解く方程式は det(M2)=0 です。しかし、w=0 は解の 1 つです。最初のゼロ以外の解が必要なので、次のように書きました。

delta = @(w) det(M2);

S(i,j)=fzero(delta,500);

次に、プログラムを実行すると、Matlab は次のように言います。

??? Error using ==> fzero at 235
FZERO cannot continue because user supplied function_handle ==> @(w)det(M2)
failed with the error below.

Undefined function or method 'det' for input arguments of type 'function_handle'.

Error in ==> frequencies at 57
    S(i,j)=fzero(delta,500);

subs と eval メソッドも試しましたが、どちらも機能しません。これらの場合のエラーメッセージは次のとおりです。

??? Undefined function or method 'isfinite' for input arguments of type 'sym'.

Error in ==> fzero at 323
    elseif ~isfinite(fx) || ~isreal(fx)

Error in ==> frequencies at 58
    S(i,j)=fzero(@(w) subs(delta,'w',w),500);

これは、私が推測するedioと同じエラーです。と：

??? Error using ==> fzero at 307
FZERO cannot continue because user supplied function_handle ==> @(w)eval(delta)
failed with the error below.

Undefined function or method 'eval' for input arguments of type 'function_handle'.

Error in ==> frequencies at 59
    S(i,j)=fzero(@(w)eval(delta),500);

助けてください。

4

2 に答える 2

3

あなたの問題は、匿名関数を他の匿名関数内に配置するときに、匿名関数を評価していないことのようです。たとえば、次のように関数を定義しますlambda。

lambda = @(w) ((16*rho(i)*A(i)*w^2*Lprime(i)^2)/(E(j)*I(i)))^0.25;

しかし、 で使用する場合は、次のようbetaに の入力値を使用して評価する必要があります。w

beta = @(w) lambda(w)*b(i)^0.5;
                %# ^--------------Pass w to lambda to evaluate the function

そのため、他の匿名関数は次のように定義する必要があると思います。

gamma = @(w) lambda(w)*Lprime(i)^0.5;

M2 = @(w) [besselj(4,beta(w)) bessely(4,beta(w)) besseli(4,beta(w)) ...
           besselk(4,beta(w)); ...
           besselj(3,beta(w)) bessely(3,beta(w)) besseli(3,beta(w)) ...
           -besselk(3,beta(w)); ...
           besselj(2,gamma(w)) bessely(2,gamma(w)) besseli(2,gamma(w)) ...
           besselk(2,gamma(w)); ...
           besselj(4,gamma(w)) bessely(4,gamma(w)) besseli(4,gamma(w)) ...
           besselk(4,gamma(w))];

delta = @(w) det(M2(w));


効率についてのメモ...

私がここで気づいているGLARING効率の問題があります。他のタイプの関数 (プライマリ関数、ネストされた関数、またはサブ関数) の代わりに無名関数を使用すると、同じ入力で同じ関数を複数回評価することになります。

たとえば、M2マトリックスを作成するために評価するたびに、同じ入力でbetaとgamma8 回の両方を評価することになります! M2関数に配置し、入力として渡すことwと、2 つの関数ハンドルbetaとを渡すことで改善できることに注意してくださいgamma。

function newMatrix = M2(w,betaFcn,gammaFcn)

  bw = betaFcn(w);   %# Evaluate the beta function once
  gw = gammaFcn(w);  %# Evaluate the gamma function once
  newMatrix = [besselj(4,bw) bessely(4,bw) besseli(4,bw) besselk(4,bw); ...
               besselj(3,bw) bessely(3,bw) besseli(3,bw) -besselk(3,bw); ...
               besselj(2,gw) bessely(2,gw) besseli(2,gw) besselk(2,gw); ...
               besselj(4,gw) bessely(4,gw) besseli(4,gw) besselk(4,gw)];

end

新しいdelta関数は次のようになります。

delta = @(w) det(M2(w,beta,gamma));
于 2011-05-26T18:34:16.543 に答える
0

こんにちは、ご協力ありがとうございます。それは機能しますが、明らかに最後の行を変更する必要があります（それを理解するのにまだ10分かかりました）：

lambda= @(w) ((16*rho(i)*A(i)*w^2*Lprime(i)^2)/(E(j)*I(i)))^0.25;
beta=@(w) lambda(w)*b(i)^0.5;
gamma=@(w) lambda(w)*Lprime(i)^0.5;


M2=@(w) [besselj(4,beta(w)) bessely(4,beta(w)) besseli(4,beta(w)) besselk(4,beta(w));
    besselj(3,beta(w)) bessely(3,beta(w)) besseli(3,beta(w)) -besselk(3,beta(w));
    besselj(2,gamma(w)) bessely(2,gamma(w)) besseli(2,gamma(w)) besselk(2,gamma(w));
    besselj(4,gamma(w)) bessely(4,gamma(w)) besseli(4,gamma(w)) besselk(4,gamma(w))]; 

delta = @(w) det(M2(w));
    S(i,j)=fzero(@(w) delta(w),500);

そして今では以前よりも本当に速くなりました。関数solveが解決を処理できる別のケースでは、ループごとに10秒ほどかかりましたが、今では0.06秒のようになっています。

他の解決策を試して、改善点を確認します。

大いに感謝する。

于 2011-05-27T04:11:44.407 に答える