Skip to main content

[sympy] タグが付いた質問

sympy は python で、記号計算を行うためのライブラリです。多項式同士の演算を行ったり、解の公式で方程式を解いたりなどを行うことができます。

絞り込み
並べ替え
タグ付けされた
0
1 回答
241 閲覧数

sympyのsubsの使い方を教えて下さい。

from sympy import * var('RA,RAv,RB',real=True) rep ={RA:RAv} ;print("# ",rep) rep |={RB:RA.subs(rep)} ;print("# ",rep) # {RA: RAv} # {RA: RAv, RB: RA} ...
mrrclb48z's user avatar
  • 429
0
1 回答
80 閲覧数

sympyでsolveを2回実行すると、TypeErrorがでます。

原因を教えて下さい。よろしくお願いします。 from sympy import * A=Symbol('A',real=True,positive=True) A=solve(Eq(1,A**2),A) print("# 1",A) A=solve(Eq(1,A**2),A) print("# 2",A) # Traceback (most recent ...
mrrclb46z's user avatar
0
0 回答
102 閲覧数

sympyでの3次元グラフ表示のアイデアを求めています

sympyについて質問です。結果としては画像のような図を出力したいです。 用いる式は次のようなものです。 P = -A**2 + 2*A*D + 2*A*x*math.cos(phi_0) + 2*A*y*math.sin(phi_0) - B**2 + C**2 - D**2 - 2*D*x*math.cos(phi_0) - 2*D*y*math.sin(phi_0) - x**2 - y**...
tomi's user avatar
  • 1
0
0 回答
1,101 閲覧数

sympyについて質問です「TypeError: list indices must be integers or slices, not float」

画像の赤丸の部分を(Φ+2π/3)に変更したいです。 この部分を変更しようと思っています。 B_vectors = [A_vectors[i] + B*(unit_vectors[i] * cos(thetas[i])-ez*sin(thetas[i]) ) for i in range(3)] 単純に B*(unit_vectors[i+(2.0*np.pi/3.0)] ...
tomi's user avatar
  • 1
0
1 回答
475 閲覧数

Pythonのclassの作り方を教えて下さい。「...object has no attribute...」(3Dの三角形)

3次元の三角形です。よろしくお願いします。 (参考)SymPy 1.9 documentation » SymPy Modules Reference » Geometry » Plane https://docs.sympy.org/latest/modules/geometry/plane.html (参考)2.8.1 クラスの定義<Pythonのテキスト作りました - Qiita https:...
mrrclb48z's user avatar
  • 429
3
3 回答
530 閲覧数

この指数関数についての方程式はどのように解けばいいでしょうか

以下の指数関数のxについての方程式はどのようにして解けばいいでしょうか? exp(-1/x**2) = exp(-2/x**2) - exp(-3/x**2) sympyのsolveメソッドを使ったのですが実行が無限に行われうまくいきませんでした。環境はjupyter notebookです。
taro's user avatar
  • 31
-1
3 回答
162 閲覧数

四面体(三角錐)の内接球のsympyを探しています。 vbaソルバーでやってみました。

①四面体(三角錐)の内接球のsympyを探しています。他言語でも。 ②vbaソルバーで実行してみました。vbaコードを短くする方法を教えて下さい。一例しか実行していません。 ➂以下のコードの、「Range("A7").Formula」複数行で書けませんでした。 よろしくお願いします。 オリジナルポスト (0,0,0),(1,0,0),(0,1,0),(0,0,2) 「...
mrrclb48z's user avatar
  • 429
1
1 回答
302 閲覧数

sympyのlistの使い方(matrix-->list)を教えて下さい。 3点の平面の方程式

mylist1,mylist2,mylist3を使わない方法を教えて下さい。関数だと助かります。 参考 3点 A(1,1,2),B(0,-2,1),C(3,-1,0) を通る平面の方程式を計算せよ。 https://manabitimes.jp/math/679 http://meg.aalip.jp/python/math01.html https://atatat.hatenablog.com/...
mrrclb48z's user avatar
  • 429
0
1 回答
232 閲覧数

sympyで「AttributeError: 'str' object has no attribute 'subs'」がでました。

文字列'x-y'をx-yに変換する方法を教えて下さい。 from sympy import * var('x y') print(type( x-y ),( x-y ).subs(x,0.0)) print(type('x-y')) print(type('x-y'),('x-y').subs(x,0.0)) # <class 'sympy.core.add.Add'> -y # &...
mrrclb46z's user avatar
-1
1 回答
962 閲覧数

2本の直線の交点計算をsolveでやってみたい。

直線が、x軸又はy軸に平行の時、場合分けの対応が必要ですか?何ケースに分ければいいですか? よろしくお願いします。 以下のエラーが、 例7ででます。 UnboundLocalError: local variable 'a' referenced before assignment (オリジナルポスト)2本の線分が交差しているか調べる https://qiita.com/tydesign/...
mrrclb48z's user avatar
  • 429
0
1 回答
433 閲覧数

sympyでlistとMatrixの相互変換をやってみました。Matrix→list変換を教えて下さい。

①Matrix→list変換を教えて下さい。 ②Matrix→list変換。list→Matrix変換。それぞれ1行で書けますか? よろしくお願いします。 (参考) sympyのMatrix.row_insertの使い方を教えて下さい。 from sympy import * mylist=[[0]*3]*3 mylist[0]=[-4, 1, 9] mylist[1]=[1, 2, -4] ...
mrrclb48z's user avatar
  • 429
-2
1 回答
181 閲覧数

sympyのMatrix.row_insertの使い方を教えて下さい。

以下のエラーがでました。 AttributeError: 'list' object has no attribute 'cols' 以下のように出力する方法を教えて下さい。 Matrix([[-4, 1, 9], [1, 2, -4], [2, -3, 6]]) from sympy import * var('myMatrix') mylist=[[0]*3]*3e mylist[0]=[-4,...
mrrclb48z's user avatar
  • 429
-1
1 回答
395 閲覧数

sympyで,連立1次不等式の係数を取り出してみました。

①プログラムが短くなりますか? ②x,y,zにするためのアドバイスお願いします。 4x+y <= 9,x+2y >= 4,2x-3y >= -6 ----> Ax>=b (参考) sympyのcoeffの使い方を教えて下さい from sympy import * var('x y') def my_Str2Inequality(my_Str): if ...
mrrclb48z's user avatar
  • 429
0
2 回答
407 閲覧数

sympyで行列式の結果を分数でだす方法を教えて下さい。

①sympyで行列式の結果を分数56/5でだす方法を教えて下さい。 ②x−y+2=0,3x+y−18=0,2x+3*y−12=0から、係数を取り出して、リストl1l2l3をつくる方法を教えて下さい。 ③プログラムが短くなりますか? sympyで「3直線で作られる三角形の面積計算」やってみた (抜粋)【数学ⅡB】3直線で作られる三角形の面積【明治学院大】 x−y+2=0,3x+y−18=0,2x+...
mrrclb48z's user avatar
  • 429
0
3 回答
243 閲覧数

素因数分解を約数にする。ループ,組み合わせの方法を教えて下さい。

問題点① 素因数分解の個数が3個の場合しか対応していません。 問題点② リスト内包表記の使い方がわかりませんでした。関数部分1行で書けますか。 問題点③ ( . [ . { 。使い方、変換についての、おすすめのページを教えて下さい。 よろしくお願いします。 from sympy import * import itertools var('My_Yakusu') def ...
mrrclb48z's user avatar
  • 429
0
1 回答
128 閲覧数

sympyのUnionの中のInterval範囲の個数は、いくつになりますか?

①sympyのUnionの中のInterval範囲の個数は、いくつになりますか? ②これでいいですか。「Union→Interval→Union」をしてみました。 よろしくお願いします。 from sympy import * k = Symbol('k') ineq = 'k**2-40>0' ans1=solve_univariate_inequality(sympify(ineq), ...
mrrclb48z's user avatar
  • 429
0
3 回答
482 閲覧数

sympyの1変数不等式「Inequality Solvers」について教えて下さい。

式 :k**2-40>0 条件 :kは整数,k>0 ほしい解:k>=7 プログラム①:symbol定義は、だめでした。 プログラム②:Unionの要素の取り出し方を教えて下さい。 プログラム③:solveの結果の要素の取り出し方を教えて下さい。 よろしくお願いします プログラム① from sympy import * k = Symbol('k',integer=...
mrrclb48z's user avatar
  • 429
-2
2 回答
430 閲覧数

sympyの代入方法を教えて下さい。

私の代入のどこが、おかしいですか? myV0に整数しか入らないと思っていました。 (参考) python のitertools.permutationsの使い方を教えて下さい。全部の組み合わせの作り方を教えて下さい。 from sympy import * var('V0 V1 myV0 myV1') my_list=list(range(2)) print(my_list) for v0 in ...
mrrclb48z's user avatar
  • 429
0
1 回答
1,392 閲覧数

sympyで、円周率(piとmath.pi)の使い方の違いを教えて下さい。

①計算1と計算2の使い方の違いを教えて下さい。 ②キーインを少なくしたいので、piだけでできると助かります。 ③(追加)ルートの正解は、sympyで、できるはず? sympyで入試問題に挑戦! 解説は、理解できてません。方べきの定理ですか? (参考)入試問題に挑戦! 平成17年宮城県(超難) http://www.5kaku.net/?p=9299 (参考)WolframAlphaで 28/3 π ...
mrrclb48z's user avatar
  • 429
0
1 回答
120 閲覧数

sympy3Dで、solveを使わないかっこいい方法を教えて下さい。

法線ベクトル? 中点?対称点の座標を求める【京大過去問】 (参考)法線ベクトルの求め方と空間図形への応用>記事の最後には、実際に京都大学で出題された「空間に存在する平面に対称な点」を求める問題を掲載しました。 https://linky-juku.com/normal-vector/ (参考)2006年度 京都大・文系数学 問2 https://www.densu.jp/kyoto/...
mrrclb48z's user avatar
  • 429
0
0 回答
133 閲覧数

sympyの積分で、subsで後から代入の方法を教えて下さい。

subsで後から代入。計算2で、pi/2にならなくて、x変数が残ります。いい方法がありますか? from sympy import * var("x r") f=sqrt(r**2-x**2) print("計算1:r=r 面積=",integrate(f,(x,-r,r))) print("計算2:r=1 面積=",integrate(...
mrrclb48z's user avatar
  • 429
0
2 回答
139 閲覧数

sympy に平方数を判定する is_square がありますが、factorint で一行でできますか?

短くなりますか? よろしくお願いします。 def My_is_square(x): myDic=factorint(x) i=0 for k, v in myDic.items(): if v % 2 == 0: i=i+1 if i == len(myDic): return True else: ...
mrrclb48z's user avatar
  • 429
0
2 回答
423 閲覧数

sympyで三角関数の値を分数表示したい

(参考)wolframalphaなら、できました。 入力: cos(x)=3/5,sin(x) https://www.wolframalpha.com/input/?i=cos%28x%29%3D3%2F5%2Csin%28x%29&lang=ja 代入: sin(x) = -4/5, sin(x) = 4/5 # 共通テスト2021 数学ⅠA問題 第1問〔2〕 # https://...
mrrclb48z's user avatar
  • 429
0
1 回答
161 閲覧数

平方数の判定を、sympyの素因数分解factorintでできますか?

目的は、平方数の判定を一行でできますか?よろしくお願いします。 print(factorint(9)) {3: 2} print(factorint(81)) {3: 4} 参考nが平方数かどうか調べる時 https://teratail.com/questions/200122 参考:Wolfram|Alpha 81は平方数ですか https://ja.wolframalpha.com/...
mrrclb48z's user avatar
  • 429
0
1 回答
123 閲覧数

SymPy Liveの和集合unionの正しい使い方を教えて下さい。

https://live.sympy.org/ 入力1 print( FiniteSet("A", "D", "A", "G", "B", "C", "B").union(FiniteSet())) 出力1 FiniteSet(A, B, C, D, G) ...
mrrclb48z's user avatar
  • 429
0
1 回答
224 閲覧数

sympyで、分配して約分する方法を教えてください

Jを変数として、以下を変形したいです。ゴールは、分数式ではなく正式になることです。2*(J + 18)+ 2とかになれば成功です。宜しくお願いします。 (2 + 2/(J + 18))*(J + 18) 追記 すみません、本当は行列の成分で分配を行いたいのです。コードは以下です。最終的には、行列の成分を正式だけにしたいのです。 B=sp.Matrix([[0],[0],[-e*(M*r*H/(z*...
SpaceTAKA's user avatar
  • 267
0
0 回答
144 閲覧数

行列の連立方程式を解きたいです。

参考までに、これはリカッチの代数方程式で求めた最適レギュレータのゲインについてのお話です。ですが、周辺知識がなくとも問題ありませんので読み進めていただければ幸いです。私はsympyでやりましたが、ライブラリはなんでもいいです。 求めたいのは、変数Jです。Jを求めるには、以下の二式の行列の連立方程式を解いて、成分を明らかにする必要があります。 K-inv_R@tB@P=0 A@tP+P@A-P@B@...
SpaceTAKA's user avatar
  • 267
0
0 回答
534 閲覧数

Jupyter Lab で sympy.init_session() を用いるとそれ以降のセルが動かなくなる

Sympyをimportしinit_session() 又は init_printing()を用いるとそれ以後のセルが[*]となってしまいどれだけ待っても実行されません。 Jupyter notebookを使用しても同様の問題が生じます。 コンソール上のIpython上ではinit_session()を行っても正常に動作します。 Python 3.8.5 (default, Sep 3 ...
Yuz's user avatar
  • 1
0
0 回答
112 閲覧数

symbolを入力とすることは出来ない?

以下のコードにおいて、inputにtの関数を入れれば動きますので、”symbolを入力とすることは出来る”と言えそうなのですが、typeで型を見ればわかる通り、これはstr扱いです。シンボルとして入力したいのですがどうすれば良いでしょうか import sympy as sp t=sp.symbols("t") f=input() print("f",...
SpaceTAKA's user avatar
  • 267
0
1 回答
534 閲覧数

sympyのsympy.plotはプロット出来ない?NameError: name 'Integral' is not defined

以下のコードでは行列計算を行っており、最終的に行列 X_ の成分を求めています。 成分はシンボル t の関数で、続くコードではその関数のグラフを描き出そうとしています。 しかし、どういうわけか "NameError: name 'Integral' is not defined" というエラーが出てきます。 理由として思い当たるのは X_ を表示すればわかる通り、...
SpaceTAKA's user avatar
  • 267
0
0 回答
259 閲覧数

VScodeでMicroPyhtonとSympyは同時に使えない?

タイトルがわかりにくくてすみません。プログラミング初心者なので、どこが分からないか分からないに近い状態なのですが、ひとつひとつ説明していきますので辛抱強く読んでくださると嬉しいです。 バックグラウンド 私はEV3というLEGO社が出しているロボット組み立ての知育玩具で遊んでいます。これがなかなか本格的で、色んな言語でプログラミング出来ます。メソッドやクラスは先に設定されていて、...
SpaceTAKA's user avatar
  • 267
0
1 回答
355 閲覧数

sympy の関数をmatplotlibでプロットしたい

以下のようなコードを書きました。そして、X[0]は関数なので、plt.plot(X[0](t),t)のようにプロットしたいのですが、どうやらsympyのオブジェクトは呼び出せないっぽいです。どうすれば良いでしょうか? import sympy as sp import matplotlib.pyplot as plt A = sp.Matrix([[-1, 0], [0, -1]]) ...
SpaceTAKA's user avatar
  • 267
1
2 回答
223 閲覧数

行列指数関数の極限をとりたい

このコードによって表示される行列指数関数を、tを無限に飛ばしたときにどんな行列になるのかを表示させたいです。sym.limit() は試しましたが、要素を飛ばすことはできないようです。 import sympy as sym A = sym.Matrix([ [1, 0], [0,1]]) t=sym.symbols("t") At=A*t At.exp()
SpaceTAKA's user avatar
  • 267
-1
2 回答
305 閲覧数

定義した関数を表示したい

例えば以下の様に記述しても x**2 と表示される訳ではありません。しかし僕は表示させたいです。自分で定義しても関数の形がわからないことだってあるじゃないですか?それを解決したいのです。 def f(x): return x**2 print(f(x)) 以下のような場合、返値は定義できても関数形は不明です。それを表示させたいのです。 よろしくお願いします。 A=np.matrix(...
SpaceTAKA's user avatar
  • 267
0
1 回答
176 閲覧数

リストのすべての要素を表現したい

例えば以下のような条件文が書きたいときに、どのようにコードすれば良いのでしょうか? suuji=[0,-2,3,0,5,3,7] if suuji[すべて]>0: ー補足ー 皆様回答ありがとうございます。とても学びになりましたが、目的は達成されませんでした。ですのでストレートに目的を書きます。 numpyのsolveメソッドによって多項式の解が分かります。...
SpaceTAKA's user avatar
  • 267
1
0 回答
119 閲覧数

sympy で cos(1/4π) は、√2/2 と展開されますが、cos(3/4π) は、cos(0.75π) と展開されます。 cos(3/4π) も同様に展開できないでしょうか。

sympy の version は、1.6.2 python の version は、3.7.4 OS は、windows10 64bit です。 下記、試した結果です。trigsimp が有ってもなくても、同じでした。 sympy.trigsimp(cos(0/4*sympy.pi)) -> 1 sympy.trigsimp(cos(1/4*sympy.pi)) -> √2/2 ...
uma3san's user avatar
  • 11
1
1 回答
132 閲覧数

このコードはP制御になっていますでしょうか?

Mは操作量を表します。 添字の1は一回前の操作であることを表します。 これを実行すると比例ゲインKp=0.8の時にオフセットなしでゴールまで達するので違和感です。比例制御はその仕組みのせいで必ずオフセットを生じるはずです… import numpy as np from matplotlib import pyplot as plt from numpy.random import * ...
SpaceTAKA's user avatar
  • 267
0
1 回答
790 閲覧数

倒立振子(台車型)のPID制御について

バックグラウンド 私は今、LEGO EV3というプログラミング練習のオモチャで遊んでいます。数日前から倒立振子の作成に取り掛かりました。私は制御工学初心者ですので、現代制御はまだ難しく感じ、コーディングが簡単そうな古典制御のPID制御で作ることにしました。 制御もプログラミングも初心者なので、とりあえず必要と思う情報を書きます。不足があれば教えてください。 Python3.8,...
SpaceTAKA's user avatar
  • 267
0
1 回答
536 閲覧数

PID制御のコーディングについて

Pythonで制御工学を学んでいます。そこで、このコードに疑問点があります。 サンプル方式 MVn = MVn-1 + ΔMVn ΔMVn = Kp(en-en-1) + Ki en + Kd((en-en-1) - (en-1-en-2)) #ただし MVn、MVn-1:今回、前回操作量  ΔMVn:今回操作量差分 en,en-1,en-2:今回、前回、前々回の偏差 とする。 ...
SpaceTAKA's user avatar
  • 267
0
1 回答
4,249 閲覧数

pythonのエラーについて TypeError: ufunc 'kv' not supported

Pythonに慣れておらず、エラー内容と解決方法がわからない。 これは呼び出したkvが積分に対応していないということですか? ソースコード import numpy as np import pandas as pd import matplotlib.pyplot as plt import scipy.special from scipy.special import kv from ...
asak's user avatar
  • 1
3
1 回答
630 閲覧数

ナブラの内積を計算したいです

∇(ナブラ)を使って内積 ∇・b の計算をしたいのですが、ここからどうすればうまくいくのでしょうか。3を出力したいです。 import numpy as np x = Symbol('x') y = Symbol('y') z = Symbol('z') del_x = np.gradient(,dx) del_y = np.gradient(,dy) del_z = np.gradient(,...
ara's user avatar
  • 103
4
1 回答
326 閲覧数

sympy.printing.mathmlで括弧()の扱いについて

以下のように、Pythonの文字列をsimpifyで数式に変換し、 mathmlコマンドで数式をMathML形式に変換すると from sympy import * from sympy.printing.mathml import mathml print(mathml(sympify("(2*x + 3*y + 1)*(4*x**2 - 6*x*y - 2*x + 9*y**2 - 3*y ...
pekochun's user avatar
  • 417
1
1 回答
201 閲覧数

sympifyで大文字の"S"を数式に変換するには (Python Sympy)

文字列から、MathMLのコードを生成するプログラムtest001.py を書いています。 test001.py import sys sys.path.append('/home/vagrant/.local/lib/python3.4/site-packages') from sympy import * from sympy.printing.mathml import mathml ...
pekochun's user avatar
  • 417
3
1 回答
407 閲覧数

Python3.6にアップデートしたら、Python2.7で使えてたものが動作しなくなった

Pythonのsympyモジュールのsympy.printing.mathmlを使用しています。 Python3.6で、以下のtest1.pyを実行すると、エラーが起きてしまいます。 test1.py import sys, json from sympy import * from sympy.printing.mathml import mathml print(mathml(1/6))...
pekochun's user avatar
  • 417
1
1 回答
677 閲覧数

Pythonでsortした際のエラーについて

以下の、test.pyというファイルを実行したところ、エラーが出ます。 sort()の行でエラーが出ているようですが、原因は何でしょうか? また、配列内の要素に虚数が含まれている条件でsortする方法はありますでしょうか? test.py import sys, json sys.path.append('/home/vagrant/.local/lib/python3.4/...
pekochun's user avatar
  • 417
0
1 回答
407 閲覧数

sympyのupgradeで旧versionがuninstallできない?

sympyを1.1.1から1.3に pip install --upgrade sympy でupgradeしようとしたら、その時に古いVersionのSympyがimportされていたためか、 Installing collected packages: sympy Found existing installation: sympy 1.1.1 Cannot uninstall '...
yton's user avatar
  • 3
1
0 回答
428 閲覧数

python,sympyでinit_printing(order='rev-lex')が効かない

ソースコード from sympy import Symbol, pprint, init_printing def print_series(n): # 次数の低い順から並べる init_printing(order='rev-lex') x = Symbol('x') series = x for i in range(2, n+1): ...
user31543's user avatar
0
1 回答
98 閲覧数

sympyのplotのgcaのAttributeErrorを教えて下さい。

細長く表示します。 (参考)大学入試問題で現れる世界遺産 https://qiita.com/kazukiii/items/bbb55f873730b6e31f2b from sympy import * import matplotlib.pyplot as plt var('x') F = x**4 - x**2 + 6 G = cos(2 * pi * x) / 2 + 7 / 2 H =...
mrrclb48z's user avatar
  • 429
0
1 回答
274 閲覧数

SymPyで作成したグラフの一部描画について

SymPyで作成した数式のグラフについて質問致します。 以下のサンプルコードを例にしますと、座標平面上の第一象限に一次関数と円のグラ フが計三つ作成されます。 これらのグラフにより分割される各領域に関しまして、例えば三つのグラフの&を取 る領域等、特定の領域のみを描画することは可能でしょうか。 もしお分かりでしたら、ご助言よろしくお願い致します。 from sympy import * x,...
poko's user avatar
  • 77
3
2 回答
457 閲覧数

SymPy と SciPy におけるグラフの統合

以下のサンプルコードにて、5つの座標点からなるボロノイ図をscipy.spatial.voronoi_plot_2dを使って生成します。 このグラフ上にsympy.plotting.plot_implicitを使って任意の陰関数(双曲線等)を描画したいのですが、やり方が分からず困っております。 適切な方法をご存知でしたら、ご教授願います。 import numpy as np import ...
poko's user avatar
  • 77