24小时热门版块排行榜

>论坛更新日志 (3662)
>虫友互识 (732)
>导师招生 (320)
>文献求助 (236)
>考博 (210)
>休闲灌水 (206)
>硕博家园 (151)
>论文投稿 (132)
>博后之家 (101)
>考研 (95)
>教师之家 (52)
>招聘信息布告栏 (49)
>精细化工 (49)
>公派出国 (43)
>基金申请 (41)
>找工作 (30)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

阿毛虫

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 89.2
散金: 81
帖子: 750
在线: 113小时
虫号: 3313820
注册: 2014-07-08
性别: GG
专业: 机械动力学

[求助] 关于梁的一个简单的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的值。应该怎么写程序啊？？？？

回复此楼