| 查看: 679 | 回复: 5 | |||
[求助]
为什么这两个程序差别这么大，只是将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开始的，明明取得向量不是从零开始的？？？  untitled.jpg  1.jpg |
回复此楼» 猜你喜欢
情人节自我反思：在爱情中有过遗憾吗？
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有7人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
-
过年走亲戚时感受到了所开私家车的鄙视链
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 超右脑照相记忆法，右脑中奇异的学习回路。【转载】
已经有71人回复
- 分子中电荷变化种种
已经有82人回复
- 关于towhee 回答zhang_jaj的问题
已经有44人回复
- 大家帮忙看看，双向，电压升高模式为什么是这个样子（先很慢，后很快）
已经有1人回复
- 【资源】催化剂制备中浸渍（等体积浸渍）不断改进
已经有151人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
- 【zhou2009个人文集】极性键电子密度差Δρ的解析
已经有15人回复
- 赴美HCII 2009会议回顾
已经有70人回复
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
2楼2015-11-06 13:22:15
3楼2015-11-06 14:37:51
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
4楼2015-11-06 15:05:06
5楼2015-11-06 17:04:46
mygt_hit
专家顾问 (职业作家)
-
专家经验: +362 - 应助: 438 (硕士)
- 贵宾: 0.019
- 金币: 19571.1
- 散金: 5130
- 红花: 135
- 沙发: 3
- 帖子: 4938
- 在线: 990.5小时
- 虫号: 1489764
- 注册: 2011-11-13
- 性别: GG
- 专业: 结构工程
- 管辖: 土木建筑
6楼2015-11-06 17:42:35