jieguo.m

来自「基于瑞利波及lamb波的频散曲线程序,可用于根据实验绘制Rali波及lamb波的」· M 代码 · 共 32 行

y0=[1.68206137682026 1.68194092156440 1.68157872834317 1.68097244155653 1.68010718889650 1.67898296899564 1.67758685742774 1.67589959135014... 
    1.67390299570862 1.67157766759147 1.66889356050644 1.66581731135091 1.64872482121357 1.61153663679823 1.54549533597046 1.44474073262713...
    1.23794972767620 1.10500968400776 1.03076557525075];
y1=[6283.186233486030 62.92312980782626 28.90912239606697 19.18165211154748 14.30475979861182 11.37340699124723 9.41790903692696 8.02184193479668 ...
    6.97651320715218 6.16571708128425 5.51958426156915 4.99357330343271 3.61481766588039 2.72074405077632 2.25442978143520 2.02046354627868 ...
    1.85206761772817 1.78525945284236 1.72041475253056];
y2=[12566.37064693027 125.6669639881980 31.59129127799133 21.19246565155007 16.01842799877972 12.92905106513531 10.87897897111140...
     9.42106083900265 8.33221539307083 7.48878797293523 6.81674603332060 6.26907240393502 4.82187912085487 3.84447448566908 3.26958104814514...
     2.89327925841040 2.43247897191065 2.13378948924025 1.99603201315627];
y3=[10^5 10^4 62.83837316912972 41.89769320789427 31.42900195797974 25.14911848311177 20.96365202755169 17.97500900723160 15.73439475384518...
    13.99248075995935 12.59967157479745 11.46077157335096 8.42888036266856 6.35405850214883 5.11903355569367 4.30535752642957...
     3.31525515187652 2.75751740711171 2.42124208912105];
x=[0.001 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.5 2.0 2.5 3.0 4.0 5.0 6.0];
plot(x,y0);
hold on
plot(x,y1);
plot(x,y2);
plot(x,y3);
grid on
xlabel('K*H');
ylabel('C/Cs');

%结果验证
%kh=3.0;p=sqrt(1.44484073262713^2-1)
%cd=6100;cs=3300;m=cs/cd;
%y1=-4*tan(kh/2*((p^2+1)*m^2-1)^(1/2))*p*((p^2+1)*m^2-1)^(1/2);
%y2=(p^2-1)^2*tan(kh/2*p);
%y=y1-y2
% 数据导入
%x=input('jieguo0=')
% x=a(:,1);y0=a(:,2);y1=a(:,3);y2=a(:,4);y3=a(:,5);y4=a(:,6);y5=a(:,7);plot(x,y0);
% hold on plot(x,y1);plot(x,y2);plot(x,y3);plot(x,y4);plot(x,y5);