自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(44)
  • 收藏
  • 关注

原创 matlab散点图

disp(['归一化后范围: [', num2str(min(density_normalized)), ', ', num2str(max(density_normalized)), ']'])disp(['原始密度范围: [', num2str(min(density)), ', ', num2str(max(density)), ']'])scatter(X, Y, 25, density_normalized, 'filled');title('归一化密度编码散点图');% 验证归一化结果(可选)

2025-03-08 09:22:59 265

原创 matlab相位图

角位移theta(t) 振荡衰减。% 位移y(t) 振荡衰减。title('相位图:角位移 \theta 与线性位移 y 的关系');% 模拟生成时间t,位移y(t)和角位移theta(t)title('a) 线性位移随时间');title('b) 角位移随时间');

2024-10-18 12:27:09 609

原创 free vibration

EX9.11简支梁突然卸载的瞬态动力分析(完全法)OUTRES,ALL,ALL!ALPHAD,0!注释后不考虑阻尼则为自由振动。LSWRITE,2!副除荷载,写人第2荷载步。KBC,1!LSWRITE,1!BETAD,0!PLVAR,2!创建模型(解释从略)

2024-09-30 16:34:28 244

原创 求解方程111

竖弯频率 其实是角速度。% 扭转频率 其实是角速度。I = 0.666;m = 14.92;rho = 1.25;

2024-09-24 10:14:08 359

原创 求解未知方程

使用 ode45 求解。% 初始化存储结果的数组。% 循环遍历每个 a 值。% 定义 a 的范围。

2024-08-14 09:09:41 328 1

原创 输出数据到excel中

对于MATLAB R2019a及以后版本。% 或使用 xlswrite(data, filename);对于较老版本的MATLAB。% 初始化一个单元数组来存储所有数据。% 使用循环将数据填充到单元数组中。% 将数据写入Excel文件。% 定义Excel文件名。

2024-06-14 11:03:11 1962

原创 信号相关性分析

添加噪声使其非线性。disp(['斯皮尔曼相关系数: ', num2str(rank_corr)]);% 添加噪声使其非线性。% 计算斯皮尔曼相关系数。

2024-06-14 09:40:29 322

原创 ansys有限元分析

在x=10处施加一个向下的力。选择 x=10 的节点。f,all,fy,0!deltim,0.1!set,last!nlgeom,on!瞬态加载步的时间设置。

2024-06-13 18:43:39 594

原创 matlab BP神经网络

100组输入,每组10个特征。% 100组输出,每组1个输出值。% 将数据分成训练集和测试集。% 设置训练、验证和测试数据。% 计算均方误差(MSE)% 创建前馈神经网络。

2024-06-08 12:36:30 479

原创 相关性分析

使用正确的采样频率。disp(['两个信号之间的相关系数为:', num2str(correlation_coefficient)]);ylabel('功率谱密度 (dB)');xlabel('频率 (Hz)');title('交叉功率谱密度图');% 绘制交叉功率谱密度图。% 计算交叉功率谱密度。

2024-05-07 12:42:44 572

原创 滤波信号总结

这里使用sin函数作为加速度信号的示例,你可以根据需要修改这部分代码。% 绘制加速度信号、带噪声的加速度信号和滤波后的加速度信号。% 加速度信号的频率。% 加速度信号的幅值。% 添加白噪声到加速度信号。% 生成加速度时程信号。

2024-04-07 23:21:37 475

原创 添加背景图片画图matlab

set(gca,'color','none') %这里以前自己没有注意。

2024-04-05 23:34:04 529

原创 极坐标散点图

将数据范围调整为1-2。% 设置颜色条的数据范围。% 假设 arry1 是你的变量。% 将极坐标下的角度转换为弧度。% 添加颜色条,并设置数据范围。% 定义颜色映射范围。

2024-04-02 13:56:19 372

原创 根据数据的密度进行作图

生成颜色矢量,根据点的密度给定颜色深浅。title('密度相关的散点图');xlabel('X轴');ylabel('Y轴');% 定义网格并计算密度。

2024-04-01 23:28:39 410

原创 渐变颜色作图

计算对应的 y 值(二次函数分布)% 绘制散点图,并根据密度设置颜色。% 添加一些随机噪声。

2024-04-01 21:26:14 236

原创 尝试一下线性拟合

clearclcclose allfigure(1)hold onhold offgrid on。

2024-03-31 23:47:19 328

原创 找到零值点

for i = 2:length(signal)-1 % 从第2个点开始,到倒数第2个点结束。% 存储最靠近0的点的索引。% 存储开始点的索引。% 存储结束点的索引。% 检查零点右侧的点是正还是负,以确定是开始点还是结束点。% 检查当前点及其相邻点,确保中间点是极小值。

2024-03-10 15:35:07 435

原创 滑动窗算一下rms

这里的10是平滑窗口的长度。overlap=5;%每次往前推进的数值。

2024-03-10 14:39:56 478

原创 绘制一下包络线

衰减曲线带有随机信号。% 时间向量,总时长为5秒。legend('衰减信号', '包络线');title('衰减信号及其包络线');xlabel('时间 (s)');% 生成衰减曲线带有随机信号。ylabel('幅度');% 绘制信号及其包络线。

2024-03-03 10:21:33 618

原创 处理一下异常值

5分钟,单位:秒。% 采样频率,单位:Hz。% 将不满足三西格玛准则的值替换为平均值。xlabel('时间 (秒)');% 计算当前分钟内的索引范围。% 获取当前分钟内的风速数据。title('风速时程');% 计算三西格玛准则的阈值。ylabel('风速');% 对风速数据进行处理。% 计算平均值和标准差。% 生成随机风速数据。% 更新原始风速数据。

2024-02-13 23:50:09 249

原创 处理一下异常值

用中位数替换异常值。disp('异常值索引:');

2024-02-03 19:56:12 233

原创 时间序列问题

创建时间点序列,每隔一分钟。% 创建起始时间和结束时间。% 找到特定时间点的行。

2024-01-29 17:59:59 165

原创 风速编辑一下

plot(time, pulsating_wind_speed_vertical, 'LineWidth', 1.5, 'DisplayName', '垂直脉动风速');plot(time, pulsating_wind_speed_parallel, 'LineWidth', 1.5, 'DisplayName', '平行脉动风速');plot(time, total_wind_speed, 'LineWidth', 1.5, 'DisplayName', '合风速');title('合风速时程');

2024-01-29 00:42:38 423

原创 短时傅里叶变换计算一下

clearclcclose allfs=100;figure(1)figure(2)figure(3)view(2)figure(5)view(2)

2024-01-23 21:01:40 443

原创 标记一下点

找到功率谱密度图中的峰值点。

2024-01-21 23:13:45 375

原创 MATLAB标记点

定义阈值,即前后20个点的纵坐标差值的阈值。% 在每个峰值点的位置添加文本标签。% % 生成随机时程信号。% 生成随机时程信号。

2024-01-21 18:58:22 1041

原创 功率谱密度

时间从0到300秒,步长为采样间隔。ylabel('功率谱密度 (dB/Hz)');% 使用pwelch计算功率谱密度函数。xlabel('频率 (Hz)');xlabel('时间 (秒)');title('功率谱密度函数');title('间谐时程信号');ylabel('幅值');% 绘制功率谱密度函数。

2024-01-21 15:13:59 492 1

原创 处理一下时间

clearclcclose all% 输入日期和时刻% 创建datetime对象% 显示结果。

2024-01-21 13:00:20 425 1

原创 指数函数拟合

衰减曲线带有随机信号。% 时间向量,总时长为5秒。% 生成衰减曲线带有随机信号。% 对峰值进行指数函数拟合。% 绘制原始曲线和拟合曲线。% 寻找衰减曲线的峰值。

2024-01-18 12:04:38 989 1

原创 自定义拟合

plot(number, fitCurve, 'r', 'DisplayName', '拟合曲线');scatter(x, y, 'b', 'DisplayName', '原始数据');title('自定义函数拟合(nlinfit)');legend('衰减曲线', '峰值点');title('衰减曲线及其峰值点');xlabel('时间(秒)');% 用拟合的参数生成拟合曲线。ylabel('振幅');

2024-01-18 10:17:22 373 1

原创 自定义函数拟合

plot(x, fitCurve, 'r', 'DisplayName', '拟合曲线');scatter(x, y, 'b', 'DisplayName', '原始数据');title('自定义函数拟合(nlinfit)');% 用拟合的参数生成拟合曲线。% 绘制原始数据和拟合曲线。% 定义自定义拟合函数。

2024-01-18 09:31:05 482 1

原创 寻找峰值点

衰减曲线带有随机信号。% 时间向量,总时长为5秒。legend('衰减曲线', '峰值点');title('衰减曲线及其峰值点');xlabel('时间(秒)');% 生成衰减曲线带有随机信号。ylabel('振幅');% 绘制原始曲线和峰值点。% 寻找衰减曲线的峰值。

2024-01-17 21:53:23 367

原创 衰减时程曲线

ylabel('功率谱密度 (dB/Hz)');ylabel('功率谱密度 (dB/Hz)');title('原始曲线功率谱密度');xlabel('频率 (Hz)');xlabel('频率 (Hz)');xlabel('时间 (秒)');xlabel('时间 (秒)');ylabel('振幅');

2024-01-17 20:34:02 505

原创 matlab画一张风速时程

平均风速(米/秒)ylabel('功率谱密度 (dB/Hz)');ylabel('功率谱密度 (dB/Hz)');ylabel('风速(米/秒)');ylabel('风速(米/秒)');xlabel('频率(赫兹)');xlabel('频率(赫兹)');

2024-01-03 23:56:37 984 1

原创 画一张时程图

假设我们有四个数据向量:u, v, w 和 direction。% 这些向量的长度都是 600,代表 600 秒的数据。% 假设我们要计算每个连续的 5 个数的平均值。% 绘制第四个子图:direction。% 创建一个新的图形窗口。% 绘制第一个子图:u。% 绘制第二个子图:v。% 绘制第三个子图:w。% 设置 x 轴的标签。

2024-01-02 23:51:15 468

原创 给matlab图改变一下字体等

设置X轴标签字体大小和粗细。% 设置Y轴标签字体大小和粗细。

2023-12-26 23:39:11 1538 1

原创 matlab正态分布拟合数据画图

legend(['拟合曲线 (\mu = ' num2str(meanValue) ', \sigma = ' num2str(stdValue) ')']);% 生成1000个均值为5,标准差为2的随机数。% 20表示直方图的条数。% 使用正态分布拟合工具进行拟合。title('正态分布拟合');xlabel('数据值');ylabel('频率');% 绘制直方图和拟合曲线。

2023-12-26 00:13:59 2265 1

原创 matlab计算一下风向夹角

fprintf('合风速与北方向的夹角:%.2f度\n', theta_deg);fprintf('合风速:%.2f\n', V);% 计算合风速与北方向的夹角(以度为单位)% 南北方向向北为正,东西方向向西为正。

2023-12-24 22:15:51 663 1

原创 matlab画一张极坐标图

随机生成0到360度之间的风向。% 随机生成1到10的风速。% 绘制风速和风向的极坐标图。% 按顺时针方式递增。% 将风速和风向转换为极坐标坐标系的坐标。polaraxes;% 生成随机风向和风速数据。

2023-12-23 00:20:37 2390 1

原创 对我的文件时间名称数据排个序

folderPath = 'C:\Users\wind\Desktop\新建文件夹';cd('C:\Users\Desktop\新建文件夹')% 提取文件的日期信息并转换为datetime对象。disp('按年月日排序后的文件名:');% 获取文件夹中的文件列表。% 显示排序后的文件名。% 按日期排序文件名。

2023-12-21 00:18:15 395

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除