代码部分程序:
%用ARMA模型仿真风速模型
clc;clear;close all;
%初始化
I=0.16;%湍流强度
L=600;%湍流尺度
average_v=8.11850232757234;%平均风速
Ts=1;%采样时间
T=720;%模拟时间
n=T/Ts;%采样个数
sigma_u=3.48648820169842;%风速序列的标准差
sigma_noise=0.2;%白噪声的方差
% u=yyx_noice_n(0,sigma_noise,n);%产生n个序列的白噪声,均值为0,方差为1;
u=normrnd(0,sigma_noise,n);
v=zeros(1,n);
S=zeros(1,n);
wind=zeros(2,n);
df=1/T;%频率分辨率
f=(0:n-1)*df;
length_f=length(f);
for i=1:length_f
S(i)=4*(sigma_u)^2*L/average_v/(1+70.8*(f(i)*L/average_v)^2)^(5/6);%Von Karman谱功率密度函数
end
R=real(ifft(S));%自相关与功率谱密度是傅立叶变换关系,求自相关
%求解AR模型参数
a(1)=R(2)/R(1);
a(2)=R(3)/R(1)-R(2)*R(2)/R(1)/R(1);
%求解MA模型参数
Ry=(-a(1)+a(2)*a(1))*R(1)+(1+a(1)^2+a(2)^2)*R(2)+(-a(1)+a(1)*a(2))*R(3)+(-a(2))*R(4);
b=Ry/sigma_noise;%beta=Ry/白噪声的方差
结果:


每次运行输出的风速都不一样!!!!!!!!!!!!!!!!
完整程序:
利用白噪声特性和ARMR模型,通过求解模型参数生成随机风速-电子商务文档类资源-CSDN文库
最后
以上就是文静过客最近收集整理的关于利用白噪声特性和ARMR模型,通过求解模型参数生成随机风速的全部内容,更多相关利用白噪声特性和ARMR模型内容请搜索靠谱客的其他文章。
发表评论 取消回复