/
alpha_func.m
12 lines (10 loc) · 10.6 KB
/
alpha_func.m
1
2
3
4
5
6
7
8
9
10
11
12
function R=alpha_func(t)
% a=0.4;
% tt=dt:dt:5.5;
% for n=1:length(tt)
% alpha(n)=tt(n)/(a)^2*exp(-tt(n)/(a));
% end
% alpha=alpha/max(alpha);
as=[0.0662791802745652,0.129285482965792,0.189140119526861,0.245960311115695,0.299859411745887,0.350947027788899,0.399329133932628,0.445108185698494,0.488383228616275,0.529250004153169,0.567801052490784,0.604125812241143,0.638310717190200,0.670439290154864,0.700592234037083,0.728847520156204,0.755280473938466,0.779963858040328,0.802967952980068,0.824360635350064,0.844207453680036,0.862571702019593,0.879514491306443,0.895094818584762,0.909369634136376,0.922393906585617,0.934220686037007,0.944901165303202,0.954484739279026,0.963019062515806,0.970550105048695,0.977122206528136,0.982778128705195,0.987559106319041,0.991504896433473,0.994653826268083,0.997042839568284,0.998707541557223,0.999682242511318,1,0.999692659829041,0.998790895725750,0.997324247803194,0.995321159839556,0.992809015407670,0.989814172888817,0.986361999403819,0.982476903693578,0.978182367980237,0.973500978839256,0.968454457111835,0.963063686886233,0.957348743575746,0.951328921120263,0.945022758337587,0.938448064449895,0.931621943810032,0.924560819851571,0.917280458285905,0.909795989568950,0.902121930659393,0.894272206089754,0.886260168370955,0.878098617750442,0.869799821343359,0.861375531655677,0.852837004517645,0.844195016445396,0.835459881448000,0.826641467296776,0.817749211273164,0.808792135410999,0.799778861248532,0.790717624105144,0.781616286897203,0.772482353507138,0.763322981719361,0.754144995736277,0.744954898287235,0.735758882342885,0.726562842447052,0.717372385677869,0.708192842249562,0.699029275765967,0.689886493136493,0.680769054164964,0.671681280821443,0.662627266206845,0.653610883219860,0.644635792935428,0.635705452703702,0.626823123978229,0.617991879881754,0.609214612517845,0.600494040036273,0.591832713459856,0.583233023280230,0.574697205829804,0.566227349436926,0.557825400371075,0.549493168584713,0.541232333258195,0.533044448153968,0.524930946786104,0.516893147411009,0.508932257844998,0.501049380114235,0.493245514942383,0.485521566081148,0.477878344488724,0.470316572361030,0.462836887020442,0.455439844666619,0.448125923993839,0.440895529679169,0.433748995745641,0.426686588804465,0.419708511180215,0.412814903922792,0.406005849709838,0.399281375643186,0.392641455942803,0.386086014541582,0.379614927584244,0.373228025833499,0.366925096986515,0.360705887904675,0.354570106759469,0.348517425097324,0.342547479826059,0.336659875125563,0.330854184285252,0.325129951470739,0.319486693422091,0.313923901085989,0.308441041184003,0.303037557719149,0.297712873422828,0.292466391144153,0.287297495183646,0.282205552573190,0.277189914304094,0.272249916505032,0.267384881571602,0.262594119249162,0.257876927670568,0.253232594350385,0.248660397137074,0.244159605124641,0.239729479525153,0.235369274503517,0.231078237975828,0.226855612372603,0.222700635368128,0.218612540577137,0.214590558219988,0.210633915757466,0.206741838496308,0.202913550166514,0.199148273471456,0.195445230611790,0.191803643784121,0.188222735655346,0.184701729813587,0.181239851196556,0.177836326498217,0.174490384554541,0.171201256709138,0.167968177159535,0.164790383284819,0.161667115955365,0.158597619825332,0.155581143608583,0.152616940338676,0.149704267613539,0.146842387825435,0.144030568376783,0.141268081882402,0.138554206358717,0.135888225400433,0.133269428345211,0.130697110426789,0.128170572917065,0.125689123257546,0.123252075180644,0.120858748821212,0.118508470818739,0.116200574410595,0.113934399516710,0.111709292816043,0.109524607815211,0.107379704909595,0.105273951437281,0.103206721726129,0.101177397134286,0.0991853660844416,0.0972300240921020,0.0953107737881653,0.0934270249360549,0.0915781944436709,0.0897637063704031,0.0879829919294422,0.0862354894856181,0.0845206445489824,0.0828379097643502,0.0811867448970007,0.0795666168147360,0.0779769994664842,0.0764173738576298,0.0748872280222461,0.0733860569923960,0.0719133627646650,0.0704686542640805,0.0690514473055674,0.0676612645530828,0.0662976354765697,0.0649600963068596,0.0636481899886552,0.0623614661317107,0.0610994809603327,0.0598617972613098,0.0586479843303823,0.0574576179173542,0.0562902801699481,0.0551455595764978,0.0540230509075708,0.0529223551566083,0.0518430794796662,0.0507848371343393,0.0497472474179437,0.0487299356050326,0.0477325328843162,0.0467546762950518,0.0457960086629708,0.0448561785358023,0.0439348401184536,0.0430316532079033,0.0421462831278601,0.0412784006632409,0.0404276819945128,0.0395938086319510,0.0387764673498511,0.0379753501207435,0.0371901540496454,0.0364205813083916,0.0356663390700792,0.0349271394436605,0.0342026994087168,0.0334927407504449,0.0327969899948836,0.0321151783444109,0.0314470416135344,0.0307923201650033,0.0301507588462629,0.0295221069262745,0.0289061180327211,0.0283025500896190,0.0277111652553526,0.0271317298611501,0.0265640143500164,0.0260077932161371,0.0254628449447687,0.0249289519526271,0.0244059005287874,0.0238934807761042,0.0233914865531657,0.0228997154167874,0.0224179685650575,0.0219460507809396,0.0214837703764403,0.0210309391373491,0.0205873722685556,0.0201528883399488,0.0197273092329045,0.0193104600873629,0.0189021692495018,0.0185022682200062,0.0181105916029391,0.0177269770552128,0.0173512652366645,0.0169832997607357,0.0166229271457562,0.0162699967668339,0.0159243608083492,0.0155858742170531,0.0152543946557696,0.0149297824577004,0.0146119005813293,0.0143006145659268,0.0139957924876509,0.0136973049162421,0.0134050248723110,0.0131188277852148,0.0128385914515192,0.0125641959940440,0.0122955238214871,0.0120324595886251,0.0117748901570848,0.0115227045566834,0.0112757939473318,0.0110340515814979,0.0107973727672245,0.0105656548316982,0.0103387970853640,0.0101167007865808,0.00989926910681375,0.00968640709635742,0.00947802165058580,0.00927402147672325,0.00907431706113160,0.00887882063710820,0.00868744615318947,0.00850010924195470,0.00831672718932474,0.00813721890435015,0.00796150488948323,0.00778950721132875,0.00762114947186746,0.00745635678014724,0.00729505572443613,0.00713717434483173,0.00698264210632138,0.00683138987228772,0.00668334987845383,0.00653845570726248,0.00639664226268410,0.00625784574544758,0.00612200362868878,0.00598905463401082,0.00585893870795104,0.00573159699884878,0.00560697183410891,0.00548500669785525,0.00536564620896891,0.00524883609950583,0.00513452319348842,0.00502265538606573,0.00491318162303724,0.00480605188073463,0.00470121714625659,0.00459862939805143,0.00449824158684235,0.00440000761689027,0.00430388232758925,0.00420982147538945,0.00411778171604264,0.00402772058716538,0.00393959649111514,0.00385336867817426,0.00376899723003727,0.00368644304359657,0.00360566781502213,0.00352663402413021,0.00344930491903688,0.00337364450109150,0.00329961751008589,0.00322718940973477,0.00315632637342289,0.00308699527021477,0.00301916365112261,0.00295279973562824,0.00288787239845489,0.00282435115658467,0.00276220615651769,0.00270140816176881,0.00264192854059806,0.00258373925397080,0.00252681284374371,0.00247112242107293,0.00241664165504043,0.00236334476149495,0.00231120649210396,0.00226020212361283,0.00221030744730779,0.00216149875867914,0.00211375284728118,0.00206704698678551,0.00202135892522440,0.00197666687542067,0.00193294950560112,0.00189018593019004,0.00184835570077980,0.00180743879727520,0.00176741561920876,0.00172826697722365,0.00168997408472141,0.00165251854967156,0.00161588236658004,0.00158004790861381,0.00154499791987857,0.00151071550784719,0.00147718413593567,0.00144438761622435,0.00141231010232149,0.00138093608236675,0.00135025037217192,0.00132023810849643,0.00129088474245525,0.00126217603305652,0.00123409804086680,0.00120663712180131,0.00117977992103714,0.00115351336704679,0.00112782466575016,0.00110270129478255,0.00107813099787656,0.00105410177935581,0.00103060189873835,0.00100761986544766,0.000985144433629297,0.000963164597071114,0.000941669584225186,0.000920648853329429,0.000900092087627072,0.000879989190682117,0.000860330281788975,0.000841105691474434,0.000822305957090254,0.000803921818494641,0.000785944213820856,0.000768364275331371,0.000751173325355825,0.000734362872311270,0.000717924606803051,0.000701850397804776,0.000686132288915886,0.000670762494695268,0.000655733397069477,0.000641037541814100,0.000626667635106841,0.000612616540150932,0.000598877273867511,0.000585443003655591,0.000572307044218325,0.000559462854454263,0.000546904034412311,0.000534624322309168,0.000522617591607981,0.000510877848157046,0.000499399227387333,0.000488175991567726,0.000477202527116772,0.000466473341969880,0.000455983063000839,0.000445726433496581,0.000435698310684146,0.000425893663308786,0.000416307569262219,0.000406935213260011,0.000397771884567112,0.000388812974770594,0.000380053975598625,0.000371490476784784,0.000363118163976773,0.000354932816688667,0.000346930306295802,0.000339106594071464,0.000331457729264523,0.000323979847217188,0.000316669167522093,0.000309521992217879,0.000302534704022548,0.000295703764603770,0.000289025712885438,0.000282497163389705,0.000276114804613810,0.000269875397440959,0.000263775773584583,0.000257812834065306,0.000251983547719923,0.000246284949741772,0.000240714140251829,0.000235268282899921,0.000229944603495427,0.000224740388666866,0.000219652984549789,0.000214679795502381,0.000209818282848215,0.000205065963645596,0.000200420409482948,0.000195879245299708,0.000191440148232206,0.000187100846484013,0.000182859118220246,0.000178712790485349,0.000174659738143854,0.000170697882843653,0.000166825192001308,0.000163039677808952,0.000159339396262328,0.000155722446209527,0.000152186968419995,0.000148731144673397,0.000145353196867907,0.000142051386147543,0.000138824012048128,0.000135669411661503,0.000132585958817612,0.000129572063284072,0.000126626169982884,0.000123746758223904,0.000120932340954747,0.000118181464026755,0.000115492705476720,0.000112864674824006,0.000110296012382765,0.000107785388588920,0.000105331503341619,0.000102933085358827,0.000100588891546803,9.82977063831192e-05,9.60583413129840e-05,9.38696341585530e-05,9.17304485409775e-05,8.96396733149073e-05,8.75962220151952e-05,8.55990323155375e-05,8.36470654988035e-05,8.17393059388057e-05,7.98747605932667e-05,7.80524585077484e-05,7.62714503303088e-05,7.45308078366622e-05,7.28296234656169e-05,7.11670098645733e-05,6.95420994448705e-05,6.79540439467696e-05,6.64020140138711e-05,6.48851987767636e-05,6.34028054457095e-05,6.19540589121734e-05,6.05382013590070e-05,5.91544918791033e-05,5.78022061023413e-05,5.64806358306445e-05,5.51890886809785e-05,5.39268877361197e-05,5.26933712030273e-05,5.14878920786580e-05,5.03098178230621e-05,4.91585300396076e-05,4.80334241621766e-05,4.69339091491871e-05,4.58594071842929e-05,4.48093533836183e-05,4.37831955093878e-05,4.27803936898121e-05,4.18004201450989e-05,4.08427589194537e-05,3.99069056189431e-05];
R=conv(t,as);
end