matlab实验数据处理
2024-05-19 08:53:33
%%
%------------------------------------------------------------------------三组信号一致性处理
clear all;
close all;
%22,0,[22,0,7000,2]\\\22
data1=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
data2=csvread('E:\data\600-600\9.24\0\2-5-0.csv',22);
data3=csvread('E:\data\600-600\9.24\0\3-6-0.csv',22);
t=data1(:,1);
figure(2);%t,data1(:,2),t,data2(:,2),t,data3(:,3)
plot(t,data1(:,2),t,data2(:,2),t,data3(:,2));
title('30khz,1-4,2-5,3-6');
%%
%--------------------------------------------------------------------一组信号差分处理
clear all;
close all;
%22,0,[22,0,7000,2]\\\22
data1=csvread('E:\data\600-600\9.24\0\1-5-0.csv',22);
data2=csvread('E:\data\600-600\9.24\1\1-5-1.csv',22);
t=data1(:,1);figure(1);
subplot(3,1,1);
plot(t,data1(:,2),t,data2(:,2));
title('1-5');
subplot(3,1,2);
plot(data1(:,1),data1(:,2)-data2(:,2));
subplot(3,1,3);
plot(t,data1(:,3));
%-----------------------------------------------------------------2000-8000
t1=data1(2000:6000,1);
figure(3);
subplot(2,1,1);
plot(t1,data1(2000:6000,2),t1,data2(2000:6000,2));
title('1-5');
subplot(2,1,2);
plot(data1(2000:6000,1),data1(2000:6000,2)-data2(2000:6000,2));
%%
%----------------------------------------------------------------------多组裂纹差分信号处理
clear all;
close all;data1_1_1=csvread('E:\data\600-600\9.24\0\6-1-0.csv',22);
data1_2_1=csvread('E:\data\600-600\9.24\0\6-2-0.csv',22);
data1_3_1=csvread('E:\data\600-600\9.24\0\6-3-0.csv',22);
data1_4_1=csvread('E:\data\600-600\9.24\0\6-4-0.csv',22);
data1_5_1=csvread('E:\data\600-600\9.24\0\6-5-0.csv',22);
% data1_6_1=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
% data1_7_1=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
% data1_8_1=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
data1_1_2=csvread('E:\data\600-600\9.24\1\6-1-1.csv',22);
data1_2_2=csvread('E:\data\600-600\9.24\1\6-2-1.csv',22);
data1_3_2=csvread('E:\data\600-600\9.24\1\6-3-1.csv',22);
data1_4_2=csvread('E:\data\600-600\9.24\1\6-4-1.csv',22);
data1_5_2=csvread('E:\data\600-600\9.24\1\6-5-1.csv',22);
% data1_6_2=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
% data1_7_2=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
% data1_8_2=csvread('E:\data\600-600\9.24\0\1-4-0.csv',22);
% D=[data1_1_1(:,1),data1_1_1(:,2)+data1_1_2(:,2),data1_2_1(:,2)+data1_2_2(:,2),data1_3_1(:,2)+data1_3_2(:,2),data1_4_1(:,2)+data1_4_2(:,2),data1_5_1(:,2)+data1_5_2(:,2),data1_6_1(:,2)+data1_6_2(:,2),data1_7_1(:,2)+data1_7_2(:,2),data1_8_1(:,2)+data1_8_2(:,2)];A=[data1_1_1(:,1),data1_1_1(:,2)-data1_1_2(:,2),data1_2_1(:,2)-data1_2_2(:,2),data1_3_1(:,2)-data1_3_2(:,2),data1_4_1(:,2)-data1_4_2(:,2),data1_5_1(:,2)-data1_5_2(:,2)];figure(4);
subplot(3,1,1);
plot(A(:,1),A(:,4),A(:,1),A(:,5));
title('3-3,3-4');
subplot(3,1,2);
plot(A(:,1),A(:,2),A(:,1),A(:,3));
title('3-1,3-2');
subplot(3,1,3);
plot(A(:,1),A(:,6));
title('3-5');%%
%-----------------------------------信号一致性
clear all;
close all;
data1=csvread('E:\data\600-600\10.12\0\1-4.csv',22);
data2=csvread('E:\data\600-600\10.12\0\2-5.csv',22);
data3=csvread('E:\data\600-600\10.12\0\3-6.csv',22);
t=data1(:,1);
figure(1);
plot(t,data1(:,3),t,data2(:,3),t,data3(:,3));%%
%----------------------------------单组信号、差分信号及激励信号
clear all;
close all;
data1_0=csvread('H:\11.17\chirp\1-5.csv',22);
data2_1=csvread('E:\data\600-600\11.17\35khz\1\1-5.csv',22);h = @(x) abs(hilbert(x));
t=data1_0(:,1);
figure(1);
subplot(3,1,1);
plot(t,data1_0(:,3),t,data2_1(:,3));
title('6-4');
subplot(3,1,2);
plot(data1_0(:,3)-data2_1(:,3));
subplot(3,1,3);
%plot(h(data1_0(:,2)));
plot(h(data1_0(:,2)));
%%
%-------------------------------------多组差分信号
clear all;
close all;
data1_1_1=csvread('E:\data\600-600\11.17\35khz\0\1-2.csv',22);
data1_2_1=csvread('E:\data\600-600\11.17\35khz\0\1-3.csv',22);
data1_3_1=csvread('E:\data\600-600\11.17\35khz\0\1-4.csv',22);
data1_4_1=csvread('E:\data\600-600\11.17\35khz\0\1-5.csv',22);
data1_5_1=csvread('E:\data\600-600\11.17\35khz\0\1-6.csv',22);
% data1_6_1=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);
% data1_7_1=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);
% data1_8_1=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);
data1_1_2=csvread('E:\data\600-600\11.17\35khz\1\1-2.csv',22);
data1_2_2=csvread('E:\data\600-600\11.17\35khz\1\1-3.csv',22);
data1_3_2=csvread('E:\data\600-600\11.17\35khz\1\1-4.csv',22);
data1_4_2=csvread('E:\data\600-600\11.17\35khz\1\1-5.csv',22);
data1_5_2=csvread('E:\data\600-600\11.17\35khz\1\1-6.csv',22);
% data1_6_2=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);
% data1_7_2=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);
% data1_8_2=csvread('E:\data\600-600\9.29\0\6-4-0.csv',22);A=[data1_1_1(:,1),data1_1_1(:,3)-data1_1_2(:,3),data1_2_1(:,3)-data1_2_2(:,3),data1_3_1(:,3)-data1_3_2(:,3),data1_4_1(:,3)-data1_4_2(:,3),data1_5_1(:,3)-data1_5_2(:,3)];
figure(4);
subplot(3,1,1);
plot(A(:,1),A(:,2),A(:,1),A(:,3));
title('X-1,X-2');
subplot(3,1,2);
plot(A(:,1),A(:,4),A(:,1),A(:,5));
title('X-3,X-4');
subplot(3,1,3);
plot(A(:,1),A(:,6));
title('X-6');
%%clear all;
close all;
data15=xlsread('E:\data\600-600\dp-10000\15khz\40-1109-15khz.xlsx',6,'A2:G9000');
data20=xlsread('E:\data\600-600\dp-10000\20khz\40-1030-20khz.xlsx',6,'A2:G9000');
data25=xlsread('E:\data\600-600\dp-10000\25khz\40-1030-25khz.xlsx',6,'A2:G9000');
data30=xlsread('E:\data\600-600\dp-10000\30khz\40-1109-30khz.xlsx',6,'A2:G9000');
data35=xlsread('E:\data\600-600\dp-10000\35khz\40-1109-35khz.xlsx',6,'A2:G9000');
data40=xlsread('E:\data\600-600\dp-10000\40khz\40-1030-40khz.xlsx',6,'A2:G9000');
hdata15=zeros(length(data15(:,1)),7);
hdata20=zeros(length(data15(:,1)),7);
hdata25=zeros(length(data15(:,1)),7);
hdata30=zeros(length(data15(:,1)),7);
hdata35=zeros(length(data15(:,1)),7);
hdata40=zeros(length(data15(:,1)),7);
A=zeros(7000,7);h = @(x) abs(hilbert(x));for i=1:7
hdata15(:,i)=h(data15(:,i));
hdata20(:,i)=h(data20(:,i));
hdata25(:,i)=h(data25(:,i));
hdata30(:,i)=h(data30(:,i));
hdata35(:,i)=h(data35(:,i));
hdata40(:,i)=h(data40(:,i));
end for i=1:7000A(i,:) =hdata15(i+1100,:)+hdata20(i+1650,:)+ hdata25(i+1350,:)+hdata30(i+1100,:)+hdata35(i+900,:)+hdata40(i+920,:);%A(i,:) = hdata15(i+1100,:)+hdata20(i+1650,:)+hdata25(i+1350,:)+hdata30(i+1100,:)+hdata35(i+900,:)+hdata40(i+920,:);% A(i,:) = data1(i+1,:)+data2(i+2,:)+data3(i+3,:);
% A(i,:) = hdata1(i+150,:).*hdata2(i+110,:);%A(i,:) = hdata1(i+150,:).*hdata2(i+110,:).*hdata3(i+110,:);A(i,1)=data15(i+1100,1);end·
matlab实验数据处理相关推荐
- matlab实验数据处理,实验一 数据处理方法的MATLAB实现
实验一数据处理方法的MATLAB实现 一.实验目的 学会在MATLAB环境下对已知的数据进行处理. 二.实验方法 1. 求取数据的最大值或最小值. 2. 求取向量的均值.标准方差和中间值. 3．在MA ...
- MATLAB实验8：数据处理与多项式计算
一.实验目的 1.掌握数据统计与分析的方法 2.掌握数值插值与曲线拟合的方法及其应用 3.掌握多项式的常用计算 二.实验内容 1.利用matlab提供的rand函数生成30000个符合均匀分布的随机数 ...
- 如何用matlab画烧杯,MATLAB在动力学实验数据处理中的应用
· 68· 化 学 教 育(http://www．hxjy．org) 2015年第 2期 MATLAB在动力学实验数 据处理 中的应用 高锦红 (渭南 师范学 院化学与生命科学学院 陕西渭南 7140 ...
- matlab物理实验数据处理,利用Matlab对物理实验数据进行处理
利用Matlab对物理实验数据进行处理 李凌云;王海军;韩之乃 [期刊名称]<南京晓庄学院学报> [年(卷),期]2005(021)005 [摘要]本文以热电偶的定标实验.二极管伏安特性实 ...
- matlab软件在环测试步骤,Matlab软件在牛顿环实验数据处理方法上的创新
<Matlab软件在牛顿环实验数据处理方法上的创新>由会员分享,可在线阅读,更多相关<Matlab软件在牛顿环实验数据处理方法上的创新(5页珍藏版)>请在人人文库网上搜索. 1 ...
- matlab直观数据处理,霍尔效应实验的MATLAB数据处理
[专题研讨 ] 霍尔效应实验 的 数据处理 王雪晴 , 杜豪杰 评顶 山学院 电气信息工程学院 , 河南 平顶 山 摘要 运用 语言对大学物理 实验霍尔效应进行最小 二乘法数据处理 . 与传统手动作 ...
- singft 函数 matlab,（新）信号与系统MATLAB实验全.doc
实验篇 信号与系统实验指导 实验一.MATLAB编程基础及典型实例 一.实验目的 熟悉MATLAB软件平台的使用: 熟悉MATLAB编程方法及常用语句: 掌握MATLAB的可视化绘图技术: 结合< ...
- matlab实验是啥,matlab实验心得体会
matlab实验心得体会这个就是我自己做出来的,发现用soundsc听出来还行,但是用wavwrite就变的很难听. 后来发现PB写的很好,而且还能够把前面一段没有噪音的部分给保存下来,我就看了他的代 ...
- 控制工程matlab实验报告小结,控制工程MATLAB实验报告.doc
控制工程MATLAB实验报告 浙江科技学院 机电系统开放性实验 设计题目: MATLAB控制机电系统实分析 专 业: 机 械 设 计 制 造及 其 自 动 化 班 级: 机 制 122 学生姓名: 王 ...
最新文章
- 北航计算机考博C语言真题_北京航空航天大学计算机考研：严厉的老父亲，教会你本事，本科非211以上就不要去了...
- android使用软引用构建缓存
- 使用Docker实现丝般顺滑的持续集成
- 【Android FFMPEG 开发】OpenSLES 播放音频 ( 创建引擎 | 输出混音设置 | 配置输入输出 | 创建播放器 | 获取播放/队列接口 | 回调函数 | 开始播放 | 激活回调 )
- UA MATH565C 随机微分方程V Stationary Measure
- Caffe官方教程翻译（1）：LeNet MNIST Tutorial
- voc_eval.py:41: RuntimeWarning: invalid value encountered in greater_equal if np.sum(rec = t) ==
- synchronized实现可见性对比volatile
- gprs连接中断记录_请收好！脉搏波血压计连接天天血压APP使用指南
- 【C语言】数据类型转换的问题
- python 线性规划问题_一学高数，线代就头疼？让python帮你解决(内含教程)
- 作为window对象属性的元素 多窗口和窗体
- Windows server 2016 安装补丁报错 - The update is not applicable to your computer Error:0x800f0823
- go使用redis——redigo使用HMSET存储结构体
- Android library快速构建上传jcenter
- 初级工程师该如何去学习，如何去研发开关电源？
- DeepSort轨迹跟踪方法
- 1人30天44587行代码，分享舍得网开发经验【修订版】
- 研究团队再次发现安全漏洞，微软警告Azure用户
- 计算机专业在房管局,计算机故障 昨天房产局无法办过户
热门文章
- 2017word计算机操作题,2017职称计算机考试Word操作试题及答案
- 移植uClinux内核到stm32f407
- 使用esp8266作为I2C传感器的主机
- 过拟合欠拟合模拟 || 深度学习 || Pytorch || 动手学深度学习11 || 跟李沐学AI
- 富士施乐m115b怎么连接电脑_施乐 m115b打印机怎么安装
- JWT —— 生成Token、解析Token的简单工具类
- Java、JSP社区交流站点的分析与实现
- vue3+ts Type ‘xxx[]‘ is not assignable to type ‘never[]‘
- Kindle fo PC 提示已超过授权限制
- Eclipse快捷键大全(转载)