请高手帮我看下这个简单的迭代计算程序 - 经管之家

› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › MATLAB等数学软件专版

请高手帮我看下这个简单的迭代计算程序

5125

27

收藏 2012-03-05

zs=zeros(23,1);
hs(1)=3;
>> for t=1:23;
hs(t+1)=-0.86+0.9*hs(t)+0.4*abs(zs(t)/sqrt(hs(t)))+0.03*zs(t)/sqrt(hs(t));
end;

就用这个计算数组hs的值，请帮我看下我的运算法则对不对，因为计算出来的hs不是一个24行1列的数组，而是转置过来的，我就觉得我程序有些运算符号是不是有问题！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

2012-3-5 20:36:29

其中我表达zs(t)/sqrt(hs(t))，是指zs与hs的绝对值的比！我觉得这里因为都是数组，好像这样表达是不对的，但是我是初学者，所以我不太确定，求高手指点啊，谢谢了！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-5 20:39:19

不好意思我说错了，，，其中我表达zs(t)/sqrt(hs(t))，是指zs与hs的平方根的比！我觉得这里因为都是数组，好像这样表达是不对的，但是我是初学者，所以我不太确定，求高手指点啊，谢谢了！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-5 20:40:37

zs(t)./sqrt(hs(t))

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-5 20:43:36

liuxin9023 发表于 2012-3-5 20:40
zs(t)./sqrt(hs(t))

你好，我这样修改之后试了下，hs还是一个1行24列的数组，但是我觉得应该得到24行1列的数组才对啊，而且运算的结果也非常不对。

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-5 20:45:26

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

2012-3-6 15:11:12

呜呜！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-6 15:28:36

跪求

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-6 19:51:13

我发现我的zs,hs这些都是数列啊，，，但是我搞不懂数列运算要遵循什么原则，，，跪求有心人帮我看看呀！~谢谢了，，万分感谢！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-6 20:17:25

zs = zeros(23,1);
hs = zeros(24,1);
hs(1) = 3;
for t = 1:23;
hs(t+1) = -0.86+0.9*hs(t)+0.4*abs(zs(t)/sqrt(hs(t)))+0.03*zs(t)/sqrt(hs(t));
end

祝好运！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-6 21:02:42

麻烦你转置下行吗

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-7 12:26:05

liuxin9023 发表于 2012-3-6 21:02
麻烦你转置下行吗

呵呵不是转置的问题

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-7 13:30:44

我的看法和10楼的一样！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-9 14:08:22

我已经把程序改好了，贴在十楼，楼主没看到吗？

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-9 16:38:30

lustboy 发表于 2012-3-6 20:17
zs = zeros(23,1);
hs = zeros(24,1);
hs(1) = 3;

您好，我看到了，这里您给我加了一个限制hs是一个列向量，谢谢，我最近遇到新的问题，加入log运算之后就又不行了，，哎！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-9 17:11:24

怎不行，你具体贴出来

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 12:21:36

lustboy 发表于 2012-3-9 17:11
怎不行，你具体贴出来

zf=zeros(23,1);
hf(1)=0.0001259318;
for t=1:23;
log(hf(t+1))=0.000004+0.177428*zf(t)/sqrt(hf(t))+0.808882*log(hf(t));
end

这个程序运行之后显示错误如下：
Subscript indices must either be real positive integers or logicals

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 12:23:25

lustboy 发表于 2012-3-9 17:11
怎不行，你具体贴出来

加入log之后，就显示这样的错误，我不清楚是为什么，谢谢您帮我看看！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 16:26:16

程序里的这个等号是赋值的意思，也就是说把等号右边的结果赋给左边的变量。
所以，等号左边不能出现函数操作，否则就变成了逻辑判断，而hf（t+1)还没有被赋值，计算机认为它是0.
修改后的程序为：
zf=zeros(23,1);
hf(1)=0.0001259318;
for t=1:23;
hf(t+1)=exp（0.000004+0.177428*zf(t)/sqrt(hf(t))+0.808882*log(hf(t))）;
end
楼主看来还不太理解计算机语言的基本逻辑，建议看看编程的入门教材。

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 17:06:31

lustboy 发表于 2012-3-10 16:26
程序里的这个等号是赋值的意思，也就是说把等号右边的结果赋给左边的变量。
所以，等号左边不能出现函数操 ...

明白了，谢谢你啊！~其实我就是没认真学习matlab，可惜最近论文里面有些小的运算需要用到它，又比较急，就随便看了几页书然后去编简单的程序，都怪平时不认真啊，好的谢谢你的提点。感激涕零！！！！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 17:16:46

lustboy 发表于 2012-3-10 16:26
程序里的这个等号是赋值的意思，也就是说把等号右边的结果赋给左边的变量。
所以，等号左边不能出现函数操 ...

今天自己又编写了一个小的运算程序，是这样的：
a(2)=0.2174;
for t=1:23;
x(t+2)=2.3823*br(t+1)-2.1393*(-0.000139+0.021936*br(t)-0.5)*(0.000253-0.009827*br(t)-0.5);
a(t+2)=0.3958*a(t+1)+exp(x(t+2));
end
运算结果没有什么问题，就是需要做点修改，请教下您，就是我这个a(2)实际上是个行向量，有21个分量，这里的0.2174实际上是第一个分量，也就是其实我有21个初始值，我自己计算的时候是把每个初始值都带一遍，这样特别麻烦，但是下面是个循环程序，初始值为向量的赋值方法我不会，请问您能懂我的意思么？能帮我解决一下么？万分感谢！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 18:34:30

说实话，没看懂你啥意思，你最好把原始问题也一起说清楚。
现在的程序里，br没有赋初值，x都是零。a不会变。

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 19:02:16

烈女乖乖发表于 2012-3-10 17:16
今天自己又编写了一个小的运算程序，是这样的：
a(2)=0.2174;
for t=1:23;

你好，我的br赋值了，我忘记贴上来了，是这样的：
br=zeros(24,1);
a(2)=0.2174;
for t=1:23;
x(t+2)=2.3823*br(t+1)-2.1393*(-0.000139+0.021936*br(t)-0.5)*(0.000253-0.009827*br(t)-0.5);
a(t+2)=0.3958*a(t+1)+exp(x(t+2));
end

我的意思就是，我的这个程序是为了求得a，最终a有23个值（从a(3)到a(25)均可以得到）。但是实际情况是，我每次运算这个程序，都要带入不同的a(2)的值，这里我带入的是0.2174，我一共有21个需要带入，我的做法就是每次运行程序就修改一下a(2)的值，这样做21遍，我感觉太麻烦，每一个a(2)对应求出23个a的值，那么如果能够一次带入a(2)的21个值，每个值对应得到的a，最终能够直接得到a类似于一个23行21列的矩阵，即每一个a(2)对应a的一列！
不好意思，我的表达能力也不太好，不知道您能不能看明白，嘿嘿！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 19:02:51

lustboy 发表于 2012-3-10 18:34
说实话，没看懂你啥意思，你最好把原始问题也一起说清楚。
现在的程序里，br没有赋初值，x都是零。a不会变 ...

你好，我的br赋值了，我忘记贴上来了，是这样的：
br=zeros(24,1);
a(2)=0.2174;
for t=1:23;
x(t+2)=2.3823*br(t+1)-2.1393*(-0.000139+0.021936*br(t)-0.5)*(0.000253-0.009827*br(t)-0.5);
a(t+2)=0.3958*a(t+1)+exp(x(t+2));
end

我的意思就是，我的这个程序是为了求得a，最终a有23个值（从a(3)到a(25)均可以得到）。但是实际情况是，我每次运算这个程序，都要带入不同的a(2)的值，这里我带入的是0.2174，我一共有21个需要带入，我的做法就是每次运行程序就修改一下a(2)的值，这样做21遍，我感觉太麻烦，每一个a(2)对应求出23个a的值，那么如果能够一次带入a(2)的21个值，每个值对应得到的a，最终能够直接得到a类似于一个23行21列的矩阵，即每一个a(2)对应a的一列！
不好意思，我的表达能力也不太好，不知道您能不能看明白，嘿嘿！~

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-10 23:59:25

还是没搞懂你的意思，a要是矩阵的话，需要用两层嵌套的循环控制。
另外，你的br都是零，x也都是零，因而a(t+2)=0.3958*a(t+1)+exp（0）。
似乎你想递归br，a还有x，但你的程序里br没有递归。
你最好找个明白人当面交流吧。

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-12 15:48:04

lustboy 发表于 2012-3-10 23:59
还是没搞懂你的意思，a要是矩阵的话，需要用两层嵌套的循环控制。
另外，你的br都是零，x也都是零，因而a( ...

好的，，，还是非常感谢你啊，，耽误你这么久时间为我解答问题！！！谢谢！！！

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-12 19:04:30

客气了，我能力不济，没帮到你什么，你找个人面谈，可能收获会比较直接。

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

2012-3-12 19:18:55

zs=zeros(23,1);
hs=zeros(24,1);
hs(1)=3;
for t=1:23;
zh = zs(t)/sqrt(hs(t));
hs(t+1)=-0.86+0.9*hs(t)+0.4*zh + 0.03*zh;
end;
提醒你，zs(t) =0, zh =0, 後面兩項可以合併，而且可以去掉

二维码

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

相关推荐

栏目导航

热门文章

推荐文章

扫码加好友，拉您进群

各岗位、行业、专业交流群