⛄ 内容介绍
⛄ 部分代码
clear all;close all;clc;
xx=xlsread('负荷数据.xlsx');
M = 500;%观测值个数
N = 256;%信号x的长度
% K = 10;%信号x的稀疏度
a1=0.2;a2=0.4;a3=0.6;a4=0.8;
s1=30;
x = xx(1:N,1);
Psi = eye(N);%x本身是稀疏的,定义稀疏矩阵为单位阵x=Psi*theta
Phi1 = randn(M,N)/sqrt(M);%测量矩阵为高斯矩阵
A1 = Phi1 * Psi;%传感矩阵
y1 = Phi1 * x;%得到观测向量y
Phi2 = randn(M,N)/sqrt(M);%测量矩阵为高斯矩阵
A2 = Phi2 * Psi;%传感矩阵
y2 = Phi2 * x;%得到观测向量y
Phi3 = randn(M,N)/sqrt(M);%测量矩阵为高斯矩阵
A3 = Phi3 * Psi;%传感矩阵
y3 = Phi3 * x;%得到观测向量y
Phi4 = randn(M,N)/sqrt(M);%测量矩阵为高斯矩阵
A4 = Phi4 * Psi;%传感矩阵
y4 = Phi4 * x;%得到观测向量y
%% 恢复重构信号x
tic
theta1 = CS_SWOMP( y1,A1,s1,a1);
x

最低0.47元/天 解锁文章

被折叠的 条评论
为什么被折叠?



