前提・実現したいこと
pythonに関する質問です。
楕円方程式を作成しました。
そのxとyの値をcsvファイルに出力したいのですがどのようにすればよいかわかりません。
発生している問題・エラーメッセージ
x,yの値が表示されない
該当のソースコード
コード import numpy as np import csv from matplotlib import pyplot as plt #楕円方程式 def ellipse(x,a,b):#円 範囲を表す y=b*np.sqrt(1-(np.power(x,2)/np.power(a,2)))#楕円方程式をy=にした形 sqrtは平方根 return y a=1.0 b=0.5 t=np.linspace(-1*np.pi,np.pi,100)#-π~πまでの範囲を100等分した等差数列 phase=np.pi/2#位相 x=a*np.sin(t-phase)#1周期分の正弦波 g=np.gradient(x)#Xの勾配計算 符号の変化 y=[]#全象限のy値計算 for i in range(len(g)): y=np.append(y,ellipse(x[i],a,b))#appendで末尾に要素を追加している if i==0: pass elif np.sign(g[i-1])==np.sign(g[i]): pass else: y=-1*y #もしi=0ならなにもしない、それ以外はg[i-1]とg[i]が同じならなにもしない #それ以外はy=-1*yを実行 f=open('homework1.csv','w',newline='') writer=csv.writer(f) writer.writerow(['x','y']) #グラフ描画 #フォントサイズ、種類 plt.rcParams['font.size']=14 plt.rcParams['font.family']='Times New Roman' #目盛をうちがわにする plt.rcParams['xtick.direction']='in' plt.rcParams['ytick.direction']='in' #目盛線をつける fig=plt.figure() ax1=fig.add_subplot(111) ax1.yaxis.set_ticks_position('both') ax1.xaxis.set_ticks_position('both') #軸のラベルを設定する ax1.set_xlabel('x') ax1.set_ylabel('y') #データの範囲と刻み目盛を明示する ax1.set_xticks(np.arange(-2,2,0.5)) ax1.set_yticks(np.arange(-2,2,0.5)) ax1.set_xlim(-1,1) ax1.set_ylim(-1,1) #データプロットの準備とともに、ラベルと線の太さ、凡例の設置を行う ax1.plot(x,y,label='ellipse',lw=1) ax1.set_aspect('equal') fig.tight_layout() #グラフを表示する plt.show() plt.close()
補足情報(FW/ツールのバージョンなど)
Spyder(Python3.7)
回答1件
あなたの回答
tips
プレビュー