Ошибка моделирования Matlab

Question

Я совершенно не знаком с Matlab. Я пытаюсь имитировать комбинированный процесс Wiener и Poisson.

Почему у меня возникает несоответствие размеров подзаголовка подзаголовка?

Я пытаюсь моделировать

Z(t)=lambda*W^2(t)-N(t) 

Где W это процесс Винера и N является пуассоновским процессом.

код я использую ниже:

T=500 
dt=1 
K=T/dt 
W(1)=0 
lambda=3 
t=0:dt:T 
for k=1:K 
r=randn 
W(k+1)=W(k)+sqrt(dt)*r 
N=poissrnd(lambda*dt,1,k) 
Z(k)=lambda*W.^2-N 
end 
plot(t,Z) 

Answer 1

Это правда, что некоторые индексации не хватает, но я подумайте, что вам удастся переписать ваш код более «способом Matlab». Следующий код использует тот факт, что базовые переменные Matlab являются матрицами и вычисляют результаты в векторном виде. Постарайтесь понять этот вид письма, так как это способ более эффективно использовать Matlab, наряду с написанием короче и читаемый код:

T = 500; 
dt = 1; 
K = T/dt; 
lambda = 3; 
t = 1:dt:T; 
sqdtr = sqrt(dt)*randn(K-1,1); % define sqrt(dt)*r as a vector 
N = poissrnd(lambda*dt,K,1); % define N as a vector 
W = cumsum([0; sqdtr],1); % cumulative sum instead of the loop 
Z = lambda*W.^2-N; % summing the processes element-wiesly 
plot(t,Z) 

Пример для результата:

Answer 2

вы забыли индекс

Z(k)=lambda*W.^2-N 

должно быть

Z(k)=lambda*W(k).^2-N(k)