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

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

ただいまの
回答率

89.63%

パイソンにて、等高線のプロットがしたいです 。私のプログラムでは、配列の要素数が合わないというエラーが出ています。

受付中

回答 2

投稿 編集

  • 評価
  • クリップ 1
  • VIEW 157

alpha_high

score 1

パイソンを使って等高線のグラフを書いております。

ですが、配列の数が合わないというエラーでうまく書くことができません。

下記のように、プロット関数に直接、関数を書いて出力すると、等高線がかけるのですが、

import numpy as np
import matplotlib.pyplot as plt

rosenbrockfunction = lambda x,y: (1-x)**2+100*(y-x**2)**2
X,Y=np.meshgrid(np.linspace(-5,5,100),np.linspace(-5,5,100))
Z=rosenbrockfunction(X,Y)
cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
plt.show()

ですが、
下記のように、関数を分けると、
_zは二次元行列です
_配列数が合いません
とエラーが出ます。

import math
import utils
import utils2
import numpy as np
import sys
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import sys
import csv
import os


def cont(func):

    X,Y=np.meshgrid(np.linspace(-5,5,100),np.linspace(-5,5,100))
    t=np.array([X,Y])
    Z,_=func(t)
    Z=np.matrix(Z)
    #print("length",X.shape[1],Y.shape[1],Z.shape[0],Y.shape[1])
    #print("X",X , "Y", Y, "Z", Z)
    #Z.reshape(X.shape[0],Y.shape[1])
    cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
    plt.show()


def rosenbrock(x):
    global FUNCCALLS
    FUNCCALLS += 1


    funcValue = (100*((x[1]-x[0]**2)**2) + (1-x[0])**2)
    gradient = np.matrix([(-400*x[0]*(x[1]-x[0]**2)-2*(1-x[0])), (200*(x[1]-x[0]*x[0]))]).T

    return funcValue, gradient


def uncon(func, x0, epsilon_g, options=None):


    if options is None:

        n=2
        xopt, fopt, outputs = optimize(func, x0, epsilon_g,n)

    return xopt, fopt, outputs



nTrails=10
testFuncs = [rosenbrock]

for _, func in enumerate(testFuncs):

    np.random.seed(0)

    trialScore = 0
    for _ in range(nTrails):
        FUNCCALLS = 0

        x0 = np.random.normal(size=2)*3

        #xopt, fopt, output = uncon(func, x0, 1e-6)
        cont(func)
        sys.exit(1)

このまま実行をすると、

Traceback (most recent call last):
  File "stud2.py", line 65, in <module>
    cont(func)
  File "stud2.py", line 18, in cont
    Z,_=func(t)
  File "stud2.py", line 34, in rosenbrock
    gradient = np.matrix([(-400*x[0]*(x[1]-x[0]**2)-2*(1-x[0])), (200*(x[1]-x[0]*x[0]))]).T
  File "/Users/mizuhotakayama/.pyenv/versions/3.8.1/lib/python3.8/site-packages/numpy/matrixlib/defmatrix.py", line 151, in __new__
    raise ValueError("matrix must be 2-dimensional")
ValueError: matrix must be 2-dimensional


とエラーが出てしまいます。
XとYのdimension数を合わせなければいけないと思い、

Z.reshape(X.shape[0],X.shape[1])


のように、配列数を合わせようとしてみたのですが、

    Z.reshape(X.shape[0],Y.shape[1])
ValueError: cannot reshape array of size 100 into shape (100,100)

上下のプログラムで、関数を分けて居るかどうかだけなのに、何故答えが出せないのか、と悩んでいます。
リストや行列への変換も、考えられるペアで変えてみましたが、出来ませんでした。

また、Zを100x100行列で静的行列にしてから計算しても、出力されませんでした。

Z=[[0 for I in range(100)] for j in range(100)]

私のコードの間違いを教えていただけないでしょうか。

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

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • can110

    2020/02/14 08:07

    2つ目のコードはInvalid Syntaxなどで実行できないため修正ください。
    またエラーは全文(Traceback)を記載ください。

    キャンセル

  • alpha_high

    2020/02/14 08:26

    エラーの件、本当に申し訳ありません。少し長くなってしまったのですが、私が動かしている状態で修正させて頂きました。また、エラー 文も追加させて頂きました。見ていただけると本当にありがたいです。

    キャンセル

  • hayataka2049

    2020/02/14 14:27

    エラー文、何十行かだらだらっと出てきませんか? 
    Exception in ***
    Traceback (most recent call last):
    から始まるやつです。あったらそれを含めたエラーに変えておいてください。

    キャンセル

  • alpha_high

    2020/02/14 22:55

    エラー文の載せ方を教えて頂いて、ありがとうございます。変更させて頂きました。

    キャンセル

回答 2

0

Python のエラー (Traceback) ってこういうやつです。(Jupyter だともっと長くなるけど。)

Traceback (most recent call last):
  File "contour.py", line 63, in <module>
    cont(func)
  File "contour.py", line 22, in cont
    cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
  File "/usr/local/lib/python3.7/site-packages/matplotlib/pyplot.py", line 2514, in contour
    **kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/__init__.py", line 1599, in inner
    return func(ax, *map(sanitize_sequence, args), **kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/axes/_axes.py", line 6422, in contour
    contours = mcontour.QuadContourSet(self, *args, **kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/contour.py", line 855, in __init__
    kwargs = self._process_args(*args, **kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/contour.py", line 1456, in _process_args
    x, y, z = self._contour_args(args, kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/contour.py", line 1514, in _contour_args
    x, y, z = self._check_xyz(args[:3], kwargs)
  File "/usr/local/lib/python3.7/site-packages/matplotlib/contour.py", line 1551, in _check_xyz
    raise TypeError("Input z must be at least a 2x2 array.")
TypeError: Input z must be at least a 2x2 array.

これで何がわかるかというと、まず 63 行目で cont 関数を呼び出し、

  File "contour.py", line 63, in <module>
    cont(func)

次に 22 行目で plt.contour を呼び出し、

  File "contour.py", line 22, in cont
    cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))

最終的に matplotlib の中でエラーになってます。

  File "/usr/local/lib/python3.7/site-packages/matplotlib/contour.py", line 1551, in _check_xyz
    raise TypeError("Input z must be at least a 2x2 array.")
TypeError: Input z must be at least a 2x2 array.

で、z 引数がおかしいということなので、

    Z=np.matrix(Z)
    print(Z.shape)

としたところ、次のように表示されました。

(1, 100)

これは 1x100 行列なので、「2x2 以上」ではありませんね。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/02/14 22:59

    回答本当にありがとうございます。エラーの、2x2以上という意味がわかりません。
    X、Yが100x100行列に対し、何故Zは1x100行列なのでしょうか。また、どのように出力すれば、Zも100x100となるのでしょうか。教えていただけると本当に嬉しいです。

    キャンセル

0

t,_=np.array([X,Y])でアンパックしたあとのtを渡しているから。

また、勾配の計算もうまくいかないようですが、そのコードでは使っていないようなのでさしあたってはコメントアウトしておくと良いかと思います。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/02/14 22:58

    t=np.array([X,Y])
    とする、という意味で合っていますでしょうか。
    私のコードでは、同じエラーが出ています。ZをX,Yと同じ形にするには、どうすれば良いのでしょうか。
    XとYが100x100に対して、Zが1x100しか値がない理由がわかりません。どの部分で、うまく計算ができていないか、教えていただけないでしょうか。

    キャンセル

  • 2020/02/14 23:17

    その変更を行った場合エラーが発生するのはgradientの行です。これは今回の実行には関係ないのでコメントアウトし、返り値のフォーマットを変えて呼び出し元でも対応させるか、Noneなどを返却すればよいでしょう。

    キャンセル

  • 2020/02/15 02:49

    回答本当にありがとうございます。質問に答えて頂き、本当にありがとうございます。rosenbrockの中身の関数は、他の計算の関係上消すことができないのですが、その場合だとグラフをfuncを使って表示させる事は難しいでしょうか。。個人的には、Z,_=func(t) としていれば、gradient の値は受け取らないものだと思っていたのですが、このアンダーバーがZの配列に影響しているのでしょうか。

    キャンセル

  • 2020/02/15 03:06

    重ねてすみません。rosenbrock 関数の中身のgradient にてnoneを返してみたのですが、同じエラーが発生しました。ですが、コメントアウトをしたところ、グラフが表示できました。引数が2つあると、配列の形が変わってしまうのですね...どうにかして、rosenbrockの関数を中身をそのままで、グラフ表示というのは厳しいでしょうか。。同時に他のグラフも表示させたいと思っていたので、できればfunc関数でグラフを表示させたいです。何度も本当に申し訳ありませんが、お知恵をお貸し頂けると嬉しいです。

    キャンセル

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

  • ただいまの回答率 89.63%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる
  • トップ
  • Pythonに関する質問
  • パイソンにて、等高線のプロットがしたいです 。私のプログラムでは、配列の要素数が合わないというエラーが出ています。