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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 为什么这两个程序差别这么大，只是将Omega的范围变化了，朋友看看怎么回事？？？？
61/1返回列表
查看: 861  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

阿毛虫

银虫 (正式写手)

[求助] 为什么这两个程序差别这么大，只是将Omega的范围变化了，朋友看看怎么回事？？？？ 已有1人参与

第一幅图的程序是：
clear all;clc
G=8e10;K=1.2; E=2.06e11; P=7.9e3; L=4; D=0.8; H=0.2; k=1/K; I=D*(H^3)/12; A=D*H;
  Omega2=linspace(sqrt((k*G*A)/(P*I)),inf,10000);
      b2=sqrt(P*A/(E*I)).*Omega2*L^2;    r=sqrt(I/A)/L; s=sqrt(E*I/(k*A*G))/L;
alfa2=sqrt(r^2+s^2-sqrt((r^2-s^2)^2+4./(b2.^2)))/sqrt(2);  
    beta2=sqrt(r^2+s^2+sqrt((r^2-s^2)^2+4./(b2.^2)))/sqrt(2);   
F2=2-2*cos(alfa2.*b2).*cos(beta2.*b2)+(b2./sqrt(-1+b2.^2*r^2*s^2)).*(-r^2+3*s^2+(b2.^2*s^2)*(r^2-s^2)^2).*sin(beta2.*b2).*sin(alfa2.*b2);
figure(2)
plot(Omega2,F2)
第二幅图的程序是
clear all;clc
G=8e10;K=1.2; E=2.06e11; P=7.9e3; L=4; D=0.8; H=0.2; k=1/K; I=D*(H^3)/12; A=D*H;
  Omega2=linspace(sqrt((k*G*A)/(P*I)),2e20,10000);
      b2=sqrt(P*A/(E*I)).*Omega2*L^2;    r=sqrt(I/A)/L; s=sqrt(E*I/(k*A*G))/L;
alfa2=sqrt(r^2+s^2-sqrt((r^2-s^2)^2+4./(b2.^2)))/sqrt(2);  
    beta2=sqrt(r^2+s^2+sqrt((r^2-s^2)^2+4./(b2.^2)))/sqrt(2);   
F2=2-2*cos(alfa2.*b2).*cos(beta2.*b2)+(b2./sqrt(-1+b2.^2*r^2*s^2)).*(-r^2+3*s^2+(b2.^2*s^2)*(r^2-s^2)^2).*sin(beta2.*b2).*sin(alfa2.*b2);
figure(2)
plot(Omega2,F2)
这两个程序只有Omega的范围有变化，一个是到无穷，一个是确定的一个挺大的数，为什么出现的图会有的有图像有的没有，还有一点就是为什么图像的横坐标是从0开始的，明明取得向量不是从零开始的？？？

为什么这两个程序差别这么大，只是将Omega的范围变化了，朋友看看怎么回事？？？？
untitled.jpg


为什么这两个程序差别这么大，只是将Omega的范围变化了，朋友看看怎么回事？？？？-1
1.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
学习吧
1楼 2015-11-06 09:30:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
阿毛虫: 金币+10, ★★★很有帮助, 明白 2015-11-06 14:38:17
1. 第一段代码中，Omega2的内容是：[50315.4605426628, 之后全是inf]，所以画图只画了一个点。
2. 第二段代码中，起点不是0，而是50315.4605426628，只不过横坐标乘以10^20，所以看上去像0了。
赞一下 回复此楼
知其然，知其所以然。
2楼2015-11-06 13:22:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)
引用回帖:
2楼: Originally posted by mygt_hit at 2015-11-06 13:22:15
1. 第一段代码中，Omega2的内容是：，所以画图只画了一个点。
2. 第二段代码中，起点不是0，而是50315.4605426628，只不过横坐标乘以10^20，所以看上去像0了。

谢谢，另外我还想问一下，如果就第二幅图来说，我如果想求零点，应该怎么求，具体一点可以吗？
赞一下 回复此楼
学习吧
3楼2015-11-06 14:37:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 阿毛虫 at 2015-11-06 14:37:51
谢谢，另外我还想问一下，如果就第二幅图来说，我如果想求零点，应该怎么求，具体一点可以吗？...

设置横坐标取值范围在0附近求吧，或者列方程用MATLAB求。
赞一下 回复此楼
知其然，知其所以然。
4楼2015-11-06 15:05:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)
引用回帖:
4楼: Originally posted by mygt_hit at 2015-11-06 15:05:06
设置横坐标取值范围在0附近求吧，或者列方程用MATLAB求。...

老师能具体一点吗？用fsolve？
赞一下 回复此楼
学习吧
5楼2015-11-06 17:04:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 阿毛虫 at 2015-11-06 17:04:46
老师能具体一点吗？用fsolve？...

我就知道个大概，抱歉，fsolve solve等，还有符号工具箱都可以求解吧。
赞一下 回复此楼
知其然，知其所以然。
6楼2015-11-06 17:42:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿毛虫 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 12:40 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 12:38 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-11-01 12:32 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 7/350 2026-11-01 12:27 by BCdvBAEpt3zm
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 11/550 2026-11-01 12:24 by BCdvBAEpt3zm
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 12:14 by BCdvBAEpt3zm
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 11:18 by Equinoxhua
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考研] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:56 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 10:27 by BCdvBAEpt3zm
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 07:01 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
关闭
请填处理意见
关闭