Need to write a code for this method....

Hi Everyone,

I need to write a code for the following method.I want to compute 10log s(t1)/s(1)

t1=1,2,3...,N (second)while the data (both s and t)are in one minute.

s        t
0.05     4.7
0.1      4.8   
0.2      4.9
0.1      4.10
0.3      4.11

Could you please guide me?

Thanks,

[ChristianW]

• What class are s and t? Char or Double?
• s=0.05min means 0 min and 0.5 sec?
• Does t allways do 1 sec steps?

[Ara]

No, S is column data which all represent as an integer value. If SI=10logs(t1)/s(1) I want to plot (t,SI)in second. s(t1)denote value s in the first minute. It means s computed in minutes and I want to compute SI while should be calculate in second.So the question is how can I write a code to compute SI. As I assume t1=4.7 is first minutes that S(1)=0.05 for first point then SI=10log(s/0.05)in dB. T just used for plot so I guess the value of s should change and convert to second, right?

Answer 1

[ChristianW]

Since your explanations feel like a puzzle, whats wrong with this?

N = 60; % or N = length(t)
t1 = 1:N; %[sec] timesteps
s = rand(1,N+100); %[any unit] generate random data s

SI = 10*log10(s(t1)/s(1)); %[dB]

plot(t1,SI); xlabel('time [sec]'); ylabel('SI [dB]')

[Ara]

Thanks Christian, this seems good but still some part are not match together which is not reasonable when I compared with *plot (time, s4)*in exact point that this plot increase the SI should have some fluctuate. The code you wrote is not fluctuate in exact increase points. I guess just some small things are not match yet. So I put the data here to make my explanation more clear.

S (in minutes)

0.136027416460800
0.106028827212226
0.0808465761860080
0.122400125714805
0.0790715707634546
0.0954173958143902
0.0719600620622301
0.0832786078654056
0.153006032573883
0.148275109050710
0.0840729558419353
0.0935178422281010
0.0977412421191791
0.113435555021342
0.123715155276142
0.106656997313819
0.0780320042034036
0.0887401268029295
0.0840915504970624
0.110844984473814
0.100027070310991
0.108952113357199
0.121564207890316
0.101530524927236
0.0896254707770063
0.0779490175884212
0.0836013177707146
0.0566822144944955
0.0887442145719934
0.0653785327152575
0.0594292927099086
0.0805781935451521
0.0547014669181732
0.0847473996238233
0.0986969380832050
0.0778304012259991
0.0598654550888908
0.0971791431172348
0.0806545553270242
0.108318490365219
0.0677679655294447
0.0690431141316786
0.0781563987911418
0.0754011812838499
0.0956967248290139
0.0797177641369852
0.0757496636758210
0.0521202396387430
0.0728714463353102
0.0530943210993417
0.0662592120689644
0.0670048162522665
0.0682305471178416
0.0566848991619461
0.0827616866128283
0.0767236920383789
0.0697480340368100
0.114746055134806
0.0845736635425000
0.0700862852846404
0.0585161252989294
0.0618954229729469
0.0609204173409868
0.0780314741626736
0.0785476552418976
0.0645112348354920
0.0608091506271877
0.0469426217738209
0.0726766518216132
0.0635435465251980
0.0563826720189812
0.0774753359011757
0.0624897225469917
0.0525328615915790
0.0763128218846610
0.0644131859559826
0.0730494309423421
0.0691956654350545
0.0496970730727676
0.0566116611662297
0.0549484375301063
0.0558807069837167
0.0552892668426703
0.0520698277796269
0.0715566666076614
0.0450323136758484
0.0488782526590303
0.0636818570787630
0.0576409018319457
0.0727569886333952
0.0601105503135681
0.0536190028907663
0.0603440796101821
0.0574775476164389
0.0655961759251254
0.0418031788264960
0.0585274975374823
0.0879656383765843
0.0474174822612926
0.0768734011931305
0.0648919569746513
0.0858163149698238
0.0919145331762067
0.0918352681980077
0.0832515859308398
0.0633579071545139
0.100287323939768
0.118911287538232
0.147068256374379
0.113292169919196
0.0922938186175001
0.147318654419595
0.0954798334728334
0.106412221854447
0.165544978359357
0.0990780381870776
0.126232370955314
0.160639416990351
0.0907611161070643
0.0765762210676395
0.0948238173825543
0.0957351344909485
0.0693346716441349
0.0618917081118303
0.0551739755772593
0.0632474132356415
0.0623579327832474
0.0653224095774184
0.0980328286034836
0.143193601445735
0.174040530362901
0.175798463113305
0.187343293445482
0.215046327508749
0.218472109970586
0.200057722592756
0.224213539930130
0.242444716747963
0.255535485238352
0.225274175217667
0.207567413781643
0.160459675582372
0.164074152677989
0.119507707851837
0.168734269666834
0.0999956336996771
0.143689466809506
0.201659100218661
0.222302716744533
0.202958083807963
0.283271177984983
0.341578420517163
0.249878128424638
0.275898267876042
0.302890475943038
0.287881537323949
0.215979741994475
0.256893199808792
0.114148091530257
0.142767523628450
0.132827175551541
0.0711602371904423
0.138426956677520
0.106124413642668
0.113959759621544
0.111374245330777
0.103009685214547
0.0601432885698812
0.0739848634046181
0.0869535669193622
0.0881421426106718
0.0951403368451048
0.0562084124931491
0.0971753220267368
0.0634307910324316
0.0489340310724551
0.0368984612145276
0.0386704615953831
0.0821848599378255
0.0514071885537422
0.0544800040473567
0.0762372898062359
0.0794862952715750
0.0426331763770892
0.0693364027117069
0.0761918386639934
0.0645543549499180
0.0683366720655901
0.0751418112304994
0.0509045459659548
0.0703051758620942
0.0601703851076258
0.0719294355879427
0.0837544932824502
0.0602038093811347
0.0604876036886898
0.0601144034653926
0.0889035951972697
0.0856485720896735
0.0557421294444337
0.0539211201385876
0.0622320708638239
0.0668044151460066
0.0735779984710103
0.0652813292373861
0.0601486012472443
0.0626309244143818
0.0501611872666507
0.0671311548016269
0.0680498817412639
0.0544888579069152
0.0571533887971658
0.0713121526459551
0.0688445732719726
0.0678097461136672
0.0702104390457715
0.0877548241408984
0.0961148498204102
0.0896887075166099
0.101688155962236
0.0967310900383119
0.107488544896654
0.166546305858161
0.167274992581079
0.111127209363864
0.111203658788729
0.0934729695206053
0.0944969387440673
0.0686684049618163
0.0904493112411587
0.0944077265693863
0.0664485965540281
0.0938955675471425
0.0968876774827429
0.0723829670571744
0.0806500248915027
0.110307353952490
0.0655479569551943
0.0815076263487043
0.0926289385883267
0.141704098741003
0.149582059833390
0.160607684448783
0.127762019532410
0.0791849504893449
0.0985581026349432
0.130700435408609
0.117573695914520
0.160467373621556
0.102084098698083
0.105751424765816
0.0685919012712142
0.0611823839352472
0.0702803608414186
0.0809974014151071
0.0724013305471661
0.0709120838221526
0.119490952607300
0.0988891897074701
0.110017254028630
0.0868582725190871
0.0939776473848968
0.0799663104250784
0.133606415732928
0.100033621672916
0.0891770402289737
0.100628311736807
0.0972749903315338
0.102328051417976
0.0818176539690549
0.0819921444042050
0.0757512662072391
0.0631519493285837
0.0887366777606644
0.0839079478952977
0.0886051586590758
0.0494228967989534
0.0446795733976948
0.134270086676072
0.0895184838175893
0.0853700007321073
0.122114209406604
0.103844812826641
0.152699760156983
0.108191849156949
0.0978762445284861
0.113832481915313
0.0838864364662131
0.0655057306882383
0.0745176967236642
0.0984767897120941
0.0604829734718789
0.129204322528312
0.113841383507932
0.0884852151774521
0.0734495162339413
0.0636696847172970
0.0749330267639043
0.0781856173154628
0.156127951510932
0.0790311533510678
0.0987533992073184
0.0965363128983079
0.105428925916942
0.0865622391346250
0.0757777861381553
0.0723371011653080
0.0925308165261714
0.101900404194488
0.0796133048993697
0.0867237056692114
0.0701313275434025
0.0816018777173663
0.0770187736529218
0.0810904351942940
0.0787414980362960
0.109144234552266
0.110211274922305
0.104000799809425
0.0920106489652149
0.0581346312020641
0.0899892202655407
0.116318607819214
0.0584792423343531
0.0742817612338318
0.0574369804220243
0.0682469420853418
0.0691338368673402
0.0981439993071405
0.0617780037230081
0.127856288511751
0.0868930882003857
0.125495311031927
0.100472039792173
0.0414138776740358
0.0563675669246066
0.0622408067669435
0.0633024836479581
0.0879882695533899
0.0695928523341299
0.107315504341172
0.105364751449429
0.0474076990055413
0.0497215833215315
0.0891420636736664
0.0609247443983149
0.140201992321793
0.0880514912309837
0.0767405979843785
0.117330311139108
0.0847253506808913
0.0924329887810624
0.103714924986715
0.0699069361365523
0.115809041551167
0.0549508069003541
0.116302445167761
0.116198085182158
0.0563839841089649
0.0773256907890256
0.0687467300749643
0.0531944259579893
0.0592784446489616
0.0952735671631959
0.0645479066043199
0.0682577126484619
0.0471239862490431
0.0818100913090800
0.101733544109109
0.148272743604480
0.0750750455810717
0.0489689331617506
0.0766075662646974
0.0684312932801946
0.142511389983398
0.0855197666039846
0.101694822169076
0.119026051589558
0.0625308037370383
0.0617079177415670
0.0468240803326664
0.0608489322913722
0.104897129608012
0.114168503362355
0.0925551289070465
0.0959274868637764
0.115948698017701
0.0823345886307814
0.0984457071689772
0.0804670868865029
0.0749553157221021
0.0706680656308067
0.0612686778949897
0.0884281543175023
0.0879264696436744
0.0622012636849124
0.0353116698415694
0.130059887190479
0.105221537201278
0.109186997206627
0.0771605705201821
0.0802998656848192
0.128992836874766
0.128045688677909
0.0749726470120937
0.0930697983021345
0.0631507189903646
0.108228828206721
0.0957195688247706
0.0442280895359499
0.0705827107796237
0.0712508591386799
0.0511583740945703
0.0632793589095212
0.0588747565600063
0.0601639990609002
0.0736104353675483
0.0533763934338018
0.0697533145879104
0.0852289243860322
0.0834575087394777
0.107721892227161
0.0714865100910654
0.0721221544672648
0.0602208281327980
0.0795833110143075
0.0681419955680783
0.0779963934550823
0.0865922735583262
0.0662421185651546
0.0668600400837451
0.132167961999117
0.117548100218591
0.137007652107464
0.160180073276922
0.101327554855528
0.0686428955682961
0.121749531107927
0.0806013293004526
0.109645026243784
0.115175291356263
0.0845964083457448
0.100676146742910
0.0640415608804158
0.122564459175570
0.0509603249204712
0.0824972757853252
0.110410185947674
0.0947714419010284
0.0830887769256474
0.0719014508546246
0.0780166938033137
0.0674630690081618
0.0603033730399884
0.110860562694765
0.131726839816341
0.0958620759059598
0.106795482376363
0.0652112722847822
0.0876783073285519
0.106803433147067
0.0832975041642906
0.102706884676734
0.0995309416764455
0.105235571742638
0.0535602407854931
0.0444353570931978
0.0707127443888300
0.0702726124816774
0.0747403537923122
0.0789148462648189
0.0809916721089274
0.108111716737826
0.114471047645245
0.0795278088218203
0.0898542432220093
0.0996636290980817
0.0936421742539119
0.100729867432654
0.0670220672017807
0.101864829077557
0.0870028122304101
0.178635891900256
0.112292025095284
0.131063839269266
0.0996889111987888
0.170593702322214
0.0872404122411168
0.0983317243670627
0.0753249682376302
0.144707526162256
0.154113777888286
0.166141942206055
0.110992170935611
0.0899555304970184
0.108902259356728
0.0771649728762993
0.148728755807342
0.157087250415175
0.124152931177641
0.0699408857893579
0.104134314805447
0.0876954199488206
0.0864672372578192
0.150825538805601
0.239438234630562
0.180134546159253
0.137301025790050
0.107411004203480
0.200867820568652
0.217647084067304
0.111817730454521
0.125158286968942
0.172566181901322
0.113548331203061
0.114626950792560
0.154981768631023
0.160189195762386
0.100825290636824
0.154324335086855
0.209423700425716
0.193318049462537
0.105930395751172
0.146450729171964
0.143580057915436
0.114964346251349
0.186306215033208
0.136376203686714
0.225497726711380
0.271962718025835
0.231092090431499

t (time in hours)
4.78333333333334
4.80000000000000
4.81666666666667
4.83333333333333
4.85000000000000
4.86666666666667
4.88333333333333
4.90000000000000
4.91666666666667
4.93333333333333
4.95000000000000
4.96666666666667
4.98333333333333
5
5.01666666666667
5.03333333333334
5.05000000000000
5.06666666666667
5.08333333333333
5.10000000000000
5.11666666666667
5.13333333333333
5.15000000000000
5.16666666666667
5.18333333333333
5.20000000000000
5.21666666666667
5.23333333333333
5.25000000000000
5.26666666666667
5.28333333333334
5.30000000000000
5.31666666666667
5.33333333333333
5.35000000000000
5.36666666666667
5.38333333333333
5.40000000000000
5.41666666666667
5.43333333333333
5.45000000000000
5.46666666666667
5.48333333333333
5.50000000000000
5.51666666666667
5.53333333333334
5.55000000000000
5.56666666666667
5.58333333333333
5.60000000000000
5.61666666666667
5.63333333333333
5.65000000000000
5.66666666666667
5.68333333333333
5.70000000000000
5.71666666666667
5.73333333333333
5.75000000000000
5.76666666666667
5.78333333333334
5.80000000000000
5.81666666666667
5.83333333333333
5.85000000000000
5.86666666666667
5.88333333333333
5.90000000000000
5.91666666666667
5.93333333333333
5.95000000000000
5.96666666666667
5.98333333333333
6
6.01666666666667
6.03333333333334
6.05000000000000
6.06666666666667
6.08333333333333
6.10000000000000
6.11666666666667
6.13333333333333
6.15000000000000
6.16666666666667
6.18333333333333
6.20000000000000
6.21666666666667
6.23333333333333
6.25000000000000
6.26666666666667
6.28333333333334
6.30000000000000
6.31666666666667
6.33333333333333
6.35000000000000
6.36666666666667
6.38333333333333
6.40000000000000
6.41666666666667
6.43333333333333
6.45000000000000
6.46666666666667
6.48333333333333
6.50000000000000
6.51666666666667
6.53333333333334
6.55000000000000
6.56666666666667
6.58333333333333
6.60000000000000
6.61666666666667
6.63333333333333
6.65000000000000
6.66666666666667
6.68333333333333
6.70000000000000
6.71666666666667
6.73333333333333
6.75000000000000
6.76666666666667
6.78333333333334
6.80000000000000
6.81666666666667
6.83333333333333
6.85000000000000
6.86666666666667
6.88333333333333
6.90000000000000
6.91666666666667
6.93333333333333
6.95000000000000
6.96666666666667
6.98333333333333
7
7.01666666666667
7.03333333333334
7.05000000000000
7.06666666666667
7.08333333333333
7.10000000000000
7.11666666666667
7.13333333333333
7.15000000000000
7.16666666666667
7.18333333333333
7.20000000000000
7.21666666666667
7.23333333333333
7.25000000000000
7.26666666666667
7.28333333333334
7.30000000000000
7.31666666666667
7.33333333333333
7.35000000000000
7.36666666666667
7.38333333333333
7.40000000000000
7.41666666666667
7.43333333333333
7.45000000000000
7.46666666666667
7.48333333333333
7.50000000000000
7.51666666666667
7.53333333333334
7.55000000000000
7.56666666666667
7.58333333333333
7.60000000000000
7.61666666666667
7.63333333333333
7.65000000000000
7.66666666666667
7.68333333333333
7.70000000000000
7.71666666666667
7.73333333333333
7.75000000000000
7.76666666666667
7.78333333333334
7.80000000000000
7.81666666666667
7.83333333333333
7.85000000000000
7.86666666666667
7.88333333333333
7.90000000000000
7.91666666666667
7.93333333333333
7.95000000000000
7.96666666666667
7.98333333333333
8
8.01666666666667
8.03333333333333
8.05000000000000
8.06666666666667
8.08333333333334
8.10000000000000
8.11666666666667
8.13333333333333
8.15000000000000
8.16666666666666
8.18333333333333
8.20000000000000
8.21666666666667
8.23333333333333
8.25000000000000
8.26666666666667
8.28333333333333
8.30000000000000
8.31666666666667
8.33333333333334
8.35000000000000
8.36666666666667
8.38333333333333
8.40000000000000
8.41666666666666
8.43333333333333
8.45000000000000
8.46666666666667
8.48333333333333
8.50000000000000
8.51666666666667
8.53333333333333
8.55000000000000
8.56666666666667
8.58333333333334
8.60000000000000
8.61666666666667
8.63333333333333
8.65000000000000
8.66666666666666
8.68333333333333
8.70000000000000
8.71666666666667
8.73333333333333
8.75000000000000
8.76666666666667
8.78333333333333
8.80000000000000
8.81666666666667
8.83333333333334
8.85000000000000
8.86666666666667
8.88333333333333
8.90000000000000
8.91666666666666
8.93333333333333
8.95000000000000
8.96666666666667
8.98333333333333
9
9.01666666666667
9.03333333333333
9.05000000000000
9.06666666666667
9.08333333333334
9.10000000000000
9.11666666666667
9.13333333333333
9.15000000000000
9.16666666666666
9.18333333333333
9.20000000000000
9.21666666666667
9.23333333333333
9.25000000000000
9.26666666666667
9.28333333333333
9.30000000000000
9.31666666666667
9.33333333333334
9.35000000000000
9.36666666666667
9.38333333333333
9.40000000000000
9.41666666666666
9.43333333333333
9.45000000000000
9.46666666666667
9.48333333333333
9.50000000000000
9.51666666666667
9.53333333333333
9.55000000000000
9.56666666666667
9.58333333333334
9.60000000000000
9.61666666666667
9.63333333333333
9.65000000000000
9.66666666666666
9.68333333333333
9.70000000000000
9.71666666666667
9.73333333333333
9.75000000000000
9.76666666666667
9.78333333333333
9.80000000000000
9.81666666666667
9.83333333333334
9.85000000000000
9.86666666666667
9.88333333333333
9.90000000000000
9.91666666666666
9.93333333333333
9.95000000000000
9.96666666666667
9.98333333333333
10
10.0166666666667
10.0333333333333
10.0500000000000
10.0666666666667
10.0833333333333
10.1000000000000
10.1166666666667
10.1333333333333
10.1500000000000
10.1666666666667
10.1833333333333
10.2000000000000
10.2166666666667
10.2333333333333
10.2500000000000
10.2666666666667
10.2833333333333
10.3000000000000
10.3166666666667
10.3333333333333
10.3500000000000
10.3666666666667
10.3833333333333
10.4000000000000
10.4166666666667
10.4333333333333
10.4500000000000
10.4666666666667
10.4833333333333
10.5000000000000
10.5166666666667
10.5333333333333
10.5500000000000
10.5666666666667
10.5833333333333
10.6000000000000
10.6166666666667
10.6333333333333
10.6500000000000
10.6666666666667
10.6833333333333
10.7000000000000
10.7166666666667
10.7333333333333
10.7500000000000
10.7666666666667
10.7833333333333
10.8000000000000
10.8166666666667
10.8333333333333
10.8500000000000
10.8666666666667
10.8833333333333
10.9000000000000
10.9166666666667
10.9333333333333
10.9500000000000
10.9666666666667
10.9833333333333
11
11.0166666666667
11.0333333333333
11.0500000000000
11.0666666666667
11.0833333333333
11.1000000000000
11.1166666666667
11.1333333333333
11.1500000000000
11.1666666666667
11.1833333333333
11.2000000000000
11.2166666666667
11.2333333333333
11.2500000000000
11.2666666666667
11.2833333333333
11.3000000000000
11.3166666666667
11.3333333333333
11.3500000000000
11.3666666666667
11.3833333333333
11.4000000000000
11.4166666666667
11.4333333333333
11.4500000000000
11.4666666666667
11.4833333333333
11.5000000000000
11.5166666666667
11.5333333333333
11.5500000000000
11.5666666666667
11.5833333333333
11.6000000000000
11.6166666666667
11.6333333333333
11.6500000000000
11.6666666666667
11.6833333333333
11.7000000000000
11.7166666666667
11.7333333333333
11.7500000000000
11.7666666666667
11.7833333333333
11.8000000000000
11.8166666666667
11.8333333333333
11.8500000000000
11.8666666666667
11.8833333333333
11.9000000000000
11.9166666666667
11.9333333333333
11.9500000000000
11.9666666666667
11.9833333333333
12
12.0166666666667
12.0333333333333
12.0500000000000
12.0666666666667
12.0833333333333
12.1000000000000
12.1166666666667
12.1333333333333
12.1500000000000
12.1666666666667
12.1833333333333
12.2000000000000
12.2166666666667
12.2333333333333
12.2500000000000
12.2666666666667
12.2833333333333
12.3000000000000
12.3166666666667
12.3333333333333
12.3500000000000
12.3666666666667
12.3833333333333
12.4000000000000
12.4166666666667
12.4333333333333
12.4500000000000
12.4666666666667
12.4833333333333
12.5000000000000
12.5166666666667
12.5333333333333
12.5500000000000
12.5666666666667
12.5833333333333
12.6000000000000
12.6166666666667
12.6333333333333
12.6500000000000
12.6666666666667
12.6833333333333
12.7000000000000
12.7166666666667
12.7333333333333
12.7500000000000
12.7666666666667
12.7833333333333
12.8000000000000
12.8166666666667
12.8333333333333
12.8500000000000
12.8666666666667
12.8833333333333
12.9000000000000
12.9166666666667
12.9333333333333
12.9500000000000
12.9666666666667
12.9833333333333
13
13.0166666666667
13.0333333333333
13.0500000000000
13.0666666666667
13.0833333333333
13.1000000000000
13.1166666666667
13.1333333333333
13.1500000000000
13.1666666666667
13.1833333333333
13.2000000000000
13.2166666666667
13.2333333333333
13.2500000000000
13.2666666666667
13.2833333333333
13.3000000000000
13.3166666666667
13.3333333333333
13.3500000000000
13.3666666666667
13.3833333333333
13.4000000000000
13.4166666666667
13.4333333333333
13.4500000000000
13.4666666666667
13.4833333333333
13.5000000000000
13.5166666666667
13.5333333333333
13.5500000000000
13.5666666666667
13.5833333333333
13.6000000000000
13.6166666666667

[Ara]

X Coordinate axes should be time and Y coordinate axes should be in dB. S data are in minutes and time I guess no change but SI should be compute in second. Is it possible? I think s should be convert to second.

[ChristianW]

The unit of S is irrelevant. Unit of SI is dB. With the conversion to dB, SI describes a relativ change to your initial value (S(1)).

What is the unit of T? First you said minutes, but in the data its hours.

Try this

SI = 10*log10(S/S(1)); %[dB] !NOTE that T and unit[S] are irrelevant for SI

Tmin = T-T(1);         %[min] Tmin starts at 0 min
Tsec = round(Tmin*60); %[sec] MinToSec (round for numerical errors)

figure('units','norm','outerpos',[0 0 1 1])
subplot(411);
plot(T,S); 
xlabel('time - T [min]'); 
ylabel('S'); axis tight

subplot(412);
plot(Tmin,S); 
xlabel('time - Tmin - T with start at 0 [min]'); 
ylabel('S'); axis tight

subplot(413);
plot(Tsec,S); 
xlabel('time - Tsec - T with start at 0 in seconds [sec]'); 
ylabel('S'); axis tight

subplot(414);
plot(Tsec,SI);  
xlabel('time [sec]'); 
ylabel('SI [dB]'); axis tight

[Ara]

Thanks for the code. But the previous code was closer to the plot that I wish. I do not understand why you use SI = 10*log*10(S/S(1)); * instead of 10logS(t)/S(1) . So sorry to make you confuse. Actually, receiver provide the value of data (s) in one minute. It means in one minute we have one (s) data and times comes out exactly as I put later. But the problem is sampling rate of SI is 1 Hz means one second so for calculating SI I should have able to compute by using s data that are in minute. I do not know how to match these.

(t=1, 2, 3, ..., N)

If S(1), S(2), S(3), ..., S(N), then can plot SI in dB based on s column S(1)is first s data in one minute, S(2) is second minutes and so on.

I really appreciate your help.

[ChristianW]

I use log10(), the decadic logarithm with base 10. That's how the unit Bel is defined. You use log(), the natural logarithm with base e.

Regarding the 1 Hz sample rate of SI, you should check the bottom plot from my last code. It shows SI with 1Hz sample rate over the full range of T (530 sec).

[Ara]

Thanks for your description. Could you please run this code? I put the data here: http://www.sendspace.com/file/o35lrv The figure should be look like this, http://www.sendspace.com/file/78hxof

clear;clc;close all;
data=xlsread('1608_1.xls');
data_nozero=find(data(:,3)==2);
newdata=data(data_nozero,:);
data_filterr=find(newdata(:,25)>60);
data_filtered=newdata(data_filterr,:);
elev_cutof20=find(data_filtered(:,6)>=15);
data_cutoff15=data_filtered(elev_cutof20,:);
r=data_cutoff15(:,2);
time=(r./3600)-24*1;
s4r=data_cutoff15(:,8);
s4cor=data_cutoff15(:,9);
s4=sqrt(s4r.^2-s4cor.^2);
PRN=data_cutoff15(:,3);
Tmin = time-time(1);         %[min] Tmin starts at 0 min
Tsec = round(Tmin*60); 
SI = 10*log10(s4/s4(1));
subplot(211);
plot (time+8, s4);
axis([6 18 0 0.4])
xlabel('Time(LT)');
ylabel('S4 Index');
title('1 November 2010');
subplot(212);plot(time,SI);
xlabel('Time(LT)');
ylabel('SI (dB)');

[ChristianW]

Your definition for SI=f(S) doesn't fit to the picture.

[Ara]

Could you please help me how to plot like that picture?

As they explained me, If S(1), S(2), S(3), ..., S(N) denote signal intensity observations during s4 period, you can find a value of signal intensity around which signal intensity fluctuates slightly when s4 do not occur. This value is defined as Sa. Then compute 10logS(t)/Sa (t=1, 2, 3, ..., N), you can plot SI. But I can't plot SI

[ChristianW]

Is S(2) == s4(2)?

Where do you calculate Sa? Why is Sa = S(1)?

[Ara]

Exactly, s(2)==s4(2) and s(a)=s(1) because the figure supposed to show the signal that look like s4 so the point where s4 are not increase the SI should show the value of around 0 dB with slightly fluctuations and where s4 is increase the SI should show large fluctuations. Another question is, How to convert S4 to second first? and then calculate it by using that relation.

[Ara]

Thanks, Christian. Your code works properly for the other types of my data. Sorry for the late response.