x=shenshi;
y=hushi;
a=mean(x);
b=mean(y)
s1=std(x)
s2=std(y)
day_num=20;
sData=zeros(1,day_num);
sData(1)=5503.93;
rad=normrnd(0,1,1,day_num);
for ii=1:19;
sData(ii+1)=sData(ii)*exp(-0.00084828+0.0003497*rad(ii))
end
(2)
clc
clear
num=200;
hushiprice1=hushiprice(1:200);
L_hushiprice1=length(hushiprice1);
Tcut=24;
sData=zeros(1,Tcut+1);
sData(1)=hushiprice1(end);
Mnum=10000;
for ii=1:num-1;
VarS(ii)=(hushiprice1(ii)-hushiprice1(end))^2;
end
VarSS=(sum(VarS)/(num-1))^0.5;
clear VarS
for kk=1:Mnum
rad=normrnd(0,1,1,Tcut);
for jj=1:Tcut
sData(jj+1)=sData(jj)+rad(jj)*VarSS/Tcut;
end
ST(kk)=sData(end);
end
ST_arrange=sort(ST);
ST_50=ST_arrange(Mnum*0.95+1);
VaR=sData(1)-ST_50;