微信:HuangL1121
QQ号:1224848052
团队:MATLAB代做|MATLAB专业代做|python代做|matlab代写|FPGA代做|FPGA专业代做|MATLAB项目合作|FPGA项目合作|matlab代码|matlab算法源码|python代做|C代做|C++代做
地址:北京市-朝阳区-双柳北街18号院
当前位置:首页 > 技术资讯 > MATLAB代做|FPGA代做|simulink代做——在MATLAB上用三角级数法对路面谱进行仿真,实现不同级数的路面谱的再现

MATLAB代做|FPGA代做|simulink代做——在MATLAB上用三角级数法对路面谱进行仿真,实现不同级数的路面谱的再现

Gq0=256e-6; 
n0=0.1;    % �ο��ռ�Ƶ��n0 
f_min=0.3; 
f_max=15; 
N=300;                     % �������� 
delta_f=(f_max-f_min)/N;    %����Ƶ�ʼ�� 
T=1/delta_f;               %�ܲ���ʱ�� 
V=20;                   %���� 
 
time=linspace(0,T,N);     %ʱ������������ 
 
 
sita=unifrnd(0,2*pi,1,N);    % ����[0,2*pi]������� 
 
 
f=zeros(size(sita)); 
w=zeros(size(sita)); 
Gq=zeros(size(sita)); 
A=zeros(size(sita)); 
Z=zeros(size(sita)); 
 
for i=1:N 
    f(i)=f_min+(i-1/2)*delta_f;  % ��i��С����Ƶ�� 
    w(i)=2*pi*delta_f*(i-1/2);         % ��Ƶ�ʣ�w=2*pi*f   
  
end 
 
for l=1:length(time) 
    for i=1:N 
        Gq(i)=Gq0.*(f(i)./n0).^(-2)*V;    %Gq(f)=Gq0.*(f/n0).^(-w)*V.^(w-1),����w=2; 
        A(i)=sqrt(2)*sqrt(Gq(i)*delta_f)*sin(w(i)*time(l)+sita(i)); 
    end 
    Z(l)=sum(A); 
end 
 
  
figure(1) 
subplot(211); 
plot(time,Z,'r');    % ʱ���ź� 
xlabel('ʱ��t');ylabel('·�治ƽ��Z');grid on;   
 
subplot(212); 
loglog(f,Gq,'r');    % �����׺��� 
xlabel('Ƶ�� Hz');ylabel('�������ܶ� ');hold on;grid on; 
y=(abs(fft(Z))).^2*delta_t/N; %ʱ���ź�FFT�任�õ��Ĺ����׺���

loglog(f,y,'b');

联系:highspeedlogic

QQ :1224848052

微信:HuangL1121

邮箱:1224848052@qq.com

网站:http://www.mat7lab.com/

网站:http://www.hslogic.com/

微信扫一扫:


上一篇:MATLAB代做|FPGA代做|simulink代做——启发式算法—蝙蝠算法

版权所有: MATLAB代做|MATLAB专业代做|python代做|matlab代写|FPGA代做|FPGA专业代做|MATLAB项目合作|FPGA项目合作|matlab代码|matlab算法源码|python代做|C代做|C++代做 Copyright © 2021 All rights reserved.
微信:HuangL1121 QQ:1224848052
地址:北京市-朝阳区-双柳北街18号院

扫一扫,关注我们