24小时热门版块排行榜

返回列表

1375481653

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 736.5
散金: 12
帖子: 133
在线: 105.2小时
虫号: 3439317
注册: 2014-09-24
专业: 机械结构强度学

[求助] 计算仿真已有1人参与

程序出错，高手帮忙指点下！！！谢谢！！公式见附件。
clear;
syms x Z;
w=0:2:10;
R=zeros(1,length(w));
for i=1:length(w)
  for a=1:10
   Z=Z*0.5*(1+erf((200-(a+1)/75)/(40*sqrt(2))));
  end
  y=int(((1/(60*sqrt(2*pi)))*exp(-(x-600)^2/7200))*Z,x,-inf,inf);
  R(i)=vpa(y,10);
end
plot(w,R,'r.')

回复此楼

» 本帖附件资源列表

欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
附件 1 : 新建_Microsoft_Word_文档_(2).doc

2016-01-23 14:58:37, 19 K

» 猜你喜欢

1楼 2016-01-23 14:58:40

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wuyu__

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 245.6
红花: 1
帖子: 108
在线: 38.2小时
虫号: 1841336
注册: 2012-05-31
性别: GG
专业: 固体力学

【答案】应助回帖

感谢参与，应助指数 +1

word文档打不开，
...
R(i)=y;
end
R=double('R');
plot(w,R,'r.')
没有类型错误，但是出不了结果，算法有误

赞一下

回复此楼

Good Good Study,Day Day UP

2楼2016-01-23 23:36:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wuyu__

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 245.6
红花: 1
帖子: 108
在线: 38.2小时
虫号: 1841336
注册: 2012-05-31
性别: GG
专业: 固体力学

【答案】应助回帖

syms Z x
w=0:2:10;
R=zeros(1,length(w))*sym('1');
for i=1:length(w)
  for a=1:10
   Z=Z.*0.5.*(1+erf((200-(a+1)/75)/(40.*sqrt(2))));
  end
R(i)=int(((1/(60.*sqrt(2*pi))).*exp(-(x-600)^2/7200)).*Z,x,-inf,inf);
end
R=double('R');
plot(w,R,'r.')

赞一下(1人)

回复此楼

Good Good Study,Day Day UP

3楼2016-01-23 23:39:50

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖