版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

風(fēng)雪孤客

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 507.8
帖子: 29
在線: 1.2小時
蟲號: 550327
注冊: 2008-04-23
性別: GG
專業(yè): 信號理論與信號處理

[求助] 請教蘭姆波頻散曲線問題

學(xué)過固體的超聲波的高手，幫幫忙，你們誰有蘭姆波頻散曲線計算的程序，能夠傳給我看看嗎？很著急。dzxukai@sina.cn 謝謝！

非高級版內(nèi)容，請勿放入高級版中

[ Last edited by 華麗的飄過 on 2012-9-22 at 21:55 ]

回復(fù)此樓

» 猜你喜歡

286求調(diào)劑已經(jīng)有9人回復(fù)
一志愿南昌大學(xué)324求調(diào)劑已經(jīng)有6人回復(fù)
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可已經(jīng)有10人回復(fù)
291求調(diào)劑已經(jīng)有5人回復(fù)
308求調(diào)劑已經(jīng)有7人回復(fù)
085404求調(diào)劑，總分309，本科經(jīng)歷較為豐富已經(jīng)有4人回復(fù)
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！已經(jīng)有3人回復(fù)
一志愿中南大學(xué)化學(xué)0703總分337求調(diào)劑已經(jīng)有4人回復(fù)
083000學(xué)碩274求調(diào)劑已經(jīng)有6人回復(fù)
304求調(diào)劑已經(jīng)有4人回復(fù)

1樓 2011-06-28 21:06:40

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wilde2540

木蟲 (正式寫手)

應(yīng)助: 3 (幼兒園)
金幣: 3736.8
散金: 7
紅花: 4
帖子: 590
在線: 188.9小時
蟲號: 835104
注冊: 2009-08-27
專業(yè): 聲學(xué)

【答案】應(yīng)助回帖

引用回帖:

1樓: Originally posted by 風(fēng)雪孤客 at 2011-06-28 21:06:40:
學(xué)過固體的超聲波的高手，幫幫忙，你們誰有蘭姆波頻散曲線計算的程序，能夠傳給我看看嗎？很著急。dzxukai@sina.cn 謝謝！

#include
#include
#include
#include
#include
#include
#include

#define PI 3.1415926
// double cd=5.8e3,ct=3.1e3; //鋼
// double cd=6.42e3,ct=3.04e3; //鋁
double cd=2.68e3,ct=1.1e3; //有機玻璃

int BisectRoot(double a,double b,double h,double eps,double fb,double* x,int n,int *m);
double Func(double cp,double fb);

void main()
{
int i,j,n=20,m,N=300;
double a,b,h,eps,*x,*cp,*dcp,*cg,fb,D=1e-3,step=0.02e3;
FILE *fp;
char str[10],str1[10],cgfile[50][10];
printf("Input the cg file name,please:\n"

;
scanf("%s",str1);

for(j=0;j {
sprintf(str,"%d.dat",j);
strcpy(cgfile[j],str1);
strcat(cgfile[j],str);
}

a=1.0e0;
b=15.0e3;
h=0.7;
eps=1e-6;

for (j=0;j<10;j++)
{
if((fp=fopen(cgfile[j],"w"

)==NULL)
{ printf("cg file %s can't open.\n",cgfile[j]); exit(0);}
else printf("file %s is written.\n",cgfile[j]);

      cp=(double*)calloc(N,sizeof(double));
if(cp==NULL) exit(1);
dcp=(double*)calloc(N,sizeof(double));
if(dcp==NULL) exit(1);
      cg=(double*)calloc(N,sizeof(double));
if(cg==NULL) exit(1);

m=0; i=0; fb=0.0;
do
{
x=(double*)calloc(n,sizeof(double));
         if(x==NULL) exit(1);

        BisectRoot(a,b,h,eps,fb,x,n,&m);
cp=x[j];

if(cp>0 && i>2)
{
dcp=(cp-cp[i-1])/step;
cg=cp/(1.0-fb*dcp/cp);
   fprintf(fp,"%6.3f %6.3f\n ",2*(fb-step)*D, cg*D); //fb*D is frequency with unit kHz
                                                         // c[j] is velocity with unit m/s
                                                         // and c[j]*D is velocity with unit km/s

}
         i++;
        fb+=step;          // fb is frequency with unit Hz, step=50
         free(x);
}while(i
   printf("m=%d\n\n",m);
   fclose(fp);
   free(cp);free(dcp);free(cg);
}
}

double Func(double cp,double fb)
{
double ab,bt,kk,at,kd,kt;
double fs,fa;
ab=sqrt(fabs(cp*cp/(cd*cd)-1.0));
bt=sqrt(fabs(cp*cp/(ct*ct)-1.0));
kk=(cp*cp/(ct*ct)-2.0)*(cp*cp/(ct*ct)-2.0);
at=4*ab*bt;
kd=ab*PI*2*fb/cp;
kt=bt*PI*2*fb/cp;

if(cp>=cd)
{
fs=kk*sin(kt)*cos(kd)+at*sin(kd)*cos(kt);
fa=kk*sin(kd)*cos(kt)+at*sin(kt)*cos(kd);
}
if(cp>ct && cp {
      fs=kk*sin(kt)*cosh(kd)-at*sinh(kd)*cos(kt);
fa=kk*sinh(kd)*cos(kt)+at*sin(kt)*cosh(kd);
}
if(cp {
fs=kk*tanh(kt)-at*tanh(kd);
fa=kk*tanh(kd)-at*tanh(kt);
}
// return (fs);          // 對稱模式
return (fa);          // 反對稱模式
}

int BisectRoot(double a,double b,double h,double eps,double fb,double* x,int n,int *m)
{
double z,z0,z1,y,y0,y1;

*m=0;
z=a;
y=Func(z,fb);
while(1)
{
if((z>b+h/2)||(*m==n))
return(1);

if(fabs(y) {
*m+=1;
x[*m-1]=z;
z+=h/2;
y=Func(z,fb);
continue;
}

z1=z+h;
y1=Func(z1,fb);
if(fabs(y1) {
         *m+=1;
x[*m-1]=z1;
z=z1+h/2;
y=Func(z,fb);
continue;
}

   if(y*y1>0)
{
y=y1;
z=z1;
continue;
}

while(1)
{
if(fabs(z1-z) {
            *m+=1;
   x[*m-1]=(z1+z)/2;
   z=z1+h/2;
   y=Func(z,fb);
   break;
}

z0=(z1+z)/2;
y0=Func(z0,fb);
if(fabs(y0) {
            *m+=1;
   x[*m-1]=z0;
   z=z0+h/2;
   y=Func(z,fb);
   break;
}

   if(y*y0<0)
{
z1=z0;
y1=y0;
}
else
{
z=z0;
y=y0;
}
}
}
}

N年之前的代碼，不負責(zé)對程序進行解釋，僅供參考喔。圖片也是N之前計算的一個例子

贊一下

回復(fù)此樓

2樓2011-10-31 15:21:11