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

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

ただいまの
回答率

90.84%

  • Python

    6341questions

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

  • Python 3.x

    4823questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

量子ウォークの別の書き方にチャレンジしており、少してこずっております。。。

解決済

回答 1

投稿 編集

  • 評価
  • クリップ 0
  • VIEW 125

jackson5

score 29

1次元量子ウォークの別の書き方をしています。
設定】
イメージ説明
上のスピンを[1,0],下のスピンは[0,1]とします。まず初期状態[1,0]の内部スピンをquantumcoinといわれるもので変更します。
その計算した内部スピンを上図のように上ならx+1,下はx-1に移動させます。この状態を時間発展させていきたいのですが、少しエラーがでてしまい一日奮闘したのですが、あと一歩のところで、上手くいきません。

import numpy as np
import matplotlib.pyplot as plt
import math

n=4
m=3
theta=3*(math.pi)/12

phi_list=[]
init_phi=[]
new_phi=[]

for i in range(0,2*n+1):
    a=[0,0]
    new_phi.append(a)
    phi_list.append(a)



#内部自由度をかえるコイン↓
def quantumcoin(theta):
    C = [[np.cos(theta),-np.sin(theta)],[np.sin(theta),np.cos(theta)]]
    return C
#print(quantumcoin(3*(math.pi)/12))

#t=0の時のphi状態
def initPositionPhi():
    phi = [1,0]
    othersphi=[0,0]
    for i in range(0,2*n+1):
        if i == n:
            init_phi.append(phi)
        else:
            init_phi.append(othersphi)
    return init_phi
#print(initPositionPhi())

def CoinOperator(x,new_phi):  #内部自由度を変える処理
    new_phi[x]=np.dot(quantumcoin(theta),new_phi[x])
    return new_phi[x]
print(CoinOperator(n,initPositionPhi()))


#def ShiftOperator(steps,phi):  
for t in range(1,m+1):      #内部自由度を変えた物に合わせて座標を移動(x=nがスタート)
    #if t==0: #怪しい
        #initPositionPhi()
    if t==1:
        new_phi[n-t][1] = (CoinOperator(n,initPositionPhi()))[1]
        new_phi[n+t][0] = (CoinOperator(n,initPositionPhi()))[0]
        print(t,new_phi)
    else:
        new_phi[n-t] = (CoinOperator(n-(t-1),new_phi))[1] + (CoinOperator(n-(t+1),new_phi))[0]
        new_phi[n+t] = (CoinOperator(n+(t-1),new_phi))[0] + (CoinOperator(n+(t+1),new_phi))[1]
    #print(t,new_phi)


結果は

[0.70710678 0.70710678]
1 [[0, 0], [0, 0], [0, 0], [0, 1.0000000000000002], [0, 0], [-0.7071067811865477, 0], [0, 0], [0, 0], [0, 0]]
1 [[0, 0], [0, 0], [0, 0], [0, 1.0000000000000002], [0, 0], [-0.7071067811865477, 0], [0, 0], [0, 0], [0, 0]]
2 [[0, 0], array([0., 0.]), 0.7071067811865477, array([-0.70710678,  0.70710678]), [0, 0], array([-0.5, -0.5]), -0.5000000000000001, array([0., 0.]), [0, 0]]
3 [array([0., 0.]), array([0.5, 0.5]), array([[ 0.5, -0.5],
       [ 0.5,  0.5]]), array([-0.70710678,  0.70710678]), [0, 0], array([-0.5, -0.5]), array([[-0.35355339,  0.35355339],
       [-0.35355339, -0.35355339]]), array([-0.35355339,  0.35355339]), array([0., 0.])]


と、結果が少しおかしくなってしまいます。
というのも、内部状態を変えた後の下スピンが反映されません。
式の定義が少しおかしいのかもしれませんが、解決できませんでした。
ご指摘、アドバイスよろしくお願いします。
以下のような答えになります。
R=[1 0],L=[0,1]としています。よろしくお願いします。
イメージ説明

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

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • hayataka2049

    2018/06/12 20:33 編集

    回答する人的には、理想の出力的なものを示していただけると多少協力しやすくなるかもしれません

    キャンセル

  • jackson5

    2018/06/12 20:34

    わかりました、少々お待ちください。

    キャンセル

回答 1

checkベストアンサー

+2

位置の移動を行なう時にnew_phiに直接計算結果をいれず,例えばtemp_phiなどのリストに計算結果を入れておき,全ての座標での計算が終わった後にnew_phiを更新する方がバグが少なくなる気がします.
また,[0,0]のリストをつくる時にa=[0,0]を代入してますが,aを定義せず[0,0]を直接入れた方が個人的には安心です(私はリストのコピー時の振る舞いをきちんと追えてないので).

座標を辞書のkeyとして管理するverを書いてみました.辞書であれば新しい座標が増えた時も対応が楽なので,速度が気にならないなら試してみてください.
問題設定を誤解していたら申し訳ありません.
参考になれば幸いです(根本的に間違っていたら申し訳ありません).

 コード

import math as m
import numpy as np

nstep = 4
theta = m.pi/4
# States {step:{pos:[R,L]}}
QStates = {0:{0:[1,0]}} # initial state

# Spinを回転
def quantumcoin(theta):
    C = [[np.cos(theta),-np.sin(theta)],
        [np.sin(theta),np.cos(theta)]]
    return C

# 状態を確認し,移動させる.
def movState(temp_state,dotted_state):
    pos = dotted_state.keys()
    for s in sorted(pos):
        for i, pos_ in enumerate([s+1,s-1]):
            if pos_ not in temp_state: temp_state[pos_]=[0,0]
            temp_state[pos_][i] += dotted_state[s][i]
    return temp_state


if __name__=="__main__":

    # main loop
    for s in range(nstep):
        temp = {}
        coined = {i:np.dot(quantumcoin(theta),QStates[s][i]) for i in sorted(QStates[s].keys())}
        QStates[s+1] = movState(temp,coined)

    # 結果を出力
    print("Step|"+"|".join([str(i) for i in range(-nstep,nstep+1)]))
    print("---|"+"|".join(["---" for i in range(-nstep,nstep+1)]))
    for s in sorted(QStates.keys()):
        text = str(s)
        for i in range(-nstep,nstep+1):
            if i in QStates[s]:
                text += "|({0:.3f}, {1:.3f})".format(QStates[s][i][0],QStates[s][i][1])
            else:
                text += "|(0,0)"
        print(text)

 出力結果

Step|-4|-3|-2|-1|0|1|2|3|4
---|---|---|---|---|---|---|---|---|---
0|(0,0)|(0,0)|(0,0)|(0,0)|(1.000, 0.000)|(0,0)|(0,0)|(0,0)|(0,0)
1|(0,0)|(0,0)|(0,0)|(0.000, 0.707)|(0,0)|(0.707, 0.000)|(0,0)|(0,0)|(0,0)
2|(0,0)|(0,0)|(0.000, 0.500)|(0,0)|(-0.500, 0.500)|(0,0)|(0.500, 0.000)|(0,0)|(0,0)
3|(0,0)|(0.000, 0.354)|(0,0)|(-0.354, 0.000)|(0,0)|(-0.707, 0.354)|(0,0)|(0.354, 0.000)|(0,0)
4|(0.000, 0.250)|(0,0)|(-0.250, -0.250)|(0,0)|(-0.250, -0.250)|(0,0)|(-0.750, 0.250)|(0,0)|(0.250, 0.000)

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/06/14 14:48 編集

    ありがとうございます,少しコードが難しいですね(~~'')。
    頑張ってみます

    キャンセル

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

  • ただいまの回答率 90.84%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

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

  • Python

    6341questions

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

  • Python 3.x

    4823questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。