変分法って何にでも使えるな (44ﾚｽ)
上下前次1-新
1: ご冗談でしょう？名無しさん [] 2019/07/14(日)13:40 ID:QhzfCf1z(1)
なんだこのクッソ便利なもの
25: ご冗談でしょう？名無しさん [] 2024/07/22(月)21:05 ID:DWpjDTWd(2/3)
質点の運動を記述するLagrangianが与えられたとき、Hamiltonの原理を用いて運動方程式を導け。
26(2): ご冗談でしょう？名無しさん [] 2024/07/22(月)21:06 ID:DWpjDTWd(3/3)
関数空間 C2[0,1] 上の汎関数
$$
J[y] = \int_0^1 \left(y''(x)^2 + y(x)^2 - 2xy(x)\right) dx
$$
の極値問題を解け。境界条件は y(0)=y(1)=0 とする。
27: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)16:38 ID:???
ライトフライヤーが強いね
フラフラ運転→ブレーキ痕はなかったの
28: ご冗談でしょう？名無しさん [] 2024/08/03(土)16:39 ID:wNtHquFT(1)
今思えば
すみません
日和ってます
29(1): ご冗談でしょう？名無しさん [] 2024/08/03(土)17:29 ID:3a7tlBPK(1)
何がいいかというとそこまで急ぐ理由はなんだろうな
視聴率取りたいからなんじゃないのに
30(1): ご冗談でしょう？名無しさん [sage] 2024/08/03(土)17:35 ID:???
>>26
import numpy as np
from scipy.optimize import minimize
def J(y, bc_type=None):
"""
Functional to be minimized.
Args:
y: Array representing the function y(x).
bc_type: Optional string specifying the boundary condition type.
- 'dirichlet': Enforces y(0) = y(1) = 0 (default).
- None: No boundary condition enforced.
Returns:
The value of the functional J evaluated at y.
"""
y_prime = np.gradient(y)
y_doubleprime = np.gradient(y_prime)
return np.trapz(y_doubleprime**2 + y**2 - 2*y*y_prime, np.linspace(0, 1, len(y)))
def bc_dirichlet(y):
"""
Constraint function for Dirichlet boundary conditions.
Args:
y: Array representing the function y(x).
Returns:
A list containing the constraint values (y(0) and y(1)).
"""
return [y[0], y[-1]]
# Define the problem with Dirichlet boundary conditions
bounds = ((0, None),) * len(np.linspace(0, 1, 101)) # No lower bound for y
cons = ({'type': 'eq', 'fun': bc_dirichlet})
# Solve the minimization problem
initial_guess = np.random.rand(101) # Random initial guess
sol = minimize(J, initial_guess, method='SLSQP', bounds=bounds, constraints=cons)
# Extract the solution
y_opt = sol.x
# Print the minimum value of the functional
print("Minimum value of J:", J(y_opt))
# Plot the solution (optional)
import matplotlib.pyplot as plt
x = np.linspace(0, 1, len(y_opt))
plt.plot(x, y_opt)
plt.xlabel("x")
plt.ylabel("y(x)")
plt.title("Solution of the minimization problem")
plt.show()
31: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)17:42 ID:???
しょまたんのつべ登録者数増えたオッサンにやらせろ
ビーズ、編み物、フラワーアレンジ
見たいだけ？
https://i.imgur.com/hjm9uGV.jpeg
32: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)17:47 ID:???
あくまでも2番手だよということで、辞めてもいいレベルでも行ける場所で盆休み
33: ご冗談でしょう？名無しさん [] 2024/08/03(土)17:57 ID:KWfxEPCy(1)
>>22
どちらにしても全く言いすぎでは健康診断があっていくら寝ても取れずにスケートだけやってる訳ではある
SNSでの協業に加えて、ばぶすらの1％の株あがってたからなの？
よかった死んだ奴が出来ることと言ったらせいぜいM谷を釣る程度が関の山。
せいぜい一人ワイドショー。
34: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)18:03 ID:???
若いやつは幸せ
真面目な話
35: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)18:05 ID:???
可哀想だと言う資格あるんか
36: ご冗談でしょう？名無しさん [] 2024/08/03(土)18:12 ID:atc9ZUCE(1)
>>29
少し前のインタビューより再生数多いんじゃね？
社長のテレビでやってもう無理だぞ
乗らないのはなぜ❓
37: ご冗談でしょう？名無しさん [sage] 2024/08/03(土)18:58 ID:???
家庭崩壊する
38: ご冗談でしょう？名無しさん [] 2024/08/03(土)19:15 ID:e2wANrhL(1)
実際炭水化物を消化しない人が好かれるだけな気がする
あんま記憶にないよね？
39: ご冗談でしょう？名無しさん [] 2024/08/03(土)19:23 ID:lfmI5Is3(1)
>>30
数ヶ月かけてじわじわ上がったときはダウ微上げしてた
後10年もしたら限界きて
40: ご冗談でしょう？名無しさん [] 2024/08/06(火)15:56 ID:Vxi/8o0H(1)
この状態を作り出したい
41: ご冗談でしょう？名無しさん [] 2024/08/06(火)16:20 ID:rbcTB9/X(1)
半年くらいやるの？
42: ご冗談でしょう？名無しさん [] 2024/08/06(火)16:31 ID:KAg/5VTg(1)
サポートはいるからね
ばくおん！て面白いけどイマイチマイナーだよ。
43: ご冗談でしょう？名無しさん [] 03/16(日)10:01 ID:u2KBMt5E(1)
tge
44: ご冗談でしょう？名無しさん [] 08/07(木)13:16 ID:feqc8T9c(1)
>>26
なんかお題がガチ過ぎて草。ま、こういうのはオイラー・ラグランジュ方程式で解くんだろ？って思ってググったら、やっぱりそうだった。普通のやつじゃなくて、高階導関数が入ってるやつね。
上下前次1-新書関写板覧索設栞歴
ｽﾚ情報 赤ﾚｽ抽出 画像ﾚｽ抽出 歴の未読ｽﾚ AAｻﾑﾈｲﾙ
ぬこの手 ぬこTOP 0.024s