質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.50%
Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Q&A

解決済

2回答

35745閲覧

【python】numpy ブロードキャストエラー

mtnow

総合スコア13

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

0グッド

0クリップ

投稿2019/09/18 07:44

編集2019/09/18 12:52

numpyブロードキャストエラーの原因について教えてください.
下記のようなエラーです.
ValueError: operands could not be broadcast together with shapes (21,) (51,)

python

1import numpy as np 2import math 3 4Pn = 2200 5z_depth = 5 6π = math.pi 7 8########ここでエラーがでていると思います############ 9################################################ 10y_c = np.arange(-10, 11) / 10 #初項-10,終点11の等差数列 11c=0.5 12y = y_c*c 13################################################ 14 15print('平面応力状態(半無限平板)') 16 17z = np.arange(0,(z_depth+0.1)*10)/10 #z決定 18 19 20σy_heimen = (-2*Pn*(y**2)*z)/(π*(((y**2)+(z**2)))**2) 21σz_heimen = (-2*Pn*(z**3))/(π*(((y**2)+(z**2)))**2) 22τyz_heimen = (-2*Pn*(z**2)*y)/(π*(((y**2)+(z**2)))**2) 23for σy_heimen1,σz_heimen1,τyz_heimen1,z1 in zip(σy_heimen,σz_heimen,τyz_heimen,z): 24 print(f'σy={σy_heimen1:.6f},σz={σz_heimen1:.6f},τyz={τyz_heimen1:.6f},z={z1:.6f}') 25 26 27input() 28

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

tiitoi

2019/09/18 07:58

質問のコードを実行してもエラーが再現しません。 Jupyter Notebook などをお使いであれば、再起動して再度実行していみてください。
guest

回答2

0

自己解決

エラーの出た原因はわかりました.
yとzの数値の個数が違ったため,エラーが出たのだと思います.
改善策は改めて考え,質問させていただくかもしれません.
回答していただいた方々,ありがとうございました.

投稿2019/09/18 13:58

mtnow

総合スコア13

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

jeanbiego

2019/09/19 00:21

y_c = np.arange(-10, 10.4, 0.4) / 10 #初項-10,終点11の等差数列 とすることで、yのsizeをzに合わせてみると動くことは動きました。(求めている結果にふさわしいかというのは別として)
guest

0

numpyのimportが誤字?していたのと、数式用の各定数が定義されていなかったので下記のように修正したところ普通に計算されました。

python3

1import numpy as np 2import math 3 4Pn = 2200 5z_depth = 5 6π = math.pi 7y=1 8 9z = np.arange(0,(z_depth+0.1)*10)/10 #z決定 10 11σy_heimen = (-2*Pn*(z**3))/(π*(((y**2)+(z**2)))**2) 12 13for σy_heimen1,z1 in zip(σy_heimen,z): 14 print(f'σy={σy_heimen1:.6f},z={z1:.6f}')

ちなみに、Broadcastとは、サイズ違いのアレイ同士を計算する際に、自動的に補正してくれる機能のようです。
その際、各次元のサイズが同じか、あるいは片方が1でないといけないようです。
各定数を、別の所で変に代入したままだったのではないでしょうか…

参考URL
Array Broadcasting in Numpy

投稿2019/09/18 08:18

jeanbiego

総合スコア3966

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

mtnow

2019/09/18 12:48

エラーが起こる箇所をみつけたので,もしよければ見ていただけないでしょうか. 質問修正します.
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.50%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問