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代做|FPGA代做|simulink代做——基于桥面不平顺公路梁桥车桥耦合随机振动研究