基于Davenport风速谱两点时程模拟
转自:http://blog.sina.com.cn/s/blog_a1d5b9ba0102wfrm.html
clear all;
clc;
N=500; %采样点数
omegaup=4*2*pi; %截断频率
dm=omegaup/N; %频率步长
dt=3.2*pi/(2*omegaup); %时间步长0.2
k=0.00464; %地面粗糙度系数地面粗糙度等级A B C D:K= 0.00129、0.00215、0.00464、0.01291
d=0.001;
f=d:d:10; %时间从0.001到10s,步进值为0.001
v10=26.8; %设计风速 26.8m/s——50年一遇十分钟平均风速最大值
x=1200*f/v10;
s=4*k*v10*v10.*x.^2./f./(1+x.^2).^(4/3); �venport风速谱
%x=1200*f/v10;
%Sv_Davenport=4*x.^2./(6*(1+x.^2).^(4/3)); �venport风速谱
%x1=z*f/(v10*(z/10)^2);
%Sv_Simiu=200*x1./(6*(1+50*x1).^(5/3)); %Simiu风速谱
%Sv_Kaimal=200*fn./(6*(1+50*fn).^(5/3)); %Kaimal风速谱
%x2=1800*f/v10;
%Sv_Harris=4*x2./(6.677*(2+x2.^2).^(5/6)); %Harris风速谱
%Lu=100*(z/30)^(1/2);
%x3=Lu*f/v0;
%Sv_Karman=4*x3./((1+70.8*x3.^2).^(5/6)); %Karman风速谱
z1=10; %取第一点为10米高度
z2=50; %取第二点为50米高度
r=0.22; %考虑地表粗糙度影响的无量纲幂指数,按中国规范取0.22-c类
v5=v10*(z2/z1)^(r); %计算n米高处的平均风速——52.8m处平均风速期望值
Cz=10; %Cz表示z方向的指数衰减系数(取经验值)
v1=zeros(2*N,1); %产生一个全零矩阵
v2=zeros(2*N,1); %产生一个全零矩阵
thta1=rand(N,1); %随机矩阵
thta2=rand(N,1); %随机矩阵
node=1;
for K=1:node
for j=1:2*N
sum1=0;
sum2=0;
for l=1:N
m1=l*dm-0.5*dm;
m2=l*dm;
x1=1200*m1/(2*pi*v10);
s11=2*pi*4*k*v10*v10*x1*x1./m1./(1+x1*x1).^(4/3); %1点自互功率谱
x2=1200*m2/(2*pi*v10);
s22=2*pi*4*k*v5*v5*x2*x2./m2./(1+x2*x2).^(4/3); %2点自互功率谱
s12=sqrt(s11*s22).*exp(-2*m2*Cz*abs(z1-z2)./(2*pi*(v10+v5))); %仅考虑竖向相关性互功率谱
s21=sqrt(s11*s22).*exp(-2*m1*Cz*abs(z1-z2)./(2*pi*(v10+v5))); %仅考虑竖向相关性互功率谱
S=[s11 s12;s21 s22]; %谱矩阵
H=chol(S); %丘拉斯基分解-因式分解
a1=abs(H(1,1));
H1=H';
a21=abs(H1(2,1));
a22=abs(H1(2,2));
b1=cos((m1*dt*(j-1))+2*pi*thta1(l,1));
b2=cos((m2*dt*(j-1))+2*pi*thta2(l,1));
c1=a1*b1;
c21=a21*b1;
c22=a22*b2;
d1=(dm).^0.5*c1;
d2=(dm).^0.5*(c21+c22);
sum1=sum1+d1;
sum2=sum2+d2;
end
sum1=0.8*sum1;
sum2=0.8*sum2;
v1(j,K)=sum1;
v2(j,K)=sum2;
end
end
u1=v1+ v5;
u2=v2+ v5;
t= (0:2*N-1)*dt;
subplot(2,2,1);
plot(t,u1,'b-');
xlabel('t(s)');
ylabel('v(t)');
axis([-1 180 0 80]);
subplot(2,2,3);
plot(t,u2,'r-');
xlabel('t(s)');
ylabel('v(t)');
axis([-1 180 0 80]);
Y=fft(v2); %对数值解作傅立叶变换
Y(1)=[]; %去掉零频量
m=length(Y)/2; %计算频率个数;
power=150*abs(Y(1:m)).^2/(length(Y).^2); %计算功率谱
freq=5*(1:m)/length(Y); %计算频率,因为时间步长为0.125,而不是1,故乘以8
subplot(2,2,2);
loglog(freq,power,'r-',f,s,'g-'); %对数显示,比较plot(t,v2,'r-');
axis([-50 15 -10 1000]); %画出坐标轴比例axis([xmin xmax ymin ymax])
set(gca,'xtick',[0.1 1 10]); %自动定义刻度
set(gca,'ytick',[0.1 1 10]);
grid on
xlabel('频率');
ylabel('功率');
基于Davenport风速谱两点时程模拟相关推荐
- 基于Davenport风速谱单点时程模拟
转自:http://blog.sina.com.cn/s/blog_a1d5b9ba0102wfrw.html %%%%%%%%%%%%%%%%main procedure%%%%%%%%%%%%%% ...
- 脉动风速谱matlab程序,1.脉动风速谱.PDF
1. 脉动风速谱 根据统计分析可知,任意一点的风速 V ( z, t) 为平稳Gauss随机过程,可以分成平均风速 [5] v (z) 与零均值脉动风速v(z, t) .平均风速的高度变化律可以采用对 ...
- java+postgis实现根据两点生成模拟轨迹gps数据
java+postgis实现根据两点生成模拟轨迹gps数据 文章目录 java+postgis实现根据两点生成模拟轨迹gps数据 前言 一.实现流程 1.请求参数 2.功能流程 3.postgis重要 ...
- Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动
Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...
- 不相关子查询的工作方式是_课题组工作|Nucleic Acids Research|基于表达密度谱的特征子空间分离及相关单细胞转录组分群新算法...
大家好!为大家分享本课题组近期发表在Nucleic Acids Research的文章,题目为 "Entropy subspace separation-based clustering f ...
- 四六级考试中的计算机类文章,毕业设计论文-基于web的英语四六级模拟在线考试系统.doc...
毕业设计论文-基于web的英语四六级模拟在线考试系统 基于web的英语四六级模拟在线考试系统 摘要 随着科学技术的不断提高,计算机科学技术日渐成熟,其强大的功能已为人们深刻认识,它已进入人类社会的各个 ...
- 51单片机 外部时钟_基于51单片机的LCD12864显示模拟时钟
今天要分享的仿真和程序是基于51单片机的LCD12864显示模拟时钟.公众号之前也分享过数字时钟,但是那些时钟都是采用数字显示的形式来显示时间,而今天分享一个采用LCD12864显示的模拟指针式时钟. ...
- 基于共识算法和区块链模拟实现超级账本
基于共识算法和区块链模拟实现超级账本 实验语言:GO 实验环境:GoLand 2022.1 : go 1.13.4.widows-amd64.msi :curl-7.83.1 实验中使用的开源包: h ...
- matlab中eacf函数,基于MATLAB的切比雪夫I型模拟低通滤波器设计
基于MATLAB的切比雪夫I型模拟低通滤波器设计课程设计名称:数字信号处理课程设计专业班级 : 电信0604 学生姓 名 : 学号 : 20064300430 指导教师 : 课程设计时间:2009.6 ...
最新文章
- 教程 | OpenCV深度神经网络实现人体姿态评估
- 数据结构(C语言实现)勘误
- 图分区技术基本概念【1】
- Python-jsonpath使用和json转换
- re匹配正则字符串中的起始和结束元字符的使用方法
- mysql使用jtable_jtable 的简单使用
- 基于TableStore的物联网元数据管理
- boost.asio基础篇 小白入门注解
- Linux的sys接口和proc接口创建
- CSE 110: Principles of Programming Languages
- 人脸识别全国产化电子元件推荐方案
- 浅谈Linux管线命令grep
- 【神经网络】权重衰减(weight-decay)
- 拉流与推流的区别_RTSP协议交互流程之推流分析
- 【小白笔记】目标跟踪Deep Regression Tracking with Shrinkage Loss
- python的roc曲线与阈值_python 使用sklearn绘制roc曲线选取合适的分类阈值
- Centos如何联网
- Skiplist(跳表)实现
- 盘点世界十大著名黑客攻击事件
- 学组合数学心得与题解(一)——组合计数
热门文章
- 更新系统后iphone6s无服务器,目瞪口呆!iPhone6S升级iOS11后,居然秒变iPhoneX
- vrrp路由器三种状态_虚拟路由器冗余协议(VRRP)
- 爱奇艺2015校园招聘产品笔试题
- PGN与SPN的区别
- win7 xp双系统启动引导修复指令
- YOLOv7: Trainable bag-of-freebies sets new state-of-the-art for real-time object detectors
- 2021年安全员-B证考试题及安全员-B证模拟试题
- uni-app使用setStorageSync存储历史搜索记录,与setStorage区别
- 相干波和杨氏双缝实验(大学物理笔记)
- 基于Java的物流快递行业仓储管理系统