| 查看: 944 | 回复: 3 | |||
阿毛虫银虫 (正式写手)
|
[求助]
关于梁的一个简单的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的值。应该怎么写程序啊？？？？ |
回复此楼» 猜你喜欢
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
机械工程论文润色/翻译怎么收费?
已经有155人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”学硕与专硕有缺额
已经有0人回复
-
太原科技大学材料科学与工程学院“03高性能轻质合金及新能源材料”方向有调剂名额
已经有0人回复
-
辽宁工程技术大学材料科学与工程、储能科学与工程、材料与化工专业接收调剂，名额充足
已经有0人回复
-
辽宁工程技术大学材料科学与工程接收调剂，目前还有名额，欢迎报考
已经有0人回复
卡尔俾斯
铁杆木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 11486.5
- 散金: 220
- 红花: 6
- 帖子: 2187
- 在线: 247.7小时
- 虫号: 2181524
- 注册: 2012-12-12
- 性别: GG
- 专业: 制造系统与自动化
2楼2015-11-07 07:52:11
阿毛虫
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 136.7
- 散金: 81
- 帖子: 752
- 在线: 113.8小时
- 虫号: 3313820
- 注册: 2014-07-08
- 性别: GG
- 专业: 机械动力学
3楼2015-11-07 15:24:25
卡尔俾斯
铁杆木虫 (著名写手)
- 应助: 10 (幼儿园)
- 金币: 11486.5
- 散金: 220
- 红花: 6
- 帖子: 2187
- 在线: 247.7小时
- 虫号: 2181524
- 注册: 2012-12-12
- 性别: GG
- 专业: 制造系统与自动化
4楼2015-11-07 16:36:49
5