• 回答数

    4

  • 浏览数

    2158

  • 收藏数

    0

作者:貪戀妳的溫柔 发表于 2018-11-20 16:32:26
跳转到指定楼层
求助:赋值具有的非单一 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)

分享:
回复

使用道具

该用户从未签到

新手上路

Rank: 1

积分
41
极客币
109
主题
15
帖子
50
注册时间
2018-10-19
在线时间
2 小时
性别
保密
发表于 2018-11-20 16:32:43 | 显示全部楼层
z1,z2都是240行,1列
在进行V(:,i+1)=...计算时,行列数能对上吗?
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
43
极客币
104
主题
10
帖子
48
注册时间
2018-10-19
在线时间
2 小时
性别
保密
 楼主| 发表于 2018-11-20 16:33:00 | 显示全部楼层
。。。。。我我知道对不上,可是我不知道应该怎么改
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
41
极客币
109
主题
15
帖子
50
注册时间
2018-10-19
在线时间
2 小时
性别
保密
发表于 2018-11-20 16:33:29 | 显示全部楼层
貪戀妳的溫柔 发表于 2018-11-20 16:33
。。。。。我我知道对不上,可是我不知道应该怎么改

如果没有公式,别人很难知道你想干什么,即便改的能执行了,也未必是你想要的。
回复

使用道具 举报

该用户从未签到

新手上路

Rank: 1

积分
43
极客币
104
主题
10
帖子
48
注册时间
2018-10-19
在线时间
2 小时
性别
保密
 楼主| 发表于 2018-11-20 16:33:42 | 显示全部楼层
回复

使用道具 举报

高级模式 评论
您需要登录后才可以回帖 登录 | 立即注册 微信登录