将以下代码转化为matlab代码表示:import xlrd import sympy import numpy as np from scipy import linalg #%% queue = [ 0, 29, 17, 2, 1, 20, 19, 26, 18, 25, 14, 6, 11, 7, 15, 9, 8, 12, 27, 16, 10, 13, 5, 4, 3, 22, 28, 24, 23, 21, 0] def read_data_model(): data = xlrd.open_workbook("/Users/lzs/Downloads/2020szcupc/data/C2.xlsx") table = data.sheet_by_name("Sheet1") rowNum = table.nrows colNum = table.ncols consumes = [] for i in range(1, rowNum): # 忽略DC的消耗 if i == 1: pass else: consumes.append(0 if table.cell_value(i, 3) == '/' else table.cell_value(i, 3)) return consumes #%% 获得矩阵A def get_A_matrix(data): A = np.ones([29,29], dtype = float) diagonal = np.eye(29) for i in range(29): for j in range(29): A[i][j] = data['consumes'][j] / data['r'] A = A - diagonal return A #%% def get_b_maatrix(data): b = np.ones([29,1], dtype=float) for i in range(29): b[i][0] = -data['dst']*data['consumes'][i]/data['velocity']+data['f'] for j in range(29): b[i][0] = b[i][0] + data['f']*data['consumes'][i]/data['r'] return b #%% 数值解 def numerical(data): data['velocity'] = 50 data['dst'] = 11469 data['r'] = 200 data['f'] = 10 A = get_A_matrix(data) b = get_b_maatrix(data) x = linalg.solve(A, b) return x #%% 符号解决方案 def symbolic(data): data['velocity'] = sympy.symbols("v", integer = True) data['dst'] = 12100 data['r'] = sympy.symbols("r", integer = True) data['f'] = sympy.symbols("f", integer = True) # 获取矩阵A并转移到符号矩阵M A = np.ones([29,29], dtype = float).tolist() diagonal = np.eye(29).tolist() for i in range(29): for j in range(29): A[i][j] = data['consumes'][j] / data['r'] - diagonal[i][j] M = sympy.Matrix(A) # 得到矩阵b并转移到符号矩阵b b = np.ones([29,1], dtype=float).tolist() for i in range(29): b[i][0] = -data['dst']*data['consumes'][i]/data['velocity']+data['f'] for j in range(29): b[i][0] = b[i][0] + data['f']*data['consumes'][i]/data['r'] b = sympy.Matrix(b) # LU solver x = M.LUsolve(b) return x #%% 主功能 if name == 'main': data = {} data['consumes'] = read_data_model() options = {"numerical":1, "symbolic":2} option = 1 if option == options['numerical']: x = numerical(data) print(x) elif option == options['symbolic']: x = symbolic(data) print(x) else: print("WARN!!!")

时间: 2023-10-09 16:15:01 浏览: 182
以下是将给定代码转化为 Matlab 代码的结果: ```matlab %% 导入所需库 import matlab.io.* import sympy.* import numpy as np from scipy import linalg %% 读取数据模型 function consumes = read_data_model() data = xlsread('/Users/lzs/Downloads/2020szcupc/data/C2.xlsx'); consumes = []; for i = 2:size(data,1) if i == 2 continue % 忽略DC的消耗 end consumes = [consumes, 0]; if data(i,4) ~= '/' consumes(end) = data(i,4); end end end %% 获取矩阵A function A = get_A_matrix(data) A = ones(29,29); diagonal = eye(29); for i = 1:29 for j = 1:29 A(i,j) = data.consumes(j) / data.r; end A(i,i) = A(i,i) - 1; end end %% 获取矩阵b function b = get_b_maatrix(data) b = ones(29,1); for i = 1:29 b(i) = -data.dst*data.consumes(i)/data.velocity + data.f; for j = 1:29 b(i) = b(i) + data.f*data.consumes(j)/data.r; end end end %% 数值解 function x = numerical(data) data.velocity = 50; data.dst = 11469; data.r = 200; data.f = 10; A = get_A_matrix(data); b = get_b_maatrix(data); x = linalg.solve(A, b); end %% 符号解决方案 function x = symbolic(data) syms v r f data.velocity = v; data.dst = 12100; data.r = r; data.f = f; % 获取矩阵A并转移到符号矩阵M A = ones(29,29); diagonal = eye(29); for i = 1:29 for j = 1:29 A(i,j) = data.consumes(j) / data.r; end A(i,i) = A(i,i) - 1; end M = sym(A); % 得到矩阵b并转移到符号矩阵b b = ones(29,1); for i = 1:29 b(i) = -data.dst*data.consumes(i)/data.velocity + data.f; for j = 1:29 b(i) = b(i) + data.f*data.consumes(j)/data.r; end end b = sym(b); % LU求解 x = M\LUsolve(b); end %% 主功能 data = struct(); data.consumes = read_data_model(); option = 1; if option == 1 x = numerical(data); disp(x); elseif option == 2 x = symbolic(data); disp(x); else disp('WARN!!!'); end ``` 需要注意的是,Matlab 中没有 `xlrd` 库,需要使用 `xlsread` 函数来读取 Excel 文件。另外,Matlab 中的 `eye` 函数可以直接得到单位矩阵,无需转化为列表再使用 `diag` 函数。
阅读全文

相关推荐

import numpy as np import xlrd import matplotlib.pyplot as plt from sklearn.feature_selection import RFE from sklearn.ensemble import RandomForestClassifier from sklearn.svm import SVC from sklearn.linear_model import LogisticRegression from sklearn.model_selection import cross_val_score def excel2m(path):#读excel数据转为矩阵函数 data = xlrd.open_workbook(path) table = data.sheets()[0] # 获取excel中第一个sheet表 nrows = table.nrows # 行数 ncols = table.ncols # 列数 datamatrix = np.zeros((nrows, ncols)) for x in range(ncols): cols = table.col_values(x) cols1 = np.matrix(cols) # 把list转换为矩阵进行矩阵操作 datamatrix[:, x] = cols1 # 把数据进行存储 return datamatrix x=excel2m("factors.xlsx") x=np.matrix(x) y=excel2m("RON.xlsx") y=np.matrix(y) rfc=RandomForestClassifier(n_estimators=10,random_state=0) score=[] for i in range(1,200,10): rfe = RFE(estimator=rfc, n_features_to_select=i, step=10).fit(x, y.astype('int')) rfe.support_.sum() rfe.ranking_ x_wrapper=rfe.transform(x) once=cross_val_score(rfc,x_wrapper,y.astype('int'),cv=5).mean() score.append(once) plt.figure(figsize=[20,5]) plt.plot(range(1,200,10),score) plt.xticks(range(1,200,10)) plt.show() np.savetxt('score.csv', score, delimiter = ',') # 确定选择特征数量后,看各个特征得分排名 # 每个特征的得分排名,特征得分越低(1最好),表示特征越好 #print(rfe.ranking_) #np.savetxt('ranking.csv', rfe.ranking_, delimiter = ',') # 每次交叉迭代各个特征得分 #print(rfe.grid_scores_) #np.savetxt('grid_scores.csv', rfe.grid_scores_, delimiter = ',')

最新推荐

recommend-type

ABAQUS2016常用库函数Numpy,Scipy,Xlrd及Xlwt,Reportlab安装.docx

为此,我们需要安装一些常用的库函数,包括 Numpy、Scipy、Xlrd、Xlwt 和 Reportlab 等。以下是安装这些库函数的详细步骤。 首先,我们需要查看 ABAQUS2016 中安装的是哪个 Python 版本。可以通过点击 ABAQUS ...
recommend-type

使用matlab或python将txt文件转为excel表格

这段代码首先通过`importdata`读取名为"data.txt"的文本文件,然后使用`xlswrite`函数将读取到的数据写入一个名为"data.xls"的Excel文件中。 接下来,我们转向Python。Python在数据处理方面有pandas库,它提供了一...
recommend-type

在pycharm中导入xlrd和xlwt.模块具体操作.docx

当涉及到Excel数据处理时,`xlrd`和`xlwt`是两个不可或缺的库,它们分别用于读取和写入Excel文件。这篇文档详细介绍了如何在PyCharm中安装和使用这两个模块。 首先,你需要确保已经安装了`xlrd`和`xlwt`。在命令行...
recommend-type

java计算器源码.zip

java毕业设计源码,可供参考
recommend-type

PHP集成Autoprefixer让CSS自动添加供应商前缀

标题和描述中提到的知识点主要包括:Autoprefixer、CSS预处理器、Node.js 应用程序、PHP 集成以及开源。 首先,让我们来详细解析 Autoprefixer。 Autoprefixer 是一个流行的 CSS 预处理器工具,它能够自动将 CSS3 属性添加浏览器特定的前缀。开发者在编写样式表时,不再需要手动添加如 -webkit-, -moz-, -ms- 等前缀,因为 Autoprefixer 能够根据各种浏览器的使用情况以及官方的浏览器版本兼容性数据来添加相应的前缀。这样可以大大减少开发和维护的工作量,并保证样式在不同浏览器中的一致性。 Autoprefixer 的核心功能是读取 CSS 并分析 CSS 规则,找到需要添加前缀的属性。它依赖于浏览器的兼容性数据,这一数据通常来源于 Can I Use 网站。开发者可以通过配置文件来指定哪些浏览器版本需要支持,Autoprefixer 就会自动添加这些浏览器的前缀。 接下来,我们看看 PHP 与 Node.js 应用程序的集成。 Node.js 是一个基于 Chrome V8 引擎的 JavaScript 运行时环境,它使得 JavaScript 可以在服务器端运行。Node.js 的主要特点是高性能、异步事件驱动的架构,这使得它非常适合处理高并发的网络应用,比如实时通讯应用和 Web 应用。 而 PHP 是一种广泛用于服务器端编程的脚本语言,它的优势在于简单易学,且与 HTML 集成度高,非常适合快速开发动态网站和网页应用。 在一些项目中,开发者可能会根据需求,希望把 Node.js 和 PHP 集成在一起使用。比如,可能使用 Node.js 处理某些实时或者异步任务,同时又依赖 PHP 来处理后端的业务逻辑。要实现这种集成,通常需要借助一些工具或者中间件来桥接两者之间的通信。 在这个标题中提到的 "autoprefixer-php",可能是一个 PHP 库或工具,它的作用是把 Autoprefixer 功能集成到 PHP 环境中,从而使得在使用 PHP 开发的 Node.js 应用程序时,能够利用 Autoprefixer 自动处理 CSS 前缀的功能。 关于开源,它指的是一个项目或软件的源代码是开放的,允许任何个人或组织查看、修改和分发原始代码。开源项目的好处在于社区可以一起参与项目的改进和维护,这样可以加速创新和解决问题的速度,也有助于提高软件的可靠性和安全性。开源项目通常遵循特定的开源许可证,比如 MIT 许可证、GNU 通用公共许可证等。 最后,我们看到提到的文件名称 "autoprefixer-php-master"。这个文件名表明,该压缩包可能包含一个 PHP 项目或库的主分支的源代码。"master" 通常是源代码管理系统(如 Git)中默认的主要分支名称,它代表项目的稳定版本或开发的主线。 综上所述,我们可以得知,这个 "autoprefixer-php" 工具允许开发者在 PHP 环境中使用 Node.js 的 Autoprefixer 功能,自动为 CSS 规则添加浏览器特定的前缀,从而使得开发者可以更专注于内容的编写而不必担心浏览器兼容性问题。
recommend-type

揭秘数字音频编码的奥秘:非均匀量化A律13折线的全面解析

# 摘要 数字音频编码技术是现代音频处理和传输的基础,本文首先介绍数字音频编码的基础知识,然后深入探讨非均匀量化技术,特别是A律压缩技术的原理与实现。通过A律13折线模型的理论分析和实际应用,本文阐述了其在保证音频信号质量的同时,如何有效地降低数据传输和存储需求。此外,本文还对A律13折线的优化策略和未来发展趋势进行了展望,包括误差控制、算法健壮性的提升,以及与新兴音频技术融合的可能性。 # 关键字 数字音频编码;非均匀量化;A律压缩;13折线模型;编码与解码;音频信号质量优化 参考资源链接:[模拟信号数字化:A律13折线非均匀量化解析](https://wenku.csdn.net/do
recommend-type

arduino PAJ7620U2

### Arduino PAJ7620U2 手势传感器 教程 #### 示例代码与连接方法 对于Arduino开发PAJ7620U2手势识别传感器而言,在Arduino IDE中的项目—加载库—库管理里找到Paj7620并下载安装,完成后能在示例里找到“Gesture PAJ7620”,其中含有两个示例脚本分别用于9种和15种手势检测[^1]。 关于连线部分,仅需连接四根线至Arduino UNO开发板上的对应位置即可实现基本功能。具体来说,这四条线路分别为电源正极(VCC),接地(GND),串行时钟(SCL)以及串行数据(SDA)[^1]。 以下是基于上述描述的一个简单实例程序展示如
recommend-type

网站啄木鸟:深入分析SQL注入工具的效率与限制

网站啄木鸟是一个指的是一类可以自动扫描网站漏洞的软件工具。在这个文件提供的描述中,提到了网站啄木鸟在发现注入漏洞方面的功能,特别是在SQL注入方面。SQL注入是一种常见的攻击技术,攻击者通过在Web表单输入或直接在URL中输入恶意的SQL语句,来欺骗服务器执行非法的SQL命令。其主要目的是绕过认证,获取未授权的数据库访问权限,或者操纵数据库中的数据。 在这个文件中,所描述的网站啄木鸟工具在进行SQL注入攻击时,构造的攻击载荷是十分基础的,例如 "and 1=1--" 和 "and 1>1--" 等。这说明它的攻击能力可能相对有限。"and 1=1--" 是一个典型的SQL注入载荷示例,通过在查询语句的末尾添加这个表达式,如果服务器没有对SQL注入攻击进行适当的防护,这个表达式将导致查询返回真值,从而使得原本条件为假的查询条件变为真,攻击者便可以绕过安全检查。类似地,"and 1>1--" 则会检查其后的语句是否为假,如果查询条件为假,则后面的SQL代码执行时会被忽略,从而达到注入的目的。 描述中还提到网站啄木鸟在发现漏洞后,利用查询MS-sql和Oracle的user table来获取用户表名的能力不强。这表明该工具可能无法有效地探测数据库的结构信息或敏感数据,从而对数据库进行进一步的攻击。 关于实际测试结果的描述中,列出了8个不同的URL,它们是针对几个不同的Web应用漏洞扫描工具(Sqlmap、网站啄木鸟、SqliX)进行测试的结果。这些结果表明,针对提供的URL,Sqlmap和SqliX能够发现注入漏洞,而网站啄木鸟在多数情况下无法识别漏洞,这可能意味着它在漏洞检测的准确性和深度上不如其他工具。例如,Sqlmap在针对 "http://www.2cto.com/news.php?id=92" 和 "http://www.2cto.com/article.asp?ID=102&title=Fast food marketing for children is on the rise" 的URL上均能发现SQL注入漏洞,而网站啄木鸟则没有成功。这可能意味着网站啄木鸟的检测逻辑较为简单,对复杂或隐蔽的注入漏洞识别能力不足。 从这个描述中,我们也可以了解到,在Web安全测试中,工具的多样性选择是十分重要的。不同的安全工具可能对不同的漏洞和环境有不同的探测能力,因此在实际的漏洞扫描过程中,安全测试人员需要选择合适的工具组合,以尽可能地全面地检测出应用中存在的漏洞。 在标签中指明了这是关于“sql注入”的知识,这表明了文件主题的核心所在。SQL注入是一种常见的网络攻击方式,安全测试人员、开发人员和网络管理员都需要对此有所了解,以便进行有效的防御和检测。 最后,提到了压缩包子文件的文件名称列表,其中包含了三个文件:setup.exe、MD5.exe、说明_Readme.html。这里提供的信息有限,但可以推断setup.exe可能是一个安装程序,MD5.exe可能是一个计算文件MD5散列值的工具,而说明_Readme.html通常包含的是软件的使用说明或者版本信息等。这些文件名暗示了在进行网站安全测试时,可能涉及到安装相关的软件工具,以及进行文件的校验和阅读相应的使用说明。然而,这些内容与文件主要描述的web安全漏洞检测主题不是直接相关的。
recommend-type

【GPStoolbox使用技巧大全】:20个实用技巧助你精通GPS数据处理

# 摘要 GPStoolbox是一个广泛应用于GPS数据处理的软件工具箱,它提供了从数据导入、预处理、基本分析到高级应用和自动化脚本编写的全套功能。本文介绍了GPStoolbox的基本概况、安装流程以及核心功能,探讨了如何
recommend-type

spring boot怎么配置maven

### 如何在 Spring Boot 项目中正确配置 Maven #### pom.xml 文件设置 `pom.xml` 是 Maven 项目的核心配置文件,在 Spring Boot 中尤为重要,因为其不仅管理着所有的依赖关系还控制着项目的构建流程。对于 `pom.xml` 的基本结构而言,通常包含如下几个部分: - **Project Information**: 定义了关于项目的元数据,比如模型版本、组ID、工件ID和版本号等基本信息[^1]。 ```xml <project xmlns="http://maven.apache.org/POM/4.0.0