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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 机械 » 动力机械 » 关于梁的一个简单的matlab程序，我想求出这个方程的零点，求大神指点
41/1返回列表
查看: 1088  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

阿毛虫

银虫 (正式写手)

[求助] 关于梁的一个简单的matlab程序，我想求出这个方程的零点，求大神指点 已有1人参与

global Omega1;
global b1;
global r;
global s;
global alfa1;
global beta1;
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;
Omega1=linspace( 0,sqrt((k*G*A)/(P*I)),10000);
b1=sqrt(P*A/(E*I)).*Omega1;
r=sqrt(I/A);
s=sqrt(E*I/(k*A*G));
alfa1=b1.*sqrt(-(r^2+s^2)+sqrt((r^2-s^2)^2+4./(b1.^2)))/sqrt(2);
beta1=b1.*sqrt(r^2+s^2+sqrt((r^2-s^2)^2+4./(b1.^2)))/sqrt(2);  
  clamped_clamped1=2-2*cosh(alfa1.*L).*cos(beta1.*L)-(b1./sqrt(1-b1.^2*r^2*s^2)).*(r^2-3*s^2-(b1.^2*s^2)*(r^2-s^2)^2).*sin(beta1.*L).*sinh(alfa1.*L);%定义函数表达式
figure(1)
plot(Omega1,clamped_clamped1,'b')
title('clamped_clamped beam1')
hold on
plot(xlim,[0 0],'r')
我就想求出clamped_clamped1=0时候，Omega1的值。应该怎么写程序啊？？？？
回复此楼

» 猜你喜欢
学习吧
1楼 2015-11-06 20:05:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡尔俾斯

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
这个搜根求数值解

发自小木虫Android客户端
赞一下 回复此楼
2楼2015-11-07 07:52:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)
引用回帖:
2楼: Originally posted by 卡尔俾斯 at 2015-11-07 07:52:11
这个搜根求数值解

具体一点把，刚学这个东西
赞一下 回复此楼
学习吧
3楼2015-11-07 15:24:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡尔俾斯

铁杆木虫 (著名写手)

【答案】应助回帖

你找个数值算法的书，看看二分法搜根，可以用二分法搜根

发自小木虫Android客户端
赞一下 回复此楼
4楼2015-11-07 16:36:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 阿毛虫 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭