查看: 692|回复: 6|关注: 0

[已解决] matlab 向量中的值加起来的和等于一个值后输出序号

[复制链接]

新手

11 麦片

财富积分


050


7

主题

31

帖子

0

最佳答案
%对1000个U进行数列求和,从第一个U值开始加,USUM=U1+U2+...U(N),当USUM=569.0349,输出N,我只遍了前面产生U的程序,后面的怎么弄
%都是标量计算
clear all
clc
close
PI=3.1415926
PPI=10;
POR=0.9;
AL=0.01;
A1PL=(1.0-0.5)*AL;
DI2=0.0254/PPI;
DI1=DI2*1.18*sqrt((1.0-POR)/(3.0*PI))*(1.0/(1.0-exp(-(1.0-POR)/0.04)));
K=0.00073*(1-POR)^(-0.224)*(DI1/DI2)^(-1.11)*DI2^2;
DA=K/AL^2.0;
H=A1PL/AL;
S=(POR/DA)^0.5;


    RO=1.205;
    CF=1005.0;
    VI=0.00001506;
    UF=0.0000181;
    RETOTAL=1511.381;
UM=RETOTAL*VI/(2.0*AL);

C3=(S*DA*exp(S-S*H)-S*DA+S*H^2/2.0+H*exp(S-S*H))/(H*(exp(S*H)+exp(2*S-S*H))-S*DA*(exp(S*H)-exp(2*S-S*H)));
C4=exp(S)-C3*exp(2*S);
C1=(C3*exp(S*H)-C4*exp(-S*H)-S*H)/(S*DA);


P=1.0/(1.0/(6.0*DA)*H^3+1.0/2.0*C1*H^2+S^2*DA*(C3*(exp(S)-exp(S*H))/S-C4*(exp(-S)-exp(-S*H))/S-1.0+H));
M=1000;
for X=(1:M)/1000
    if X>H
        U(X*1000)=P*(C3*exp(S*X)+C4*exp(0-S*X)-1.0)
    else
        U(X*1000)=P*(1.0/(2.0*DA)*X^2+C1*X)
    end

end

X = (1:M)/1000;

plot(X,U,'r')
grid on

论坛优秀回答者

权威

8043 麦片

财富积分



5

主题

9478

帖子

1828

最佳答案
  • 关注者: 417
发表于 2018-6-11 07:35:08 | 显示全部楼层 |此回复为最佳答案
  1. CU=cumsum(U)
  2. find(abs(CU-569.0349)<=1e-4)
复制代码

但事实上,没有累加到569.0349这个数,相近的数是 569.5368

新手

11 麦片

财富积分


050


7

主题

31

帖子

0

最佳答案
 楼主| 发表于 2018-6-11 07:38:11 | 显示全部楼层
本帖最后由 邢打铁 于 2018-6-11 07:39 编辑
jingzhaos 发表于 2018-6-11 07:35
但事实上,没有累加到569.0349这个数,相近的数是 569.5368

兄弟,非常感谢,那么如何输出此时的序号(就是它是这1000个数加到第几项的时候得到的)

论坛优秀回答者

权威

8043 麦片

财富积分



5

主题

9478

帖子

1828

最佳答案
  • 关注者: 417
发表于 2018-6-11 07:45:03 | 显示全部楼层
邢打铁 发表于 2018-6-11 07:38
兄弟,非常感谢,那么如何输出此时的序号(就是它是这1000个数加到第几项的时候得到的) ...

你没试吗?find出来的就是序号,你赋值给N不就行了。

新手

11 麦片

财富积分


050


7

主题

31

帖子

0

最佳答案
 楼主| 发表于 2018-6-11 07:47:19 | 显示全部楼层
本帖最后由 邢打铁 于 2018-6-11 07:49 编辑
jingzhaos 发表于 2018-6-11 07:45
你没试吗?find出来的就是序号,你赋值给N不就行了。

试了,是个空矩阵

如图

如图

论坛优秀回答者

权威

8043 麦片

财富积分



5

主题

9478

帖子

1828

最佳答案
  • 关注者: 417
发表于 2018-6-11 07:59:38 | 显示全部楼层
本帖最后由 jingzhaos 于 2018-6-11 08:00 编辑
邢打铁 发表于 2018-6-11 07:47
试了,是个空矩阵

空阵?按照你的数,当然是空阵,前面已经说了,没有这个数
快照1.jpg



新手

11 麦片

财富积分


050


7

主题

31

帖子

0

最佳答案
 楼主| 发表于 2018-6-11 08:05:33 | 显示全部楼层
jingzhaos 发表于 2018-6-11 07:59
空阵?按照你的数,当然是空阵,前面已经说了,没有这个数

了解了,谢谢大哥,大哥洪福齐天
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表