編集履歴

質問編集履歴

補足説明の追加、プログラムの修正

2020/10/21 05:03

投稿

Yu_a

スコア3

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -145,3 +145,421 @@
 初めての投稿でおかしな部分があるとは思いますがご容赦ください。
 https://www.jstage.jst.go.jp/article/jvsj2/56/6/56_13-LC-019/_pdf
+[Si2p data]
+x	y
+106.95	100
+106.9	200
+106.85	329
+106.8	315
+106.75	323
+106.7	309
+106.65	333
+106.6	307
+106.55	315
+106.5	311
+106.45	290
+106.4	317
+106.35	346
+106.3	317
+106.25	293
+106.2	287
+106.15	295
+106.1	306
+106.05	291
+106	306
+105.95	312
+105.9	302
+105.85	272
+105.8	309
+105.75	284
+105.7	270
+105.65	286
+105.6	266
+105.55	281
+105.5	256
+105.45	272
+105.4	288
+105.35	295
+105.3	299
+105.25	292
+105.2	311
+105.15	307
+105.1	305
+105.05	291
+105	324
+104.95	333
+104.9	350
+104.85	383
+104.8	424
+104.75	430
+104.7	469
+104.65	526
+104.6	548
+104.55	608
+104.5	638
+104.45	700
+104.4	793
+104.35	887
+104.3	983
+104.25	1041
+104.2	1146
+104.15	1222
+104.1	1312
+104.05	1389
+104	1558
+103.95	1602
+103.9	1700
+103.85	1840
+103.8	1983
+103.75	2073
+103.7	2147
+103.65	2257
+103.6	2347
+103.55	2370
+103.5	2432
+103.45	2484
+103.4	2487
+103.35	2473
+103.3	2508
+103.25	2497
+103.2	2508
+103.15	2454
+103.1	2327
+103.05	2289
+103	2326
+102.95	2241
+102.9	2127
+102.85	1871
+102.8	1855
+102.75	1685
+102.7	1602
+102.65	1453
+102.6	1411
+102.55	1300
+102.5	1154
+102.45	1067
+102.4	987
+102.35	934
+102.3	842
+102.25	790
+102.2	742
+102.15	683
+102.1	644
+102.05	629
+102	573
+101.95	569
+101.9	579
+101.85	567
+101.8	536
+101.75	524
+101.7	521
+101.65	508
+101.6	517
+101.55	514
+101.5	478
+101.45	440
+101.4	440
+101.35	485
+101.3	467
+101.25	460
+101.2	468
+101.15	483
+101.1	485
+101.05	509
+101	494
+100.95	480
+100.9	490
+100.85	487
+100.8	485
+100.75	473
+100.7	476
+100.65	456
+100.6	436
+100.55	462
+100.5	506
+100.45	559
+100.4	661
+100.35	729
+100.3	872
+100.25	998
+100.2	1199
+100.15	1579
+100.1	1870
+100.05	2273
+100	2656
+99.95	3019
+99.9	3227
+99.85	3339
+99.8	3207
+99.75	3017
+99.7	2852
+99.65	2683
+99.6	2671
+99.55	2893
+99.5	3389
+99.45	3916
+99.4	4621
+99.35	5264
+99.3	5484
+99.25	5662
+99.2	5299
+99.15	4596
+99.1	3847
+99.05	2929
+99	2234
+98.95	1572
+98.9	1045
+98.85	670
+98.8	485
+98.75	357
+98.7	299
+98.65	226
+98.6	170
+98.55	174
+98.5	187
+98.45	149
+98.4	117
+98.35	131
+98.3	142
+98.25	147
+98.2	154
+98.15	130
+98.1	120
+98.05	125
+98	116
+97.95	128
+97.9	130
+97.85	122
+97.8	130
+97.75	110
+97.7	121
+97.65	125
+97.6	125
+97.55	135
+97.5	120
+97.45	127
+97.4	127
+97.35	120
+97.3	114
+97.25	115
+97.2	100
+97.15	106
+97.1	122
+97.05	117
+97	119
+96.95	115
+96.9	135
+96.85	113
+![イメージ説明](bfe6c03380ecd6dfc16ac08bc84a257b.png)

指摘いただいた点を修正して再びプログラムを作りました。自分が実現したい事をまとめてPDFにしたので、そちらを参照していただけるとありがたいです。Si2pのデータについても同様に上げておきます。

2020/10/21 05:03

投稿

Yu_a

スコア3

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -26,9 +26,9 @@
+```Python
-```Python
 import numpy as np
@@ -50,49 +50,51 @@
-class function:
-    def calculation(self):
-        k=0
+k = 1
-        while True:
+Flist =[]
-             self.R = h * (y1 - 0.5*(y[1] - y[-1]) + (y[self.x] - y[-1]))
+while True:
-             self.Q = h * (y1 - 0.5*(y[self.x] - y[-1] + y[202] - y[-1]))
+    for i in range(203):
-             self.B = (y[1] - y[202])*(self.Q/self.R) + y[-1]
+        R1 = h * (y1 - 0.5*(y[1] - y[-1]) + (y[k] - y[-1]))
-             self.F = y[self.x] - self.B
+        Q1 = h * (y1 - 0.5*(y[k] - y[-1] + y[203] - y[-1]))
-             self.Blist.append(self.B)
+        B1 = (y[1] - y[203])*(Q1/R1) + y[-1]
-             while True:
+        F = y[k] - B1
-                 self.k = self.B
+        k += 1
-                 self.Q = h * (y1 - 0.5*(y[self.F] - y[-1] + y[202] - y[-1]))
+        Blist.append(B1)
-                 self.R = h * (y1 - 0.5*(y[1] - y[-1]) + (y[self.F] - y[-1]))
+        Flist.append(F)
-                 self.B = (y[1] - y[202])*(self.Q/self.R) + y[-1]
+    for i in range(100):
-                 self.Blist.append(self.B)
+        t = 1
-                 if round(self.k,5) == round(self.B,5):
+        z = Flist[t]
-                     break
+        Q = h * (y1 - 0.5*(z - y[-1] + y[203] - y[-1]))
-             self.x += 1
+        R = h * (y1 - 0.5*(y[1] - y[-1]) + (z - y[-1]))
+        B = (y[1] - y[203])*(Q/R) + y[-1]
+        Blist.append(B)
-             yield self.F
+        Flist.append(F)
+        t += 1
-classInstance = function()
-print(function())
+print(Blist)
@@ -117,8 +119,6 @@
 plt.ylabel("Intensity[a.u.]",fontname="MS Gothic")
 ax.invert_xaxis()
-#plt.legend(["name"], loc=2, title="LABEL NAME",fontsize = 13,prop = {"size":15} )

ご指摘いただいた通り変更してみたつもりです。ご確認よろしくお願いします。

2020/10/21 04:28

投稿

Yu_a

スコア3

test CHANGED Viewed

File without changes

test CHANGED Viewed

@@ -17,6 +17,18 @@
 ### 該当のソースコード
+```Python
 import numpy as np
@@ -94,6 +106,8 @@
 x_min, x_max = ax.get_xlim()
+#ax.set_xlim(97,102)
 plt.title("Shirley",fontname="MS Gothic")
 plt.plot(x,y,color="black")
@@ -104,17 +118,13 @@
 ax.invert_xaxis()
+#plt.legend(["name"], loc=2, title="LABEL NAME",fontsize = 13,prop = {"size":15} )
 plt.show()
-```ここに言語名を入力
-Python3.7
 ```