| 查看: 755 | 回复: 0 | |||
[交流]
【求助】高斯分布下，求半宽度
|
|
最近在做高斯分布的振幅，半高宽NL=30，我只是验证一下我的程序能否算出我输入的高斯分布的振幅的半高宽，可是做出的程序算出来的半高宽是50，请教高手，我错在哪？ clc,clear,close all; t=200;%时间 m=600;%时间份数 stept=t/m; Y2=zeros(m,1); Ip0=zeros(m,1); NL=30; for i=1:m%for t t2=i*stept-t/2; Ep0(i,1)=exp(-(t2*t2)/(NL*NL)); Y2(i,1)=i; end%for i plot(Y2(:,1),Ep0(:,1),'b-'); %计算 temp=max(Ep0(:,1)); Ztemp=-temp+Ep0(:,1);%Φ-Φm center=Center(m,Ztemp);%center是编号 Ztemp1=Ep0(:,1)-(temp/2);%Φ-Φm/2 Ztemp2=abs(Ztemp1); Ztemp3=submatrix(1,center,Ztemp2); Ztemp4=min(Ztemp3); Ztemp5=Ztemp3-Ztemp4; T1=Center(center,Ztemp5); Tleft=T1*stept; Ztemp12=submatrix(center+1,m,Ztemp2); Ztemp13=min(Ztemp12); Ztemp15=Ztemp12-Ztemp13; kkk=m-center; T2=Center(kkk,Ztemp15); Tright=(T2+center)*stept; ma=Tright-Tleft |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复