これは、VBで作ってみたPICコードなのですが、このコードではSin波形が出力されます。この波形を変えたいのですが、どこを変えればいいのでしょう? どなたか分かる方教えてください。 (つづきがあります。)
'一次元静電粒子コード
M = 3:IM = 2 ^ (M + 1) + 1:IM2 = IM * 2 - 1
Dim RO(18):Dim PHI(18):Dim A(18, 4):Dim E(18):Dim X(161):Dim V(161)
' X:電子の位置(X座標) V:電子の速度 RO:電荷密度
' A(I,J):ポアソン・マトリックスエレメント PHI:静電ポテンシャル
' E:静電界(X方向) NP:1セル当りの粒子数 XMAX:全メッシュ数
' KK:最大計算回数 DT:タイムステップ XMAX:X方向の最大値点
' DX: メッシュ間隔
NP = 10:NPT = NP * (IM - 1):PW = 1 / NP
XMAX = IM - 1:DX = XMAX / (IM - 1):DXP = DX / NP:DT = 0.1
' 粒子の初期設定
For IP = 1 To NPT
X(IP) = DXP * (IP - 0.5):V(IP) = 0
Next IP
' ポテンシャルの初期値
For I = 1 To IM
PHI(I) = 5! * Sin(2 * 3.14159 * (I - 1) / (IM - 1))
Next I
' 計算開始
For KK = 1 To 120
' 電場の計算
For I = 1 To IM - 1
E(I) = -(PHI(I + 1) - PHI(I)) / DX
Next I
E(IM) = E(1):E(0) = E(IM - 1):E(IM + 1) = E(2)
' 粒子の運動
For IP = 1 To NPT
XNM = X(IP) / DX + 0.5:IXP = Int(XNM):IXM1 = IXP - 1
If IXM1 <= 0 Then IXM1 = IM - 1 - IXM1: B = XNM - IXP
If IXP = 0 Then IXP = IM - 1:EP = E(IXP) * B + E(IXM1) * (1 - B)
V(IP) = V(IP) + Q * EP * DT: X(IP) = X(IP) + V(IP) * DT
If X(IP) > XMAX Then X(IP) = X(IP) - XMAX
If X(IP) < 0 Then X(IP) = XMAX + X(IP)
Next IP
つづく・・・
A 回答 (1件)
- 最新から表示
- 回答順に表示
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
おすすめ情報
- ・漫画をレンタルでお得に読める!
- ・人生のプチ美学を教えてください!!
- ・10秒目をつむったら…
- ・あなたの習慣について教えてください!!
- ・牛、豚、鶏、どれか一つ食べられなくなるとしたら?
- ・【大喜利】【投稿~9/18】 おとぎ話『桃太郎』の知られざるエピソード
- ・街中で見かけて「グッときた人」の思い出
- ・「一気に最後まで読んだ」本、教えて下さい!
- ・幼稚園時代「何組」でしたか?
- ・激凹みから立ち直る方法
- ・1つだけ過去を変えられるとしたら?
- ・【あるあるbot連動企画】あるあるbotに投稿したけど採用されなかったあるある募集
- ・【あるあるbot連動企画】フォロワー20万人のアカウントであなたのあるあるを披露してみませんか?
- ・映画のエンドロール観る派?観ない派?
- ・海外旅行から帰ってきたら、まず何を食べる?
- ・誕生日にもらった意外なもの
- ・天使と悪魔選手権
- ・ちょっと先の未来クイズ第2問
- ・【大喜利】【投稿~9/7】 ロボットの住む世界で流行ってる罰ゲームとは?
- ・推しミネラルウォーターはありますか?
- ・都道府県穴埋めゲーム
- ・この人頭いいなと思ったエピソード
- ・準・究極の選択
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
フェルミエネルギー
-
量子力学的に光子のスピンは何...
-
Lennard-Jones pairwise potent...
-
固相焼結と液相焼結について
-
光の波動と粒子の二重性
-
シュレディンガーの猫(3)
-
ダイラタンシー現象と液状化現...
-
統計力学(量子統計) 演習問題
-
原子内の電子について.2
-
スピンヘリシティとスピンカイ...
-
低圧側から高圧側への気体の移動
-
酢酸ビニルの乳化重合について。
-
4成分でないDirac方程式
-
負のエネルギー、反粒子
-
量子論が意味不明すぎですw
-
基本粒子の大きさについて
-
素粒子
-
光は波ですか?それとも直進で...
-
波動関数の二乗は確率か確率密度か
-
化学ポテンシャルと粒子数─μが...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報