请大家帮我看看此程序是否matlab程序，若是，错误在哪里？-经管之家官网！

免费学术公开课,扫码加入

%先定义方差协方差矩阵sigma

A=cholesky(sigma);

Iterations=10000;

theta=zeros(Iteration,5);

y=zeros(5,1);cgpi=zeros(6,1);ex=zeros(7,1);nr=zeros(6,1);spi=zeros(6,1);

rgdp=zeros(6,1);

y(1)=-4.1254;

cgpi(1:2)=[104.4,110.60];

ex(1:3)=[9.2033,10.8357,11.8377];

nr(1:2)=[2.25,2.25];

rgdp(1:2)=[8.9,10.7];rgdp(6)=6;

spi(1:2)=[28.94,30.94];

for i=1:Iterations

for t=1:3

r=normrnd(0,1,6,1);

e=transpose(A)*r;

rgdp(t+2)=3.083+1.3917*rgdp(t+1)-0.6946*rgdp(t)+e(5);

end

for t=1:4

r=normrnd(0,1,6,1)

e=transpose(A)*r; （此处的transpose表示矩阵的转置）

cgpi(t+2)=33.1521+1.2533*cgpi(t+1)-0.5483*cgpi(t)+e(2);

ex(t+3)=1.1657*ex(t+1)-0.7189*ex(t+1)+0.5885*ex(t)+e(3);

nr(t+2)=1.4657*nr(t+1)-0.476*nr(t)+e(4);

spi(t+2)=4.6564+1.5136*spi(t+1)-0.6989*spi(t)+e(6);

y(t+1)=-4.5088+0.046*cgpi(t+2)-0.009*ex(t+3)+0.2429*nr(t+2)-0.1539*rgdp(t+2)-0.0138*spi(t+2) +0.732*y(t)+e(1);

end

theta(i,:)=[cgpi(6),ex(7),nr(6),spi(6),y(5)];

end

p=exp(theta(:,5))/(1+exp(theta(:,5)));

pquantile=sort(p);

m=mean(p)

pq=[pquantile(9000),pquantile(9500),pquantile(9900)]

subplot(2,2,1);

hist(theta(:,5),-6:0.11:0);

subplot(2,2,2);
hist(p,0:0.01:0.3);

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP