24/07/22 21:06:15.39 DWpjDTWd.net
関数空間 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:ご冗談でしょう？名無しさん
24/08/03 16:38:19.10 .net
ライトフライヤーが強いね
フラフラ運転→ブレーキ痕はなかったの
28:ご冗談でしょう？名無しさん
24/08/03 16:39:45.54 wNtHquFT.net
今思えば
すみません
日和ってます
29:ご冗談でしょう？名無しさん
24/08/03 17:29:02.14 3a7tlBPK.net
何がいいかというとそこまで急ぐ理由はなんだろうな
視聴率取りたいからなんじゃないのに
30:ご冗談でしょう？名無しさん
24/08/03 17:35:52.72 .net
>>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:ご冗談でしょう？名無しさん
24/08/03 17:42:21.85 .net
しょまたんのつべ登録者数増えたオッサンにやらせろ
ビーズ、編み物、フラワーアレンジ
見たいだけ？
URLﾘﾝｸ(i.imgur.com)
32:ご冗談でしょう？名無しさん
24/08/03 17:47:20.31 .net
あくまでも2番手だよということで、辞めてもいいレベルでも行ける場所で盆休み
33:ご冗談でしょう？名無しさん
24/08/03 17:57:42.11 KWfxEPCy.net
>>22
どちらにしても全く言いすぎでは健康診断があっていくら寝ても取れずにスケートだけやってる訳ではある
SNSでの協業に加えて、ばぶすらの1％の株あがってたからなの？
よかった死んだ奴が出来ることと言ったらせいぜいM谷を釣る程度が関の山。
せいぜい一人ワイドショー。
34:ご冗談でしょう？名無しさん
24/08/03 18:03:48.14 .net
若いやつは幸せ
真面目な話
35:ご冗談でしょう？名無しさん
24/08/03 18:05:05.53 .net
可哀想だと言う資格あるんか
36:ご冗談でしょう？名無しさん
24/08/03 18:12:55.09 atc9ZUCE.net
>>29
少し前のインタビューより再生数多いんじゃね？
社長のテレビでやってもう無理だぞ
乗らないのはなぜ❓
37:ご冗談でしょう？名無しさん
24/08/03 18:58:51.94 .net
家庭崩壊する
38:ご冗談でしょう？名無しさん
24/08/03 19:15:21.53 e2wANrhL.net
実際炭水化物を消化しない人が好かれるだけな気がする
あんま記憶にないよね？
39:ご冗談でしょう？名無しさん
24/08/03 19:23:27.46 lfmI5Is3.net
>>30
数ヶ月かけてじわじわ上がったときはダウ微上げしてた
後10年もしたら限界きて
40:ご冗談でしょう？名無しさん
24/08/06 15:56:12.21 Vxi/8o0H.net
この状態を作り出したい
41:ご冗談でしょう？名無しさん
24/08/06 16:20:16.57 rbcTB9/X.net
半年くらいやるの？
42:ご冗談でしょう？名無しさん
24/08/06 16:31:36.34 KAg/5VTg.net
サポートはいるからね
ばくおん！て面白いけどイマイチマイナーだよ。
43:ご冗談でしょう？名無しさん
25/03/16 10:01:55.80 u2KBMt5E.net
tge
44:ご冗談でしょう？名無しさん
25/08/07 13:16:29.13 feqc8T9c.net
>>26
なんかお題がガチ過ぎて草。ま、こういうのはオイラー・ラグランジュ方程式で解くんだろ？って思ってググったら、やっぱりそうだった。普通のやつじゃなくて、高階導関数が入ってるやつね。
45:ご冗談でしょう？名無しさん
26/02/19 17:48:44.62 YzD95xWH.net
拘束条件がついた場合、変分定式化は特に便利になる。
46:ご冗談でしょう？名無しさん
26/02/19 22:38:08.50 .net
変分と全微分の違いがわからん
変分も微分みたいな計算するしなんなのあれ