如果将c变为c=1:15,将此范围内所有max_non_diag_MAC 的值写入一个excel文件内,matlab

时间: 2024-04-12 14:30:40 浏览: 190
如果你想将范围内所有的`max_non_diag_MAC`值写入一个Excel文件中,你可以使用MATLAB中的`xlswrite`函数。下面是一个示例代码,演示如何将每个`c`值对应的`max_non_diag_MAC`值写入Excel文件: ```matlab [data1, text, raw] = xlsread('振型','Sheet1'); ZX = data1; m = size(ZX, 1); c_range = 1:15; max_non_diag_MAC_values = zeros(1, numel(c_range)); for i = 1:numel(c_range) c = c_range(i); d = [1:m]; for j = 1:m-c E = ZX * (ZX' * ZX)^(-1) * ZX'; EA = diag(E); [a, b] = min(EA); ZX(b, :) = []; d(:, b) = []; end M = ZX; [m, n] = size(M); MAC = zeros(n, n); for k = 1:n for l = 1:n a = M(:, k); b = M(:, l); MAC(k, l) = (a' * b)^2 / ((a' * a) * (b' * b)); end end diag_MAC = diag(MAC); non_diag_MAC = MAC - diag(diag_MAC); max_non_diag_MAC = max(non_diag_MAC(:)); max_non_diag_MAC_values(i) = max_non_diag_MAC; end xlswrite('max_non_diag_MAC_values.xlsx', max_non_diag_MAC_values'); ``` 在修改后的代码中,我们首先定义了变量`c_range`,它是一个包含1到15的范围。然后,我们使用`for`循环遍历每个`c`的值。在每次迭代中,我们都重新初始化`d`向量,并根据当前的`c`值进行循环迭代删除样本点的操作。 然后,我们计算了`max_non_diag_MAC`值,并将其存储在`max_non_diag_MAC_values`向量中。 最后,我们使用`xlswrite`函数将`max_non_diag_MAC_values`向量写入名为"max_non_diag_MAC_values.xlsx"的Excel文件中。 请确保在运行代码之前已经安装了MATLAB的Excel支持工具箱,并将Excel文件保存在当前工作目录或指定正确的文件路径。 如果你有任何进一步的问题,请随时提问。
阅读全文

相关推荐

clear all; close all; clc; tic bits_options = [0,1,2]; noise_option = 1; b = 4; NT = 2; SNRdBs =[0:2:20]; sq05=sqrt(0.5); nobe_target = 500; BER_target = 1e-3; raw_bit_len = 2592-6; interleaving_num = 72; deinterleaving_num = 72; N_frame = 1e8; for i_bits=1:length(bits_options) bits_option=bits_options(i_bits); BER=zeros(size(SNRdBs)); for i_SNR=1:length(SNRdBs) sig_power=NT; SNRdB=SNRdBs(i_SNR); sigma2=sig_power10^(-SNRdB/10)noise_option; sigma1=sqrt(sigma2/2); nobe = 0; Viterbi_init for i_frame=1:1:N_frame switch (bits_option) case {0}, bits=zeros(1,raw_bit_len); case {1}, bits=ones(1,raw_bit_len); case {2}, bits=randi(1,raw_bit_len,[0,1]); end encoding_bits = convolution_encoder(bits); interleaved=[]; for i=1:interleaving_num interleaved=[interleaved encoding_bits([i:interleaving_num:end])]; end temp_bit =[]; for tx_time=1:648 tx_bits=interleaved(1:8); interleaved(1:8)=[]; QAM16_symbol = QAM16_mod(tx_bits, 2); x(1,1) = QAM16_symbol(1); x(2,1) = QAM16_symbol(2); if rem(tx_time-1,81)==0 H = sq05(randn(2,2)+jrandn(2,2)); end y = Hx; if noise_option==1 noise = sqrt(sigma2/2)(randn(2,1)+j*randn(2,1)); y = y + noise; end W = inv(H'H+sigma2diag(ones(1,2)))H'; X_tilde = Wy; X_hat = QAM16_slicer(X_tilde, 2); temp_bit = [temp_bit QAM16_demapper(X_hat, 2)]; end deinterleaved=[]; for i=1:deinterleaving_num deinterleaved=[deinterleaved temp_bit([i:deinterleaving_num:end])]; end received_bit=Viterbi_decode(deinterleaved); for EC_dummy=1:1:raw_bit_len, if bits(EC_dummy)~=received_bit(EC_dummy), nobe=nobe+1; end if nobe>=nobe_target, break; end end if (nobe>=nobe_target) break; end end = BER(i_SNR) = nobe/((i_frame-1)*raw_bit_len+EC_dummy); fprintf('bits_option:%d,SNR:%d dB,BER:%1.4f\n',bits_option,SNRdB,BER(i_SNR)); end figure; semilogy(SNRdBs,BER); xlabel('SNR(dB)'); ylabel('BER'); title(['Bits_option:',num2str(bits_option)]); grid on; end将这段代码改为有噪声的情况

clear all; close all; clc; tic bits_options = [0,1,2]; noise_option = 1; b = 4; NT = 2; SNRdBs =[0:2:20]; sq05=sqrt(0.5); nobe_target = 500; BER_target = 1e-3; raw_bit_len = 2592-6; interleaving_num = 72; deinterleaving_num = 72; N_frame = 1e8; for i_bits=1:length(bits_options) bits_option=bits_options(i_bits); BER=zeros(size(SNRdBs)); for i_SNR=1:length(SNRdBs) sig_power=NT; SNRdB=SNRdBs(i_SNR); sigma2=sig_power*10^(-SNRdB/10)*noise_option; sigma1=sqrt(sigma2/2); nobe = 0; Viterbi_init for i_frame=1:1:N_frame switch (bits_option) case {0}, bits=zeros(1,raw_bit_len); case {1}, bits=ones(1,raw_bit_len); case {2}, bits=randi(1,raw_bit_len,[0,1]); end encoding_bits = convolution_encoder(bits); interleaved=[]; for i=1:interleaving_num interleaved=[interleaved encoding_bits([i:interleaving_num:end])]; end temp_bit =[]; for tx_time=1:648 tx_bits=interleaved(1:8); interleaved(1:8)=[]; QAM16_symbol = QAM16_mod(tx_bits, 2); x(1,1) = QAM16_symbol(1); x(2,1) = QAM16_symbol(2); if rem(tx_time-1,81)==0 H = sq05*(randn(2,2)+j*randn(2,2)); end y = H*x; if noise_option==1 noise = sqrt(sigma2/2)*(randn(2,1)+j*randn(2,1)); y = y + noise; end W = inv(H'*H+sigma2*diag(ones(1,2)))*H'; X_tilde = W*y; X_hat = QAM16_slicer(X_tilde, 2); temp_bit = [temp_bit QAM16_demapper(X_hat, 2)]; end deinterleaved=[]; for i=1:deinterleaving_num deinterleaved=[deinterleaved temp_bit([i:deinterleaving_num:end])]; end received_bit=Viterbi_decode(deinterleaved); for EC_dummy=1:1:raw_bit_len, if bits(EC_dummy)~=received_bit(EC_dummy), nobe=nobe+1; end if nobe>=nobe_target, break; end end if (nobe>=nobe_target) break; end end = BER(i_SNR) = nobe/((i_frame-1)*raw_bit_len+EC_dummy); fprintf('bits_option:%d,SNR:%d dB,BER:%1.4f\n',bits_option,SNRdB,BER(i_SNR)); end figure; semilogy(SNRdBs,BER); xlabel('SNR(dB)'); ylabel('BER'); title(['Bits_option:',num2str(bits_option)]); grid on; end

clear all; close all; clc;ticits_option = 2;noise_option = 1;raw_bit_len = 2592-6;interleaving_num = 72;deinterleaving_num = 72;N_frame = 1e4;SNRdBs = [0:2:20];sq05 = sqrt(0.5);bits_options = [0, 1, 2]; % 三种bits-option情况obe_target = 500;BER_target = 1e-3;for i_bits = 1:length(bits_options) bits_option = bits_options(i_bits); BER = zeros(size(SNRdBs)); for i_SNR = 1:length(SNRdBs) sig_power = 1; SNRdB = SNRdBs(i_SNR); sigma2 = sig_power * 10^(-SNRdB/10); sigma = sqrt(sigma2/2); nobe = 0; for i_frame = 1:N_frame switch bits_option case 0 bits = zeros(1, raw_bit_len); case 1 bits = ones(1, raw_bit_len); case 2 bits = randi([0,1], 1, raw_bit_len); end encoding_bits = convolution_encoder(bits); interleaved = []; for i = 1:interleaving_num interleaved = [interleaved encoding_bits([i:interleaving_num:end])]; end temp_bit = []; for tx_time = 1:648 tx_bits = interleaved(1:8); interleaved(1:8) = []; QAM16_symbol = QAM16_mod(tx_bits, 2); x(1,1) = QAM16_symbol(1); x(2,1) = QAM16_symbol(2); if rem(tx_time - 1, 81) == 0 H = sq05 * (randn(2,2) + j * randn(2,2)); end y = H * x; if noise_option == 1 noise = sigma * (randn(2,1) + j * randn(2,1)); y = y + noise; end W = inv(H' * H + sigma2 * diag(ones(1,2))) * H'; K_tilde = W * y; x_hat = QAM16_slicer(K_tilde, 2); temp_bit = [temp_bit QAM16_demapper(x_hat, 2)]; end deinterleaved = []; for i = 1:deinterleaving_num deinterleaved = [deinterleaved temp_bit([i:deinterleaving_num:end])]; end received_bit = Viterbi_decode(deinterleaved); for EC_dummy = 1:1:raw_bit_len if nobe >= obe_target break; end if received_bit(EC_dummy) ~= bits(EC_dummy) nobe = nobe + 1; end end if nobe >= obe_target break; end end BER(i_SNR) = nobe / (i_frame * raw_bit_len); fprintf('bits-option: %d, SNR: %d dB, BER: %1.4f\n', bits_option, SNRdB, BER(i_SNR)); end figure; semilogy(SNRdBs, BER); xlabel('SNR (dB)'); ylabel('BER'); title(['Bits-Option: ', num2str(bits_option)]); grid on;end注释这段matlab代码

最新推荐

recommend-type

1基于蓝牙的项目开发--蓝牙温度监测器.docx

1基于蓝牙的项目开发--蓝牙温度监测器.docx
recommend-type

AppDynamics:性能瓶颈识别与优化.docx

AppDynamics:性能瓶颈识别与优化
recommend-type

Haskell编写的C-Minus编译器针对TM架构实现

资源摘要信息:"cminus-compiler是一个用Haskell语言编写的C-Minus编程语言的编译器项目。C-Minus是一种简化版的C语言,通常作为教学工具使用,帮助学生了解编程语言和编译器的基本原理。该编译器的目标平台是虚构的称为TM的体系结构,尽管它并不对应真实存在的处理器架构,但这样的设计可以专注于编译器的逻辑而不受特定硬件细节的限制。作者提到这个编译器是其编译器课程的作业,并指出代码可以在多个方面进行重构,尽管如此,他对于编译器的完成度表示了自豪。 在编译器项目的文档方面,作者提供了名为doc/report1.pdf的文件,其中可能包含了关于编译器设计和实现的详细描述,以及如何构建和使用该编译器的步骤。'make'命令在简单的使用情况下应该能够完成所有必要的构建工作,这意味着项目已经设置好了Makefile文件来自动化编译过程,简化用户操作。 在Haskell语言方面,该编译器项目作为一个实际应用案例,可以作为学习Haskell语言特别是其在编译器设计中应用的一个很好的起点。Haskell是一种纯函数式编程语言,以其强大的类型系统和惰性求值特性而闻名。这些特性使得Haskell在处理编译器这种需要高度抽象和符号操作的领域中非常有用。" 知识点详细说明: 1. C-Minus语言:C-Minus是C语言的一个简化版本,它去掉了许多C语言中的复杂特性,保留了基本的控制结构、数据类型和语法。通常用于教学目的,以帮助学习者理解和掌握编程语言的基本原理以及编译器如何将高级语言转换为机器代码。 2. 编译器:编译器是将一种编程语言编写的源代码转换为另一种编程语言(通常为机器语言)的软件。编译器通常包括前端(解析源代码并生成中间表示)、优化器(改进中间表示的性能)和后端(将中间表示转换为目标代码)等部分。 3. TM体系结构:在这个上下文中,TM可能是一个虚构的计算机体系结构。它可能被设计来模拟真实处理器的工作原理,但不依赖于任何特定硬件平台的限制,有助于学习者专注于编译器设计本身,而不是特定硬件的技术细节。 4. Haskell编程语言:Haskell是一种高级的纯函数式编程语言,它支持多种编程范式,包括命令式、面向对象和函数式编程。Haskell的强类型系统、模式匹配、惰性求值等特性使得它在处理抽象概念如编译器设计时非常有效。 5. Make工具:Make是一种构建自动化工具,它通过读取Makefile文件来执行编译、链接和清理等任务。Makefile定义了编译项目所需的各种依赖关系和规则,使得项目构建过程更加自动化和高效。 6. 编译器开发:编译器的开发涉及语言学、计算机科学和软件工程的知识。它需要程序员具备对编程语言语法和语义的深入理解,以及对目标平台架构的了解。编译器通常需要进行详细的测试,以确保它能够正确处理各种边缘情况,并生成高效的代码。 通过这个项目,学习者可以接触到编译器从源代码到机器代码的转换过程,学习如何处理词法分析、语法分析、语义分析、中间代码生成、优化和目标代码生成等编译过程的关键步骤。同时,该项目也提供了一个了解Haskell语言在编译器开发中应用的窗口。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

【数据整理秘籍】:R语言与tidyr包的高效数据处理流程

![【数据整理秘籍】:R语言与tidyr包的高效数据处理流程](https://www.lecepe.fr/upload/fiches-formations/visuel-formation-246.jpg) # 1. 数据整理的重要性与R语言介绍 数据整理是数据科学领域的核心环节之一,对于后续的数据分析、模型构建以及决策制定起到至关重要的作用。高质量的数据整理工作,能够保证数据分析的准确性和可靠性,为数据驱动的业务决策提供坚实的数据基础。 在众多数据分析工具中,R语言因其强大的统计分析能力、丰富的数据处理包以及开放的社区支持而广受欢迎。R语言不仅仅是一种编程语言,它更是一个集数据处理、统
recommend-type

在使用STEP7编程环境为S7-300 PLC进行编程时,如何正确分配I/O接口地址并利用SM信号模板进行编址?

在西门子STEP7编程环境中,对于S7-300系列PLC的I/O接口地址分配及使用SM信号模板的编址是一个基础且至关重要的步骤。正确地进行这一过程可以确保PLC与现场设备之间的正确通信和数据交换。以下是具体的设置步骤和注意事项: 参考资源链接:[PLC STEP7编程环境:菜单栏与工具栏功能详解](https://wenku.csdn.net/doc/3329r82jy0?spm=1055.2569.3001.10343) 1. **启动SIMATIC Manager**:首先,启动STEP7软件,并通过SIMATIC Manager创建或打开一个项目。 2. **硬件配置**:在SIM
recommend-type

水电模拟工具HydroElectric开发使用Matlab

资源摘要信息:"该文件是一个使用MATLAB开发的水电模拟应用程序,旨在帮助用户理解和模拟HydroElectric实验。" 1. 水电模拟的基础知识: 水电模拟是一种利用计算机技术模拟水电站的工作过程和性能的工具。它可以模拟水电站的水力、机械和电气系统,以及这些系统的相互作用和影响。水电模拟可以帮助我们理解水电站的工作原理,预测和优化其性能,以及评估和制定运行策略。 2. MATLAB在水电模拟中的应用: MATLAB是一种高性能的数值计算和可视化软件,广泛应用于工程、科学和数学领域。在水电模拟中,MATLAB可以用于建立模型、模拟、分析和可视化水电站的性能。MATLAB提供了强大的数学函数库和图形工具箱,可以方便地进行复杂的计算和数据可视化。 3. HydroElectric实验的模拟: HydroElectric实验是一种模拟水电站工作的实验,通常包括水轮机、发电机、水道、负荷等部分。在这个实验中,我们可以模拟各种运行条件下的水电站性能,如不同水流量、不同负荷等。 4. MATLAB开发的水电模拟应用程序的使用: 使用MATLAB开发的水电模拟应用程序,用户可以方便地设置模拟参数,运行模拟,查看模拟结果。应用程序可能包括用户友好的界面,用户可以通过界面输入各种参数,如水流量、负荷等。然后,应用程序将根据输入的参数,进行计算,模拟水电站的工作过程和性能,最后将结果以图表或数据的形式展示给用户。 5. MATLAB的高级功能在水电模拟中的应用: MATLAB提供了丰富的高级功能,如优化工具箱、神经网络工具箱、符号计算等,这些功能可以进一步提高水电模拟的效果。例如,使用优化工具箱,我们可以找到最佳的工作参数,使水电站的性能最优化。使用神经网络工具箱,我们可以建立更复杂的模型,更准确地模拟水电站的工作过程。使用符号计算,我们可以处理更复杂的数学问题,如求解非线性方程。 6. 水电模拟的未来发展方向: 随着计算机技术的不断发展,水电模拟的应用前景广阔。未来,水电模拟可能会更加注重模型的精确度和复杂度,更多地运用人工智能、大数据等先进技术,以提高模拟的效率和准确性。此外,水电模拟也可能更多地应用于其他领域,如能源管理、环境影响评估等。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

【数据分析必修课】:R语言中tidyr包的终极使用指南

![【数据分析必修课】:R语言中tidyr包的终极使用指南](https://study.com/cimages/videopreview/ewh840ozgx.jpg) # 1. R语言与数据分析基础 ## 1.1 R语言简介 R语言是一种专门用于统计分析和图形表示的编程语言。它以其自由开源的特性、强大的数据处理能力以及丰富的社区支持著称。无论您是初学者还是有经验的数据分析师,R语言都提供了一个灵活的平台来探索数据,生成报告,或创建复杂的数据模型。 ## 1.2 数据分析基础 数据分析是在一系列数据上进行的系统性的研究过程,目的是提取有用信息、发现模式、验证假设,以及支持决策。数据分析通
recommend-type

在机器学习项目中,如何采用可解释性技术来提升文本分类模型的透明度,并确保模型解释性?

在机器学习项目中,尤其是在文本分类任务中,确保模型的透明度和解释性是一个重要议题。《可解释的机器学习:深入理解黑盒模型》这本书为解决这一问题提供了丰富的理论和实践指导。以下是一些关键步骤和建议: 参考资源链接:[可解释的机器学习:深入理解黑盒模型](https://wenku.csdn.net/doc/41nji7cnbf?spm=1055.2569.3001.10343) 首先,理解和掌握机器学习模型的基础知识是非常重要的。这包括模型是如何从输入数据中学习规律并作出预测的,以及与可解释性相关的术语。在文本分类任务中,理解诸如TF-IDF、词嵌入、卷积神经网络(CNN)或循环神经网络(RN