道路空间功率谱密度与时间功率谱密度(笔记)
0.随机路面
仿真时间200秒
车速15m/s
采样间隔0.25秒(4Hz)
A级路面平均振幅:0.0044米
B级路面平均振幅:0.0087米
C级路面平均振幅:0.0173米
D级路面平均振幅:0.0346米
1.先上代码其中之一
clc
clear
close all
%%
SimTime=200;
dt=0.01;%仿真步长
time=0:dt:SimTime;
sim_step = length(time);
Nt=length(time); % 采样点(可能要修改)
u=10; % m/s
df=1/(Nt*dt); % 采样频率间隔
f=0:df:1/(2*dt); % 采用频率一半(可修改)
nf=length(f); % 采样点的个数
n=f/u; % 采样空间频率,F是浮动的,随车速的改变而改变
%获得标准空间路面谱,A级
Pa=LMP(16,nf,n);
pa=Pa/u; % 转化为频谱
%获得标准空间路面谱,B级
Pb=LMP(64,nf,n);
pb=Pb/u; % 转化为频谱
%获得标准空间路面谱,C级
Pc=LMP(256,nf,n);
pc=Pc/u; % 转化为频谱
%获得标准空间路面谱,D级
Pd=LMP(1024,nf,n);
pd=Pd/u; % 转化为频谱
%%
%1输出标准路面空间频率谱G(n),结果与车速、SimTime、dt无关
figure(1)
loglog(n,Pa); %1输出标准路面空间频率谱G(n)
hold on
loglog(n,Pb);
hold on
loglog(n,Pc);
hold on
loglog(n,Pd);
xlabel('空间频率n[1/m]')
ylabel('路面不平度位移功率谱密度Gq(n)[m^3]')
legend('A','B','C','D')%2输出时间频率谱G(f),结果与车速有关
% figure(2)
% loglog(f,pa)
% hold on
% loglog(f,pb)
% hold on
% loglog(f,pc)
% hold on
% loglog(f,pd)
% xlabel('时间频率[Hz]')
% ylabel('G(f)[m^2/Hz]') % 2输出时间频率谱G(f)
% legend('A','B','C','D')
2.参考链接
链接1
链接2
链接3
3.笔记
Sine_RoadRoughnessGradeClassCal_230713.m中:
3.1. 降SimTime可以降等级,升omiga会使改线条等级范围变大,降车速可降等级
但是在时间-输入关系不变的情况下,低车速对应高空间频率路面,高车速对应低空间频率路面
代码一:SimTime=20/4;dt=0.01;omiga=2pi1.5;A=0.05;V=18;
%%
%1.正弦路面
SimTime=20/4;
dt=0.01;%仿真步长
time=0:dt:SimTime;
sim_step = length(time);
omiga=2*pi*1.5;
A=0.05;%振幅
Z1=A*sin(time*omiga);%路面高程值
Z_p1 =A*cos(time*omiga)*omiga;%路面高程值导数
x=Z1;Nt=length(x); % 采样点(可能要修改)
V=18; % km/h
结果C级路面

注意此图中的标准等级路面PSD-(1/n)函数线条与汽车理论208页图6-5一致
代码一点五:SimTime=20;dt=0.01;omiga=2*pi;A=0.05;V=18;
%%
%1.正弦路面
SimTime=20;
dt=0.01;%仿真步长
time=0:dt:SimTime;
sim_step = length(time);
omiga=2*pi;
A=0.05;%振幅
Z1=A*sin(time*omiga);%路面高程值
Z_p1 =A*cos(time*omiga)*omiga;%路面高程值导数
x=Z1;Nt=length(x); % 采样点(可能要修改)
V=18; % km/h
结果B级路面

代码二:SimTime=20/2;dt=0.01;omiga=2pi2;A=0.019;V=18;
SimTime=20/2;
dt=0.01;%仿真步长
time=0:dt:SimTime;
sim_step = length(time);
omiga=2*pi*2;
A=0.019;
Z1=A*sin(time*omiga);%路面高程值
Z_p1 =A*cos(time*omiga)*omiga;%路面高程值导数
x=Z1;
dt=0.01;
Nt=length(x); % 采样点(可能要修改)
V=18; % km/h
结果A级路面

代码二点五:SimTime=20/4;dt=0.01;omiga=2*pi;A=0.019;V=9;
SimTime=20/4;
dt=0.01;%仿真步长
time=0:dt:SimTime;
sim_step = length(time);
omiga=2*pi;
A=0.019;
Z1=A*sin(time*omiga);%路面高程值
Z_p1 =A*cos(time*omiga)*omiga;%路面高程值导数
x=Z1;
dt=0.01;
Nt=length(x); % 采样点(可能要修改)
V=9; % km/h
结果B级路面

相关文章:
道路空间功率谱密度与时间功率谱密度(笔记)
0.随机路面 仿真时间200秒 车速15m/s 采样间隔0.25秒(4Hz) A级路面平均振幅:0.0044米 B级路面平均振幅:0.0087米 C级路面平均振幅:0.0173米 D级路面平均振幅:0.0346米 1.先上代码其中之一 clc clear clo…...
JMeter接口测试之文件上传
最近用JMeter做接口测试,频繁遇到了文件上传的接口,与其他一般接口的处理方式不一样,想着分享下,希望能给测试同学一点启发。 文章将围绕三个部分进行展开: 一、用户场景 二、接口请求参数 三、JMeter脚本编写步骤…...
自动化测试需知的4项测试工具!
一般来说学自动化会建议大家先学selenium,因为最早的时候,自动化就代表selenium,进入测试行业就开始做接口测试,而且现在基本每个公司都需要接口测试。今天就和大家聊一下接口测试的工具。 一、Robot Framework 机器人框架。之所…...
【深度学习】clip-interrogator clip docker 容器启动过程
文章目录 dockerfile备忘ENTRYPOINT ["bash", "/app/startProject.sh"]常用docker指令web服务脚本访问接口文件 给一张图片,输出图片描述。 dockerfile备忘 只有从dockerfile制作的镜像才有分层结构,加速传输,故第一步…...
Linux设备驱动之gpio-keys
Linux设备驱动之gpio-keys 前两个章节介绍了Linux字符设备和platform设备的注册,他们都是比较基础的,让大家理解Linux内核的设备驱动是如何注册、使用的。但在工作中,个人认为完全手写一个字符设备驱动的机会比较少,更多的都是基…...
【vue3页面展示代码】展示代码codemirror插件
技术版本: vue 3.2.40、codemirror 6.0.1、less 4.1.3、vue-codemirror 6.1.1、 codemirror/lang-vue 0.1.2、codemirror/theme-one-dark 6.1.2 效果图: 1.安装插件 yarn add codemirror vue-codemirror codemirror/lang-vue codemirror/theme-one-dar…...
【面试必刷TOP101】链表相加 单链表的排序
目录 题目:链表相加(二)_牛客题霸_牛客网 (nowcoder.com) 题目的接口: 解题思路: 代码: 过啦!!! 题目:单链表的排序_牛客题霸_牛客网 (nowcoder.com) 题目的接口:…...
Visual Studio复制、拷贝C++项目与第三方库配置信息到新的项目中
本文介绍在Visual Studio软件中,复制一个已有的、配置过多种第三方库的C项目,将其拷贝为一个新的项目,同时使得新项目可以直接使用原有项目中配置好的各类**C**配置、第三方库等的方法。 在撰写C 代码时,如果需要用到他人撰写的第…...
rust迭代器
迭代器用来遍历容器。 迭代器就是把容器中的所有元素按照顺序一个接一个的传递给处理逻辑。 Rust中的迭代器 标准库中定义了Iterator特性 trait Iterator {type Item;fn next(&mut self) -> Option<Self::Item>; }实现了Iterator特性的结构体就是迭代器。 很多类…...
软件定制开发的优势与步骤|APP搭建|小程序
软件定制开发的优势与步骤|APP搭建|小程序 定制开发的优势: 1. 满足特定需求:定制开发可以根据客户的实际需求进行设计和开发,使得软件系统能够更好地满足客户的业务目标。 2. 优化用户体验:通过深入了解客户的需求,定…...
ERR_CONNECTION_REFUSED等非标准的HTTP错误状态码原因分析和解决办法
文章目录 一、DNS Resolution Failed1,DNS服务器故障2,DNS配置错误3,DNS劫持4,域名过期-5,其他网络问题 二、ERR_CONNECTION_REFUSED-"ERR_CONNECTION_REFUSED" 错误可能有多种原因 三、ERR_SSL_PROTOCOL_ER…...
瀑布流 - Vue3基于Grid布局简单实现一个瀑布流组件
瀑布流 - Vue3基于Grid布局简单实现一个瀑布流组件 前言 在学习Grid布局之时,我发现其是CSS中的一种强大的布局方案,它将网页划分成一个个网格,可以任意组合不同的网格,做出各种各样的布局,在刷某书和某宝首页时&…...
ES6面试题总结
1. 谈谈你对 ES6 的理解,为什么要学习es6? ES6是新一代的JS语言标准,对分JS语言核心内容做了升级优化,规范了JS使用标准,新增了JS原生方法,使得JS使用更加规范,更加优雅,更适合大型应用的开发。学习ES6是成…...
mybatisplus,jdbc 批量插入
1.测试用例 项目中遇到在做导入号码的时候我们会用到批量导入,提高入库的速度。接下来我们以10000条为测试用例。 1.1 批量执行sql语句 当需要成批插入或者更新记录时,可以采用Java的批量更新机制,这一机制允许多条语句一次性提交给数据库…...
如何使用IP归属地查询API来追踪网络活动
引言 在当今数字化世界中,了解网络活动的源头和位置对于网络安全、市场研究和用户体验至关重要。IP归属地查询API是一种强大的工具,可以帮助您追踪网络活动并获取有关IP地址的重要信息。本文将探讨如何使用IP归属地查询API来追踪网络活动,以…...
【SQL】S0 系列博文大纲
系列博文大纲 SQL 学习环境建议系列博文相关书籍系列博文大纲阶段进展 SQL 学习环境建议 对于 SQL 语言的学习,博主本地使用:MySQL DataGrip; MySQL 提供本地数据库服务; DataGrip IDE,承担编程运行测试任务…...
2023年8月体育用品行业数据分析(京东数据产品)
当前,亚运会临近,这也带动了国民对体育消费的热情,体育产品内销逐渐旺盛,“亚运经济”红利开始显现。鲸参谋数据显示,今年8月份,京东平台上体育用品行业的销量为185万,同比增长2%;销…...
国内高校镜像网站
国内各大高校开源镜像站 排名不分前后 清华大学:https://mirrors.tuna.tsinghua.edu.cn/ 北京大学:https://mirrors.pku.edu.cn/ 北京外国语大学:http:// https://mirrors.bfsu.edu.cn/ 北京理工大学:https://mirrors.bit.e…...
Linux安装kafka-manager
相关链接https://github.com/yahoo/kafka-manager/releases kafka-manager-2.0.0.2下载地址 百度云链接:https://pan.baidu.com/s/1XinGcwpXU9YBF46qkrKS_A 提取码:tzvg 一、安装部署 1.把kafka-manager-2.0.0.2.zip拷贝到目录 /opt/app/elk 2.解压…...
MYSQL索引——B+树讲解
B-/B树看 MySQL索引结构 B-树 B-树,这里的 B 表示 balance( 平衡的意思),B-树是一种多路自平衡的搜索树.它类似普通的平衡二叉树,不同的一点是B-树允许每个节点有更多的子节点。下图是 B-树的简化图. B-树有如下特点: 所有键值分布在整颗树中; 任何一…...
别再折腾了!STM32CubeMX+Keil 5+Proteus 8.9保姆级联调配置,一次搞定
STM32开发环境联调实战:从零搭建CubeMXKeilProteus高效工作流 第一次接触STM32开发时,我被各种工具链的配置折磨得焦头烂额——CubeMX生成的工程在Keil里报错、Proteus仿真时芯片毫无反应、Debug选项神秘消失...如果你也经历过这种绝望,这篇文…...
2026年全国优质化妆培训机构深度盘点
颜值经济持续升温背景下,2026年化妆行业迎来规范化、专业化发展新阶段,涵盖影视、时尚、婚庆等多个领域,市场对专业化妆人才的需求持续攀升。据相关行业数据显示,2026年美业职业培训市场预计保持11%以上增速,美业门店专…...
现代Web应用特性管理:从概念到工程实践
1. 项目概述:一个面向现代Web开发的特性管理工具 如果你和我一样,长期在Web应用开发的一线摸爬滚打,那你一定对“特性开关”这个概念不陌生。简单来说,它就像你家里电灯的总闸,可以随时控制某个功能是“亮”还是“灭”…...
别再复制粘贴了!用LabVIEW 2023实现TCP/IP通讯的保姆级教程(附完整DEMO下载)
LabVIEW 2023 TCP/IP通讯实战:从原理到健壮性设计的深度解析 在工业自动化与测试测量领域,稳定可靠的通讯系统如同设备的神经系统。许多LabVIEW开发者虽然能够通过复制粘贴完成基础通讯功能,却在真实项目中频繁遭遇数据丢失、连接不稳定等&qu…...
别再死记公式了!用Python的NumPy和Matplotlib玩转坐标转换(附象限处理代码)
用Python实战坐标转换:从数学公式到可视化应用 坐标转换是计算机图形学、机器人学和数据可视化中的基础操作。传统教学中,我们往往被要求死记硬背转换公式,却很少有机会直观理解其实际应用场景。本文将带你用NumPy和Matplotlib这两个Python利…...
边缘计算大模型部署实战:从LLaMA量化到树莓派推理优化
1. 项目概述:一个为边缘计算优化的轻量级大语言模型最近在折腾边缘设备上的AI应用,发现一个挺有意思的项目——KuiperLLama。这名字听起来就很有“边缘”感,Kuiper(柯伊伯带)是太阳系边缘的一个区域,用它来…...
Armv8-A原子操作指令解析与应用优化
1. A64原子操作指令概述在Armv8-A架构中,A64指令集提供了一组强大的原子操作指令,这些指令在多核处理器环境下对实现线程安全的并发操作至关重要。原子操作的核心特性是保证特定内存操作的不可分割性——即这些操作要么完全执行,要么完全不执…...
终极指南:如何用sndcpy将Android音频无损转发到电脑
终极指南:如何用sndcpy将Android音频无损转发到电脑 【免费下载链接】sndcpy Android audio forwarding PoC (scrcpy, but for audio) 项目地址: https://gitcode.com/gh_mirrors/sn/sndcpy 你是否曾经想在电脑上收听手机上的音乐、播客或游戏音频࿱…...
基于Kubernetes Operator的企业级区块链网络自动化部署实践
1. 项目概述:企业级区块链的云原生部署方案如果你正在寻找一个能够将企业级区块链网络快速、稳定地部署到Kubernetes集群上的成熟方案,那么ConsenSys开源的quorum-kubernetes项目绝对值得你花时间深入研究。这个项目不是一个简单的概念验证,而…...
从零构建高频无线传输系统:调幅技术实战解析
1. 调幅无线传输系统入门指南 第一次接触调幅无线传输系统时,我也被各种专业术语搞得一头雾水。简单来说,调幅(AM)就是通过改变载波信号的幅度来传递信息的技术。想象一下快递员送包裹:载波就像快递车,而我们要发送的信息就是包裹…...
