m基于遗传优化的时域声辐射模态的振动控制算法的matlab仿真

我爱C编程

已于 2023-02-18 17:34:17 修改

阅读量492

点赞数 1

分类专栏： Matlab优化算法 Matlab控制器文章标签： matlab 遗传优化时域声辐射模态振动控制

于 2023-02-18 16:33:40 首次发布

本文链接：https://blog.csdn.net/hlayumi1234567/article/details/129101343

版权

Matlab优化算法同时被 2 个专栏收录

116 篇文章 35 订阅

订阅专栏

Matlab控制器

27 篇文章 7 订阅

订阅专栏

1.算法仿真效果

matlab2013b仿真结果如下：

2.MATLAB核心程序

...........................................................................
%仿真时域声辐射模态平板的区域范围
%仿真时域声辐射模态平板的区域范围
Lx      = 0.3;
Ly      = 0.3;
%平面厚度
hs      = 0.01;
%材料的密度
density = 5900;
%结构刚度
Es      = 1e11;
As      = Lx*hs;
MOM     = hs^3*sqrt(Lx^2+Ly^2)/24;
%阻尼比
Rb      = 0.01;
%这两个值设置较小即可，否则运算量非常大
m       = 2;
n       = 2;
%平板区域分割数
MN      = m*n;
%力的区域
Xp        = [0.15,0.2];
Xcon      = [0.5 ,0.4];
Num_Force = 1;
Ms        = density*Lx*Ly*As*eye(MN);
for j = 1:MN
    a1(j) = j/5;
    b1(j) = (sinh(a1(j)) + sin(a1(j)))/(cosh(a1(j)) - cos(a1(j)));
end
Ks        = Es*MOM*Lx*diag(a1(1:MN).^2,0);
for k = 1:MN
    Fvalue(k,1:Num_Force) = b1(k)*(cos(a1(k)*mean(Xcon)));
    Pvalue(k,1)           = b1(k)*(cos(a1(k)*mean(Xp))); 
end
......................................................................

for j=1:1:NIND 
    x(j,1)  = phen(j,1);
    deltaT  = 0.0005;
    Tall    = 2000;
    T       = deltaT:deltaT:deltaT*Tall; 
    %定义一个控制输入
    %控制力频率
    f       = 50*(1+5*x(j,1));%50~200Hz之间优化
    U       = sin(2*pi*f*[1:length(T)]*deltaT); 
    [Y,X]   = lsim(Ac,Bc,Cc,Dc,U,T); 
    for t = 1:Tall
        Ws(t)= 1*(abs(Y(t)))^2; 
        tmps(t) = 10*log10(Ws(t)/Wr);
        if abs(tmps(t)) >1e10
           tmps(t) =  -1e6;
        end
    end
    J(j,1)  = min(tmps);
end

Objv  = 1./(J+eps);
gen   = 0; 

while gen < MAXGEN;   
      gen
      FitnV=ranking(Objv);    
      Selch=select('sus',Chrom,FitnV);    
      Selch=recombin('xovsp', Selch,0.9);   
      Selch=mut( Selch,0.01);   
      phen1=bs2rv(Selch,FieldD);   
      for a=1:1:NIND  
          x1(a,1)=phen1(a,1);              
            deltaT  = 0.0005;
            Tall    = 2000;
            T       = deltaT:deltaT:deltaT*Tall; 
            %定义一个控制输入
...............................................................
            JJ(j,1)  = min(tmps);
      end 
      Objvsel=1./(JJ+eps);    
      [Chrom,Objv]=reins(Chrom,Selch,1,1,Objv,Objvsel);   
      gen=gen+1; 
end 

MIN=min(Objv); 
for ttt=1:1:size(Objv)     
    if Objv(ttt)<=MIN         
       tt=ttt;         
       break;     
    end
end                 

deltaT  = 0.0005;
Tall    = 2000;
T       = deltaT:deltaT:deltaT*Tall; 
for t = 1:Tall
    f(t) = (t-1)/Tall*50*(1+5*x1(tt,1));  
end
U     = sin(2*pi*f.*[1:Tall]*deltaT);
[Y,X] = lsim(Ac,Bc,Cc,Dc,U,T); 

figure; 
plot(T,Y); 
xlabel('time(s)');
ylabel('y1');
grid on
title('第1阶伴随系数');
axis([0,1,-20,20]);
for t = 1:Tall
    Ws(t)= 1*(abs(Y(t)))^2; 
end
 
figure;
plot(deltaT*(1:Tall),10*log10(Ws/Wr),'b');
xlabel('time(s)');
ylabel('声功率级别db');
grid on
axis([0,1,20,130]);
save r3.mat Y Ws Wr
08_020_m