求助:赋值具有的非单一 rhs 维度多于非单一下标数
错误地方:
V(:,i+1)=(b*z3).^2+b*z3.*sqrt(b^2*(z3.^2)+4*V(:,i)+4*k*a*dt+4*V(:,i)*k*dt+4*k^2*a*dt^2)+2*B*(V(:,i)+k*a*dt)/2*B^2;
请问上面哪里错了
clear all;
close all;
clc;
St=10;
u=0.14;
k=4;
a=0.1;
b=0.7;
c=-0.8;
dt=1/240;
z1=random('normal',0,1,240,1);
z2=random('normal',0,1,240,1);
z3=c*z1+sqrt(1-c^2)*z2;
B=1+k*dt;
for i=1:1:240;
V(1,1)=0.15;
V(:,i+1)=(b*z3).^2+b*z3.*sqrt(b^2*(z3.^2)+4*V(:,i)+4*k*a*dt+4*V(:,i)*k*dt+4*k^2*a*dt^2)+2*B*(V(:,i)+k*a*dt)/2*B^2;
C=c*b.^2*dt^2+sqrt(b.^2*b.^2*dt^2+4*V(:,i)+4*k*a*dt+4*V(:,i)*k*dt+4*k.^2*a*dt^2)*z1;
S(1,1)=10;
S(1,i+1)=S(:,i).*exp((u*dt+dt*V(:,i)/4*B^2)+C/2*B);
end
plot(V) |
|
|