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

編集履歴

質問編集履歴

コードの編集

2020/02/14 13:54

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -74,7 +74,7 @@
     X,Y=np.meshgrid(np.linspace(-5,5,100),np.linspace(-5,5,100))
-    t,_=np.array([X,Y])
+    t=np.array([X,Y])
     Z,_=func(t)
@@ -182,7 +182,25 @@
 ```error
+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")
-TypeError: Input z must be at least a 2x2 array.
+ValueError: matrix must be 2-dimensional
 ```

文法の修正をしました。

2020/02/14 13:54

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -10,13 +10,21 @@
 ```python
+import numpy as np
+import matplotlib.pyplot as plt
-　rosenbrockfunction = lambda x,y: (1-x)**2+100*(y-x**2)**2
+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))
+X,Y=np.meshgrid(np.linspace(-5,5,100),np.linspace(-5,5,100))
-　Z=rosenbrockfunction(X,Y)
+Z=rosenbrockfunction(X,Y)
-　cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
+cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
+plt.show()
 ```

タイトルの変更

2020/02/14 08:20

投稿

alpha_high

スコア6

test CHANGED Viewed

	@@ -1 +1 @@
1	- パイソンにて、等高線のプロットがしたいです。~~１００x１００行列~~の~~X, １００x１００行列のYから、１００x１００のzの式を出したいの~~ですが、配列の要素数が合いません。
1	+ パイソンにて、等高線のプロットがしたいです。私のプログラムでは、配列の要素数が合わないというエラーが出ています。

test CHANGED Viewed

File without changes

タイトルの編集をしました。

2020/02/14 04:20

投稿

alpha_high

スコア6

test CHANGED Viewed

	@@ -1 +1 @@
1	- パイソンにて、等高線のプロットがしたいです。１００x１００行列のX, １００x１００行列のYから、１００x１００のzの式を出したいのですが、配列の要素すうが合いません。
1	+ パイソンにて、等高線のプロットがしたいです。１００x１００行列のX, １００x１００行列のYから、１００x１００のzの式を出したいのですが、配列の要素数が合いません。

test CHANGED Viewed

File without changes

隙間を消しました。

2020/02/14 01:19

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -34,9 +34,7 @@
-```uncon
+```python
 import math
@@ -170,16 +168,6 @@
-```
-```
 このまま実行をすると、

試したことを足しました。

2020/02/14 01:06

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -220,4 +220,14 @@
+また、Zを１００x１００行列で静的行列にしてから計算しても、出力されませんでした。
+```
+Z=[[0 for I in range(100)] for j in range(100)]
+```
 私のコードの間違いを教えていただけないでしょうか。

コードを動くように修正しました。申し訳ありません。

2020/02/14 01:06

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -38,11 +38,33 @@
+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))
@@ -114,6 +136,8 @@
 nTrails=10
+testFuncs = [rosenbrock]
 for _, func in enumerate(testFuncs):

追加で編集しました。

2020/02/14 00:33

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -35,6 +35,32 @@
 ```uncon
+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()
@@ -116,24 +142,6 @@
         sys.exit(1)
-        l=xopt[0,0]
-        m=xopt[0,1]
-        xopt=[l,m]
-        _, gopt =  rosenbrock(np.array(xopt))
-        if (not (np.linalg.norm(gopt, ord = np.inf) <= 1e-6)):
-            continue
-        print("---------------result------------------")
-        print('x_star',xopt)
-        print('fopt ', fopt, FUNCCALLS)
 ```
@@ -142,27 +150,7 @@
-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()

エラーの追加をしました。

2020/02/14 00:30

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -154,40 +154,48 @@
     Z=np.matrix(Z)
-    print("length",X.shape[1],Y.shape[1],Z.shape[0],Y.shape[1])
+    #print("length",X.shape[1],Y.shape[1],Z.shape[0],Y.shape[1])
-    print("X",X , "Y", Y, "Z", Z)
+    #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()
+```
+このまま実行をすると、
+```error
+TypeError: Input z must be at least a 2x2 array.
+```
+とエラーが出てしまいます。
+XとYのdimension数を合わせなければいけないと思い、
+```python
+Z.reshape(X.shape[0],X.shape[1])
+```
+のように、配列数を合わせようとしてみたのですが、
+```error
     Z.reshape(X.shape[0],Y.shape[1])
-    cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
-    plt.show()
-```
-ここで、
-```python
-Z.reshape(X.shape[0],X.shape[1])
-```
-のように、配列数を合わせようとしてみたのですが、
-```error
-    Z.reshape(X.shape[0],Y.shape[1])
 ValueError: cannot reshape array of size 100 into shape (100,100)
 ```

コードを追加しました。

2020/02/13 23:25

投稿

alpha_high

スコア6

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -34,71 +34,167 @@
+```uncon
+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
+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)
+        l=xopt[0,0]
+        m=xopt[0,1]
+        xopt=[l,m]
+        _, gopt =  rosenbrock(np.array(xopt))
+        if (not (np.linalg.norm(gopt, ord = np.inf) <= 1e-6)):
+            continue
+        print("---------------result------------------")
+        print('x_star',xopt)
+        print('fopt ', fopt, FUNCCALLS)
+```
+```
+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()
+```
+ここで、
 ```python
-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
-funcs = [rosenbrock]
-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]X.shape[1])
-    cs=plt.contour(X,Y,Z,np.logspace(0,3.5,7,base=10))
-    plt.show()
-plot(func)
-```
-ここで、
-```python
 Z.reshape(X.shape[0],X.shape[1])
 ```
 のように、配列数を合わせようとしてみたのですが、
+```error
+    Z.reshape(X.shape[0],Y.shape[1])
+ValueError: cannot reshape array of size 100 into shape (100,100)
+```
-１００のarray を１００x１００に出来ませんと出てきてしまい、上下のプログラムで、関数を分けて居るかどうかだけなのに、何故答えが出せないのか、と悩んでいます。
+上下のプログラムで、関数を分けて居るかどうかだけなのに、何故答えが出せないのか、と悩んでいます。
 リストや行列への変換も、考えられるペアで変えてみましたが、出来ませんでした。