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

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

新規登録して質問してみよう
ただいま回答率
85.54%
CSV

CSV(Comma-Separated Values)はコンマで区切られた明白なテキスト値のリストです。もしくは、そのフォーマットでひとつ以上のリストを含むファイルを指します。

JupyterLab

JupyterLabは、Jupyter notebookの後継の対話型開発環境(IDE)です。データの可視化がインタラクティブで、プラグイン作成により新しいコンポーネントの追加および既存のコンポーネントも統合可能。サーバに閉じているため、データ分析に向いています。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

2回答

565閲覧

pythonでif文の条件分岐後の結果を足し合わせたいがうまくいかない

NMKN

総合スコア8

CSV

CSV(Comma-Separated Values)はコンマで区切られた明白なテキスト値のリストです。もしくは、そのフォーマットでひとつ以上のリストを含むファイルを指します。

JupyterLab

JupyterLabは、Jupyter notebookの後継の対話型開発環境(IDE)です。データの可視化がインタラクティブで、プラグイン作成により新しいコンポーネントの追加および既存のコンポーネントも統合可能。サーバに閉じているため、データ分析に向いています。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2023/01/07 21:39

前提

pythonでフーリエベッセル関数を用いた級数展開法を行なっています。今基底関数系のデータを決め、それぞれの基底関数系に対応する係数を求めることはできたのですが、この求めた係数を対応する基底関数系に掛けて値を足し合わせるというところで苦労しています。最終的には設定したデータと合うように像を再構成できれば良いので等高線を描くためにcontourを利用しています。

実現したいこと

if文で場合分けしている関数の値を全て足し合わせたい

発生している問題・エラーメッセージ

明らかに異なる結果が出ているのでif文でミスがあるように思われる。27~34行を変える必要があると思われる。

該当のソースコード

python

1import numpy as np 2import matplotlib.pyplot as plt 3from scipy.special import jv 4import pandas as pd 5import itertools 6import csv 7 8csv_file = open('coeffication.csv') #csvファイルを開く 9csv_reader = csv.reader(csv_file) #開いたcsvファイルからreaderオブジェクトを生成 10clf = [] #抽出するデータを格納する空のリストを作る 11 12for row in csv_reader: #readerオブジェクトをループしてデータ抽出 13 if csv_reader.line_num == 1: #ヘッダー行はスキップする 14 continue #Trueになる1行目はなにもしない 15 clf.append(row[0]) #row行目の係数データを格納 16 17csv_file.close() # csvファイルを閉じる。 18 19clf_i = [float(s) for s in clf] #配列は文字列型なのでint型に変換 20 21fig = plt.figure() 22ax1 = fig.add_subplot(111) 23x = np.linspace(0, 4*2^(1+0), 40) #等間隔で値を0~4まで20個作成 lの値はここで変える 24y = np.linspace(0, 2*np.pi, 40) #等間隔で値を0~360まで20個作成 25x, y = np.meshgrid(x, y) #x,yからメッシュグリッドを作成 26 27def fz(x,y): 28 for e in range(41): 29 if e % 2 ==0: 30 a = clf_i[e]*jv(e,x)*np.cos(e*y) 31 else: 32 b = clf_i[e]*jv(e,x)*np.sin(e*y) 33 return a+b #足し合わせる 34z = fz(x,y) 35 36x1 = x*np.cos(y) 37y1= x*np.sin(y) 38 39ax = plt.contour(x1,y1,z,colors='black') #等高線表示 40ax.clabel(fmt='%1.1f', fontsize=16) #等高線の値を表示 41ax = plt.contourf(x1,y1,z,cmap='rainbow') #等高線レベルに応じて色を塗る 42ax = plt.colorbar(label="contour level") #カラーバー表示 43 44 45plt.show()

試したこと

対応する基底関数系の像をcontourで図示したので、結果もcontourで図示しないといけないのがネックになっている。一応3次元データをプロットすれば求めているようなグラフになるのだが、contourfの場合、高度の値に二次元を求められるためにこのような冗長な式になっている。

補足情報(FW/ツールのバージョンなど)

ここにより詳細な情報を記載してください。

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

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

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

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

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

guest

回答2

0

ベストアンサー

27~34行を変える必要があると思われる。

fz 関数は以下の様になるかと思います。

python

1def fz(x,y): 2 z = np.zeros(x.shape) 3 for e in range(41): 4 if e % 2: 5 z += clf_i[e]*jv(e,x)*np.sin(e*y) 6 else: 7 z += clf_i[e]*jv(e,x)*np.cos(e*y) 8 return z

投稿2023/01/08 01:12

編集2023/01/08 01:15
melian

総合スコア19427

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

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

NMKN

2023/01/08 12:44

np.zeros(x.shape)は知りませんでした。ありがとうございます!参考にさせていただきます!
guest

0

ぱっと見た感じfz関数に問題があります。
ループの最後(e==40)のaとbだけがreturn a+bされています。(e=0〜38のaとbは使われていない)

投稿2023/01/07 23:02

TaroToyotomi

総合スコア1430

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

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

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

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

ただいまの回答率
85.54%

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

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

質問する

同じタグがついた質問を見る

CSV

CSV(Comma-Separated Values)はコンマで区切られた明白なテキスト値のリストです。もしくは、そのフォーマットでひとつ以上のリストを含むファイルを指します。

JupyterLab

JupyterLabは、Jupyter notebookの後継の対話型開発環境(IDE)です。データの可視化がインタラクティブで、プラグイン作成により新しいコンポーネントの追加および既存のコンポーネントも統合可能。サーバに閉じているため、データ分析に向いています。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。