2011-04-09 4 views
0

方程式の記述に問題があります。Matlabの方程式の問題

R = 25、K = 2、R = 50:25:600、DR = 0.5:0.5:4.0

H = Rの*の[1-COS(ASIN((SQRT(2 * R * DR + r * sin(acos(rk)/ r))] - k

しかし、resaultとして、これはh = 1.9118e + 001 + 1.7545e + 002iとなる。

私はMatlabから始めます。私はあなたが書いたものから、何を得るのおかげ

+0

結果は何をすべきですか? – BlackBear

+0

私は計算機を使用したとき、結果として私は4,56を得た – user699994

答えて

2

は実際にあなたが互いことにより、2つの行ベクトルを乗算しようとしているので、正しい

??? Error using ==> mtimes 
Inner matrix dimensions must agree. 

です。実際に使用したコードを表示してください。

とにかく、それが何とか処理されているとすれば、それは1よりもはるかに大きいasinに何かを与えているかのように見えます。それはあなたに複雑な結果を与えるでしょう。あなたが渡しているものが、おそらくR^2またはDR^2などで分けられることを意味するasinにありますか?ちょっと後で、引数がacosの類似の問題があります。

は、私はまた、あなたの*^/事業者のいくつかは、実際に要素ごとのもの.*.^./べきであると思われます。

+0

まあ、多分私の質問は少し混乱しています。その結果、私はテーブル(またはマトリックス)を取得する必要があるので、最初の式では、R = 50、DR = 0.5、r = 25、k = 2を使用してhを取得する必要があります。第2の式において、R = 600、DR = 4.0、r = 25、k = 2を使用する最後の式iについて、R = 75、DR = 1.0、r = 25、k = 2 ...を使用する。 – user699994

+0

'acos(r-k)'と同じこと、たぶんacos((r-k)/ r) ' –

0

EDIT:コメントが正しく反映されるように、わずかに調整された式(ブラケッティング)。

あなたはテーブルがほしいと言いますと、それはRDRテーブルであると思います(異なる長さのベクトルが必要なので)。これを行うには、Rを列ベクトル(R'以下)として使用し、*.*ではない)を掛ける必要があります。 Rが、ones(size(R))(またはrepmatを使用)を乗算して正しい形にDRすることができない場合。 DRを要素で四角にするには、DR.^2が必要です。 acosの括弧が間違っているようですが、acosを受け取る前に必ずrで割ります。 asinにはrなどの除算が必要です(sqrtになっているため、r^2ではありません)。最後に、rによる最後の除算は、直前に同じレベルでrを掛けているため、書かれているように冗長です。とにかく、私がしなければ、次の

h= r*(1-cos(asin((sqrt(2*R'*DR+ones(size(R))'*DR.^2)/r)+sin(acos((r-k)/r)))))-k 

は私が R DRによってテーブルを取得します。小さい R,DRの結果は実数です。より高い R,DRは、最初の asinの引数が> 1であるために複雑です。必要に応じて、テーブルの最初のエントリは4.56です。

+0

と言っていただき、ありがとうございました。私はちょうどr = 25、k = 2、R = 50、DR = 0.5を取って、この式に含まれていましたh = r * [1-cos(asin((sqrt(2 * R * DR + DR。 r)+(r * sin(acos(rk)/ r))/ r)))] - kとなり、結果として4.56となった。 – user699994

0

あなたが言ったようにやろうとしている場合:

ので最初の式Iにおいて、R = 25、DR = 0.5、R = 50を使用し、K = 2、iは時間を取得する必要があります。第2の式では、R = 75、 DR = 1.0、r = 25、k = 2 ...最後に 式iが使用された。 R = 600、DR = 4.0、r = 25、k = 2。

DR及びRは、Rは、25刻みで50と600との間になった場合、DRは0.5〜0.5単位で 12.5 0.5から行く、または必要がありますので...同じ長さである必要が

0.1522刻みで4.0 ...

あなたがそれを把握たら、(例えば*または^)全ての行列乗算演算の前にピリオドを追加してください

関連する問題