?? lames0.m
字號:
%求解Lame波對稱部分的基本模態So
format long
syms c; %定義速度c為符號函數
kh=5 %k為波數,h為總的板厚
cd=5790;cs=3200; %cd為縱波波速,cs為橫波波速
y1=4*tanh(kh/2*(1-(c/cd)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cs)^2)^(1/2));
ezplot(y1,[0,3500]); %作圖求解
hold on
grid on
ezplot(y2,[0,3500]);
[xx,yy]=ginput(1);
[x(1),fval,exitflag]=fzero(inline(y1-y2),xx(1),[]);
c=x(1)
y1=4*tanh(kh/2*(1-(c/cd)^2)^(1/2))*(1-(c/cd)^2)^(1/2)*(1-(c/cs)^2)^(1/2);
y2=(2-(c/cs)^2)^2*tanh(kh/2*(1-(c/cs)^2)^(1/2));
e1=y1-y2
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -