pythonのソースコードの改変案

Question

前提・実現したいこと

下のソースコードの横軸の数値とランダムで出てきている点を自分で指定した数値にしたいのですがどのように直したらいいのかわかりません。

該当のソースコード

import numpy as np
import matplotlib.pyplot as plt

# hyperparameters
mu = 0.1
alpha = 1.5
beta = 30

# 
neuton_iter_num = 200
num_datapoints = 15

# step1
U = np.random.uniform(0, 1)
# step2
t1 = -np.log(U)/mu
t = [t1]


for k in range(num_datapoints):
    # step3
    U = np.random.uniform(0, 1)
    
    # step 4
    if k == 0:
        Sk = 1

    u = t[k] - np.log(U)/mu  # initial value of iteration
    for _ in range(neuton_iter_num):
        fu  = np.log(U) + mu*(u-t[k])  + alpha/beta*Sk*(1-np.exp(-beta*(u-t[k])))
        fu_dash = mu+alpha*Sk*np.exp(-beta*(u-t[k]))
        u -= fu/fu_dash

    # step 5
    t.append(u)
    Sk = np.exp(-beta*(t[k+1]-t[k]))*Sk + 1  # find S(k+1)

t = np.array(t)/max(t)

# plt.figure(figsize=(20, 10), dpi=50)
# plt.scatter(t, np.ones(len(t)), facecolors='none', edgecolors='red')
# plt.show()


# intensity
time = np.linspace(0,1,10000)
increase = list(map(lambda x: alpha*np.exp(-beta * x), time))
intensity = np.zeros(10000*2)
for timepoints in t:
    yuko = np.append(np.append(np.zeros(len(time[time<timepoints])), increase), np.zeros(len(time[time>=timepoints])))    
    intensity += yuko


# plt.figure(figsize=(20, 10), dpi=50)
# plt.plot(time, intensity[:10000])
# plt.show()



fig, ax1 = plt.subplots(figsize=(20, 5))
ax2 = ax1.twinx()
ax1.scatter(t, np.ones(len(t)), facecolors='none', edgecolors='red', marker='*', s=150)
ax2.plot(time, intensity[:10000], c='b', lw=1)
plt.show()

補足情報（FW/ツールのバージョンなど）

上のソースコードはこのサイトから取ってきました
https://babaye.hatenablog.com/entry/2019/10/17/214442

Comments

[fumu7]

「ランダムで出てきている点を自分で指定した数値にしたい」と書かれていますが、そもそもランダムである数値は一定ではないので、指定した数値に固定できるものではありません。
まずは、「ランダム」という言葉の意味を辞書で調べましょう！
意味も判らない用語に基づいて、まともなプログラムを作るのは無茶苦茶です。

[izm2192]

ランダムでだすというのをなくしてそこを自分が指定した変数にしたいという意味です、

Answer 1

U = np.random.uniform(0, 1)

が二か所にありますね。これを

S = input('0と1の間の数を入れてください.')

に変更しましょう。
そのあとで、Sに入った文字列が浮動小数に変換できるかどうかをチェックしましょう。
変換出来たら0 と1の間にあるかをチェックしましょう。
問題なければUに入れれば良いでしょう。

問題があった時に、プログラムを終了するか、再度入力させるかなどは、ご自身の望むようにしてください。

Comments

[izm2192]

回答ありがとうございます
今試して変更してみましたがグラフの外に点が出たりなどしてしないます。。
横軸が83
その中に15,20,34.5,60の位置に点を打ちたい場合どのように書けばいいでしょうか
質問ばかりですいません。

[ppaul]

そういうのは、質問を出すときに書くべきですね。

[izm2192]

すいません。。