导入数据

加载wecg信号进入MATLAB工作空间。采样频率为180hz。

负载wecgFs = 180;

可视化量图

开放小波时频分析仪．在分析仪选项卡上,单击进口的信号．出现一个窗口，其中列出了应用程序可以处理的所有工作区变量。选择wecg信号并点击进口．在简单的一次性初始化之后，导入信号的名称和类型将填充信号窗格。应用程序显示的标量图Scalogram-wecg情节。该应用程序生成的尺度图使用分析莫尔斯(3,60)小波和类使用默认设置函数。影响锥显示尺度图中可能受影响的区域边缘效应的构件。由于存在潜在的边缘效应，您应该将虚线白线以外的阴影区域视为可疑区域。要了解更多关于影响锥的信息，请参见边界效应和影响锥．

属性中的选项，可以隐藏或显示影响边界线的锥体，或为边界区域加阴影首选项▼菜单。如果信号是复值，还可以选择将正(逆时针)和负(顺时针)分量显示为单独的或串联的标量图。中所选择的选项首选项▼菜单在MATLAB会话中持久化。

默认情况下，由于信号不是时间表，应用程序以周期/样本为单位绘制频率，并使用样本索引作为时间轴的基础。如果导入的信号是一个时间表，应用程序将标量图绘制为频率的函数(以赫兹为单位)，并使用时间表的行时间作为时间轴的基础。

在分析仪选项卡，选择采样率单选按钮。尺度图中的轴标签使用默认采样率1hz进行更新。为对象设置适当的轴标签wecg信号，应用值180在采样率字段。不能修改时间表的采样率。

修改小波参数

若要访问用于创建标量图的CWT参数，请单击量图选项卡。类的输入参数对应类函数。参数设置为默认值。因为您设置了采样率，所以最小和最大频率限制以赫兹为单位。默认频率限制取决于小波、信号长度、采样率和每八度的声音。有关更多信息，请参见cwtfreqbounds．

为了获得莫尔斯(40,60)小波的尺度图，首先，改变的莫尔斯对称参数对称文本框改为40。然后，要输入新值，请按回车键在键盘上或在文本框外的任何地方单击鼠标。

单击，可将CWT参数恢复为默认值重新设置参数．重新设置参数将启用计算量图按钮。

注意:为了防止无效的设置，应用程序会立即验证您更改的任何参数。如果输入无效值，应用程序会自动将其替换为有效值。应用程序将以下值视为无效值:

新值可能不是所需的值。为了避免意外的结果，您应该确保您输入的任何值总是得到有效的设置。有关更多信息，请参见示例调整莫尔斯小波参数．

要应用更改并可视化新的标量图，请单击计算量图．

量图进行比较

要将当前的标量图与使用凹凸小波获得的标量图进行比较，首先单击重复的按钮。分析仪选项卡。第二个信号，wecgCopy，出现在信号窗格。副本的标量图出现在Scalogram-wecgCopy情节。然后在量图选项卡上,选择撞从小波下拉菜单。观察莫尔斯小波参数现在已禁用，频率限制已更新。要使用凹凸小波创建标量图，请单击计算量图．要与第一个标量图进行比较，请选择wecg在信号窗格。

导出的结果

您可以将结构数组导出到包含所选信号的CWT的工作空间。您还可以生成MATLAB®脚本来在工作空间中重新创建标量图。

出口结构

输出的CWTwecgCopy中选择该信号信号窗格。然后选择导出量图从出口▼菜单创建结构数组wecgCopy_scalogram在你的工作空间。结构数组有三个字段:

通过执行以下命令，您可以使用字段值在新图形中可视化标量图:

图pcolor (wecgCopy_scalogram.timeVector,.．.wecgCopy_scalogram.frequencyVector,.．.abs (wecgCopy_scalogram.coefficients))阴影平集(gca),“yscale”，“日志”)标题(“量图”)包含(“时间(s)”) ylabel (“频率(赫兹)”）

生成脚本

属性的标量图生成一个脚本wecg信号，选择wecg信号在信号窗格。然后选择生成MATLAB脚本从出口▼菜单。程序中选择的信号名称会生成脚本信号窗格。在MATLAB编辑器中打开一个未命名的脚本，其中包含以下可执行代码。

%的参数sampleRate = 180;小波参数= [40,60];计算时间向量t = 0:1/sampleRate:(长度(wecg)*1/sampleRate)-1/sampleRate;%计算类如果有必要，在下面的代码中用工作区变量名代替wecg作为cwt()函数的第一个输入运行下面不带输出参数的函数调用来绘制结果[小波变换，频率]= cwt(wecg, sampleRate，.．.小波参数=小波参数);尺度图= abs(小波变换);

保存并执行脚本。您可以使用工作区变量量图，t,频率使标量图形象化。

Subplot (2,1,1) pcolor(t，频率，标量图)阴影平集(gca),“yscale”，“日志”)标题(“量图”) ylabel (“频率(赫兹)”)包含(“时间(s)”) subplot(2,1,2) plot(t,wecg)轴紧标题(“信号”)包含(“时间(s)”）