matlab难题求解

英文试题，由于自身水平有限，简单的翻译了一下，很多术语不懂，不精确，一共是三问，一问一百个论坛币，如果有版友完成了部分小问，我会另开贴支付论坛币的，说到做到。先放英文试题，点击图片会有大图，接着是中文试题，谢谢大家了

     写一个函数 X = myARsample(c,phi,X0,N)，来定义以下变量X0;X1;…;XN，其符合一阶自相关函数：

Xn= c + φXn-1 +Zn;    n = 1;……;N;

有给定的参数c和时间长度N，随机数Z0，Z1，...，ZN是独立的标准正态分布（即Z = randn（1，N））

1）调用前面的函数编写一个script，用于直观检查的不同φ值的样本路径的属性（|φ|<1，φ=正负1，  |φ|>1）。（提示：在调用该函数的产生循环之前，预先定义向量X，设置X =zero（1,N +1））；

2）编写一个M 为上例的1000倍（即生成M个独立样本路径）和时间总计Nk，使每个样品路径超过给定阈值K。在该情况下，c=0，φ=0.5，X0 =0，K =1，创建Nk的直方图；

3）修改你的函数myARsample，直到当它变为负值时停止，即Xn <0然后重复2)步。

实话实说，这样的题不难。我把函数跟脚本都贴一起了

% script
c = 0;
phi = 0.5;
X0 = 0;
K = 1;
M = 1000;
N = 100;
X = zeros(1,N+1);

for i = 1:M;
    X = [X;myARsample_modify(c,phi,X0,N)];
end

X(1,:) = [];
Counts = sum((X>K),2);
hist(Counts);

%第一问的function
function [result, rnd] = myARsample(c,phi,X0,N)
X = zeros(1,N+1);
X(1) = X0;
Z = randn(1,N);
for i = 2:N+1
    X(i) = c + phi*X(i-1) + Z(i-1);
end
result = X;
rnd = Z;
end

%第三问，修改后的function
function [result, rnd] = myARsample_modify(c,phi,X0,N)
X = zeros(1,N+1);
X(1) = X0;
Z = randn(1,N);
for i = 2:N+1
    X(i) = c + phi*X(i-1) + Z(i-1);
    if X(i)<=0
        break
    end
end
result = X;
rnd = Z;
end

题目没有定义N，我取了100