塑胶制品东莞网站建设,梅州做网站多少钱,企业集团网站建设,什么是所见即所得的网页制作工具使用python调用matlab API接口进行数据分析找到matlab安装目录下自带的setup.py的路径我的路径如下D:\MATLAB\R2018a\extern\engines\python打开cmd进入到刚才找到的路径管理员权限执行#xff01;#xff01;#xff01; pyhton setup.py install具体代码1.将csiTool中的ma…使用python调用matlab API接口进行数据分析找到matlab安装目录下自带的setup.py的路径我的路径如下D:\MATLAB\R2018a\extern\engines\python打开cmd进入到刚才找到的路径管理员权限执行 pyhton setup.py install具体代码1.将csiTool中的matlab文件夹下的文件拷到pycharm的工作目录下我的拷出来不行于是从Github中重新下载了下来考进去成功csi_traceread_bf_file(csi.dat);csi_entrycsi_trace{1}timestamp_low: 4 (In the sample trace, timestamp_low is invalid and always 4.)bfee_count: 72Nrx: 3Ntx: 1rssi_a: 33rssi_b: 37rssi_c: 41noise: -127agc: 38perm: [3 2 1]rate: 256csi: [1x3x30 double]timestamp_lowNIC网卡1MHz时钟的低32位。它大约4300s(72min)重复一回。 (从0-2^32需要4300s)时间戳相连两包此值差单位为微秒通过验证发现100hz的发包频率此差值为10000,20hz的发包频率此差值为50000此参数可以确定出波形的横轴时间。当发包频率过大会出现接收端停止收数的情况20Hz正常200Hz和1000Hz的采样频率收端最多接收一分多钟的的数据(原因是recv函数接收数据阻塞导致)1、将包长变小后还是有同样的问题。2、由数据速率引起的它太高程序无法处理(可以尝试设置非阻塞)0x4101决定了主机仅一根天线发送数据因此此处的Ntx值为1另外采用的是OFDM模式 如果想要收端接收到发端两根天线的数据(CSI的格式为2×3×30)只需要把第14、15位均设置为1即在发送端的脚本里面设置为0xC101 注关于0x4101的具体说明如下bfee_count驱动记录并发送到用户控件的波束测量值的总数。内核和用户空间中netlink频道是有损的可以用该变量来检测被丢弃的测量值。Nrx接收端使用的天线数量。Ntx发送端使用的天线数量。rssi_a, rssi_b, rssi_c由接收端NIC测量出的RSSI值。perm展示NIC如何将3个接收天线的信号排列到3个RF链上上图中的数据表示天线A被发送到RF链A天线B被发送到RF链B天线C被发送到RF链C。rate发包频率。csiCSI值3374094MATLAB命令a[1 2 3;4 5 6] #矩阵A #A的共轭转置矩阵ones(2,3) #一矩阵b1:0.1:2 #集合间距为0.1zeros(2,3) #零矩阵rand(2,3) #随机矩阵数值介于0和1之间randn(2,3) #矩阵(一个平均值为0的高斯分布方差或者等于1的标准偏差)eye(2) #单位阵hist() #直方图size() #矩阵的尺寸size(A, 1) #矩阵第一维的尺度length(A) #由于矩阵A是一个3×2的矩阵因此最大的维度应该是3因此该命令会返回3。who # 命令能显示出工作空间中的所有变量whos #命令能更详细地进行查看load(featureX.dat) # 加载数据文件变量名等于文件名clear #命令删除所有变量加变量名删除特定save hello.mat v #这个命令会将变量v存成一个叫 hello.mat 的文件按二进制存储save hello.txt v -ascii # 这样就会把数据存成一个文本文档或者将数据的 ascii 码存成文本文档键入 A(3,2) # 这将索引到 矩阵的 (3,2) 元素。键入A(2,:) #返回第二行的所有元素A([1 3],:) #取的是A矩阵的第一行和第三行的每一列A [A, [100, 101,102]] #这样做的结果是在原矩阵的右边附加了一个新的列矩阵A(:) #这是一个很特别的语法结构意思是把 A中的所有元素放入一个单独的列向量这样我们就得到了一个 9×1 的向量A.*B #这么做将矩阵 A中的每一个元素与矩阵 B中的对应元素相乘A.^2 #这将对矩阵中每一个元素平方1./A #得到中每一个元素的倒数。同样地这里的点号还是表示对每一个元素进行操作。exp(A) #自然数e的幂次运算abs(A) #取绝对值valmax(A) #默认情况下max(A)返回的是每一列的最大值如果你想要找出整个矩阵A的最大值你可以输入max(max(A))[val, ind] max(A) #这将返回矩阵中的最大值存入以及该值对应的索引A3 #这将进行逐元素的运算所以元素小于3的返回1否则返回0。find(A3) #这将告诉我A中的哪些元素是小于3的。A magic(3) #magic 函数将返回一个矩阵称为魔方阵或幻方 (magic squares)它们具有以下这样的数学性质它们所有的行和列和对角线加起来都等于相同的值。[r,c] find(A7) #这将找出所有矩阵中大于等于7的元素r 和c分别表示行和列sum(a) #就把 a 中每列所有元素加起来了。sum(a,2) #求每行的和prod(a) #prod 意思是product(乘积)它将返回这四个元素的乘积。floor(a) #向下四舍五入因此对于 a中的元素0.5将被下舍入变成0。ceil(a) #向上四舍五入所以0.5将上舍入变为最接近的整数也就是1。type(3) #这通常得到一个3×3的矩阵max(rand(3),rand(3)) #这样做的结果是返回两个3×3的随机矩阵并且逐元素比较取最大值。max(A,[],1) #这样做会得到每一列的最大值。max(A,[],2) #这将得到每一行的最大值。sum(sum(A.*eye(9)hold on #函数的功能是将新的图像绘制在旧的之上。再加上命令xlabel(time) 来标记X轴即水平轴输入ylabel(value)来标记垂直轴的值。legend(sin,cos) #表示这两条曲线表示的内容。title(myplot) #显示这幅图的标题。print(gcf,-dpng,abc.png) #保存为png格式的图片到当前路径squeeze() #用于删除矩阵中的单一维,对二维矩阵无效。close会让这个图像关掉。键入figure(1); plot(t, y1);将显示第一张图键入figure(2); plot(t, y2); 将显示第二张图(同时显示两张图)subplot(1,2,1)它将图像分为一个1*2的格子也就是前两个参数然后它使用第一个格子也就是最后一个参数1的意思axis([0.5 1 -1 1])也就是设置了轴的范围。横轴的范围调整至0.5到1竖轴的范围为-1到1。Clf清除一幅图像。我有时用一个巧妙的方法来可视化矩阵也就是imagesc(A)命令我还可以使用函数colorbar让我用一个更复杂的命令 imagesc(A)colorbarcolormap gray 它生成了一个颜色图像一个灰度分布图并在右边也加入一个颜色条。所以这个颜色条显示不同深浅的颜色所对应的值。csi_trace read_bf_file(csi6.dat);%数据读取for l1:1 %取50个数据包的数据csiaget_scaled_csi(csi_trace{l});%提取csi矩阵for i1:1 %1个发射天线for j1:3 %2个接收天线for k1:30 %30个子载波数据B(i,j,k)csia(i,j,k);endendendplot(db(abs(squeeze(B).)))%squeeze通过移除第一个单维度将csi变成3*30的矩阵%db将线性空间变成以十为底的对数空间%.转置得到30*3的矩阵%plot(Y)如果Y是m×n的数组以1:m为X横坐标Y中的每一列元素为Y坐标绘制n条曲线hold on%当前轴及图像保持而不被刷新准备接受此后将绘制的图形多图共存endlegend(RX Antenna A, RX Antenna B,RX Antenna C, Location, SouthEast );xlabel(Subcarrier index);ylabel(SNR [dB])hold off;%get_eff_SNRs() 它接受一个CSI矩阵作为输入并返回一个线性(幂)空间中有效信噪比值的74矩阵。%这4列对应的有效信噪比使用四种802.11调制方案即BPSK/QPSK/16QAM/64QAM。%7行对应于7个可能的天线选择包括3个天线和1、2或3个空间流。%特别是前3行对应于天线A、B或c的单流传输后3行对应于天线AB、AC或BC的双流传输。最后一行对应于使用所有天线的3流传输
