应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求平均值程序出错，麻烦给我指导下。
51/1返回列表
查看: 748  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

涛涛屈

银虫 (初入文坛)

[求助] 求平均值程序出错，麻烦给我指导下。 已有1人参与

想求出一个周期fp=(f+m*a)/cos(th3)平均值，这么写程序可以吗？麻烦给指导下！谢谢了！

clear all
clc
r=200;
l=400;
e=100;
wr=10;
f=1000;
m=20;
t=0.2*pi*rand(1,1e5);
th1=wr*t;
th3=asin((r*sin(th1)+e)/l)
wl=(r*wr*cos(th1))/(l*cos(th3));
a=(r*wr.^2*cos(th1+th3)+l*wl.^2)/cos(th3);
fp=(f+m*a)/cos(th3);
mean(fp)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-01-25 16:04:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涛涛屈

银虫 (初入文坛)
应改可以吧，大家都回去过年了吗？呵呵
赞一下 回复此楼
2楼2015-01-25 16:52:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
clear all
clc
r=200;
l=400;
e=100;
wr=10;
f=1000;
m=20;
t=0.2*pi*rand(1,1e5);
th1=wr*t;
th3=asin((r*sin(th1)+e)/l);
wl=(r*wr*cos(th1))./(l*cos(th3));
a=(r*wr^2*cos(th1+th3)+l*wl.^2)./cos(th3);
fp=(f+m*a)./cos(th3);
mean(fp)
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

涛涛屈
3楼2015-01-25 16:53:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

涛涛屈

银虫 (初入文坛)
送红花一朵
引用回帖:
3楼: Originally posted by CelestialCYJ at 2015-01-25 16:53:55
clear all
clc
r=200;
l=400;
e=100;
wr=10;
f=1000;
m=20;
t=0.2*pi*rand(1,1e5);
th1=wr*t;
th3=asin((r*sin(th1)+e)/l);
wl=(r*wr*cos(th1))./(l*cos(th3));
a=(r*wr^2*cos(th1+th3)+l*wl.^2)./cos( ...

感觉不对了，因为随机抽样出来的，cos（th3）有可能为0.分母出现0了就不对了。知道怎么搞这种问题的平均值吗？
赞一下 回复此楼
4楼2015-01-25 17:09:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

CelestialCYJ

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
涛涛屈: 金币+8, ★★★很有帮助 2015-01-26 09:04:55
引用回帖:
4楼: Originally posted by 涛涛屈 at 2015-01-25 17:09:19
感觉不对了，因为随机抽样出来的，cos（th3）有可能为0.分母出现0了就不对了。知道怎么搞这种问题的平均值吗？...

分母加个小量eps
赞一下(1人) 回复此楼
5楼2015-01-25 19:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 涛涛屈 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 05:47 by xzcv0Tehd8lQ
[论文投稿] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 4/200 2026-11-02 03:50 by xzcv0Tehd8lQ
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 03:37 by xzcv0Tehd8lQ
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +5 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:05 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:24 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:05 by xzcv0Tehd8lQ
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 16:58 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 16:52 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 16:48 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
信息提示
关闭
请填处理意见
关闭