def signal_john(df, para=[20, 0.05, 20], proportion=1): para_ = para[0] df['median'] = df['close'].rolling(para_, min_periods=1).mean() df['bias'] = df['close'] / df['median'] - 1 bias_pct = float(para[1]) wr_pct = para[2] df['TYP'] = (df['high'] + df['low'] + df['close']) / 3 df['H'] = np.where(df['high'] - df['TYP'].shift(1) > 0, df['high'] - df['TYP'].shift(1), 0) df['L'] = np.where(df['TYP'].shift(1) - df['low'] > 0, df['TYP'].shift(1) - df['low'], 0) df['max_high'] = df['high'].rolling(wr_pct, min_periods=1).max() df['min_low'] = df['low'].rolling(wr_pct, min_periods=1).min() df['Wr_%s' % str(wr_pct)] = (df['max_high'] - df['close']) / (df['max_high'] - df['min_low']) * 100 # 计算均线 df['Cr_%s' % (para_)] = df['H'].rolling(para_).sum() / df['L'].rolling(para_, min_periods=1).sum() * 100 # ======= 找出做多信号 CR 上穿 200 condition1 = df['Cr_%s' % str(para_)] > 200 # 均线大于0 # condition2 = df['Cr_%s' % str(params)].shift(1) <= 200 # 上一周期的均线小于等于0 condition2 = df['Cr_%s' % str(para_)].shift() <= 200 condition3 = df['Wr_%s' % str(wr_pct)] < 20 df.loc[condition1 & condition2 & condition3, 'signal_long'] = 1 # 1代表做多 condition1 = df['Cr_%s' % str(para_)] < 50 # 均线大于0 # condition2 = df['Cr_%s' % str(params)].shift(1) <= 200 # 上一周期的均线小于等于0 condition2 = df['Cr_%s' % str(para_)].shift() >= 50 condition3 = df['Wr_%s' % str(wr_pct)] < 80 df.loc[condition1 & condition2 & condition3, 'signal_short'] = -1 # 1代表做多 # 合并做多做空信号,去除重复信号 df['signal'] = df[['signal_long', 'signal_short']].sum(axis=1, min_count=1, skipna=True) # 若你的pandas版本是最新的,请使用本行代码代替上面一行 temp = df[df['signal'].notnull()][['signal']] temp = temp[temp['signal'] != temp['signal'].shift(1)] df['signal'] = temp['signal'] # ===根据bias,修改开仓时间 df['temp'] = df['signal'] # 将原始信号做多时,当bias大于阀值,设置为空 condition1 = (df['signal'] == 1) condition2 = (df['bias'] > bias_pct) df.loc[condition1 & condition2, 'temp'] = None # 将原始信号做空时,当bias大于阀值,设置为空 condition1 = (df['signal'] == -1) condition2 = (df['bias'] < -1 * bias_pct) df.loc[condition1 & condition2, 'temp'] = None # 使用之前的信号补全原始信号 df['temp'].fillna(method='ffill', inplace=True) df['signal'] = df['temp'] # ===考察是否需要止盈止损 df = process_stop_loss_close(df, proportion) return df #### 3.2、参数生成代码 def generate_fibonacci_sequence(start, end): """ 生成斐波那契数列 :param start: 数列起始值 :param end: 数列结束值 :return: 从起始值到结束值的斐波那契数列 """ seq = [start, start + 1] while seq[-1] <= end: seq.append(seq[-1] + seq[-2]) return seq[:-1] 优化以上代码

时间: 2023-10-09 15:17:28 浏览: 40
以下是对代码的优化建议: 1. 可以使用 @numba.jit 对函数进行装饰,以提高函数的执行速度。 2. 可以使用 pandas 的 DataFrame.pipe() 方法,将多个操作串联起来,以提高代码的可读性和简洁性。 3. 可以将函数中的一些常量或变量提取出来,作为函数的参数,以增加函数的灵活性。比如,可以将 bias_pct 和 wr_pct 提取出来,作为函数的参数。 4. 可以将一些重复的代码进行封装,以减少代码的重复性。比如,可以将处理做多做空信号的代码封装成一个函数。 5. 可以将函数中的一些不必要的变量删除,以减少内存的占用。比如,可以将 df['temp'] 删除,直接使用 df['signal']。 6. 可以使用 pandas 的 rolling_apply() 方法,代替 rolling().apply() 方法,以提高代码的执行速度。 7. 可以使用 pandas 的 query() 方法,代替 loc[] 方法,以提高代码的执行速度和可读性。 8. 可以使用 numpy 的 vectorize() 方法,对函数进行向量化,以提高函数的执行速度。
相关问题

LOPA_rating_button = wx.Button(left_panel, label="SIL定级") LOPA_sizer.Add(LOPA_rating_button, flag=wx.ALL | wx.EXPAND, border=self.boder, proportion=1) 怎么给这段代码加凸起边框

你可以使用`wx.BoxSizer`来为按钮添加一个凸起的边框效果。可以按照以下步骤进行操作: 1. 首先,将按钮添加到一个新的`wx.Panel`中,以便您可以对其进行布局和样式设置。 ```python LOPA_rating_panel = wx.Panel(left_panel) LOPA_rating_sizer = wx.BoxSizer(wx.HORIZONTAL) LOPA_rating_panel.SetSizer(LOPA_rating_sizer) LOPA_rating_button = wx.Button(LOPA_rating_panel, label="SIL定级") LOPA_rating_sizer.Add(LOPA_rating_button, flag=wx.ALL | wx.EXPAND, border=self.border, proportion=1) ``` 2. 接下来,使用`wx.StaticBox`和`wx.StaticBoxSizer`来创建一个静态框,并将按钮添加到其中。 ```python LOPA_rating_staticbox = wx.StaticBox(left_panel, label="边框") LOPA_rating_staticbox_sizer = wx.StaticBoxSizer(LOPA_rating_staticbox, wx.VERTICAL) LOPA_rating_staticbox_sizer.Add(LOPA_rating_panel, flag=wx.EXPAND) LOPA_sizer.Add(LOPA_rating_staticbox_sizer, flag=wx.ALL | wx.EXPAND, border=self.border, proportion=1) ``` 请注意,我假设您已经有一个名为`LOPA_sizer`的`sizer`来管理左侧面板的布局。如果没有,请根据您的实际情况进行调整。 这样做会在按钮周围添加一个凸起的边框效果。您可以根据需要调整边框的样式和参数。

def calculate_importance(data, feature_index): unique_values = np.unique(data[:, feature_index]) # 属性取值的唯一值 importance = 0 for value in unique_values: subset = data[data[:

feature_index]] count = np.sum(subset[:, feature_index] == value) proportion = count / len(data) importance += -proportion * np.log2(proportion) return importance 这是一个计算给定数据集中某个特征的重要性的函数。它首先获取该特征的唯一值集合,然后遍历每个唯一值,计算该值在数据集中的比例,并根据比例计算重要性。最后将各个唯一值的重要性求和并返回。 你还有其他问题吗?

相关推荐

帮我解释一下 PID_TypeDef g_location_pid; /* 位置PID参数结构体*/ /** * @brief 初始化PID参数 * @param 无 * @retval 无 / void pid_init(void) { /位置环初始化/ g_location_pid.SetPoint = (float)(50PPM); /* 设定目标Desired Value*/ g_location_pid.ActualValue = 0.0; /* 期望值*/ g_location_pid.SumError = 0.0; /* 积分值*/ g_location_pid.Error = 0.0; /* Error[1]/ g_location_pid.LastError = 0.0; / Error[-1]/ g_location_pid.PrevError = 0.0; / Error[-2]/ g_location_pid.Proportion = L_KP; / 比例常数 Proportional Const*/ g_location_pid.Integral = L_KI; /* 积分常数 Integral Const*/ g_location_pid.Derivative = L_KD; /* 微分常数 Derivative Const*/ g_location_pid.IngMax = 20; g_location_pid.IngMin = -20; g_location_pid.OutMax = 150; /* 输出限制 / g_location_pid.OutMin = -150; } /* * 函数名称:位置闭环PID控制设计 * 输入参数:当前控制量 * 返 回 值:目标控制量 * 说 明:无 */ int32_t increment_pid_ctrl(PID_TypeDef PID,float Feedback_value) { PID->Error = (float)(PID->SetPoint - Feedback_value); / 偏差 / #if INCR_LOCT_SELECT PID->ActualValue += (PID->Proportion * (PID->Error - PID->LastError)) / E[k]项 / + (PID->Integral * PID->Error) / E[k-1]项 / + (PID->Derivative * (PID->Error - 2 * PID->LastError + PID->PrevError)); / E[k-2]项 / PID->PrevError = PID->LastError; / 存储误差,用于下次计算 / PID->LastError = PID->Error; #else PID->SumError += PID->Error; if(PID->SumError > PID->IngMax) { PID->SumError = PID->IngMax; } else if(PID->SumError < PID->IngMin) { PID->SumError = PID->IngMin; } PID->ActualValue = (PID->Proportion * PID->Error) / E[k]项 / + (PID->Integral * PID->SumError) / E[k-1]项 / + (PID->Derivative * (PID->Error - PID->LastError)); / E[k-2]项 / PID->LastError = PID->Error; #endif if(PID->ActualValue > PID->OutMax) { PID->ActualValue = PID->OutMax; } else if(PID->ActualValue < PID->OutMin) { PID->ActualValue = PID->OutMin; } return ((int32_t)(PID->ActualValue)); / 返回实际控制数值 */ }

import pandas as pd import numpy as np from sklearn.decomposition import PCA from sklearn.preprocessing import StandardScaler import matplotlib.pyplot as plt # 读取数据 data = pd.read_csv('D:/pythonProject/venv/BostonHousing2.csv') # 提取前13个指标的数据 X = data.iloc[:, 5:18].values # 数据标准化 scaler = StandardScaler() X_scaled = scaler.fit_transform(X) # 主成分分析 pca = PCA() X_pca = pca.fit_transform(X_scaled) # 特征值和特征向量 eigenvalues = pca.explained_variance_ eigenvectors = pca.components_.T # 碎石图 variance_explained = np.cumsum(eigenvalues / np.sum(eigenvalues)) plt.plot(range(6, 19), variance_explained, marker='o') plt.xlabel('Number of Components') plt.ylabel('Cumulative Proportion of Variance Explained') plt.title('Scree Plot') plt.show() # 选择主成分个数 n_components = np.sum(variance_explained <= 0.95) + 1 # 前2个主成分的载荷图 loadings = pd.DataFrame(eigenvectors[:, 0:2], columns=['PC1', 'PC2'], index=data.columns[0:13]) plt.figure(figsize=(10, 6)) plt.scatter(loadings['PC1'], loadings['PC2'], alpha=0.7) for i, feature in enumerate(loadings.index): plt.text(loadings['PC1'][i], loadings['PC2'][i], feature) plt.xlabel('PC1') plt.ylabel('PC2') plt.title('Loading Plot') plt.grid() plt.show() # 主成分得分图 scores = pd.DataFrame(X_pca[:, 0:n_components], columns=['PC{}'.format(i+1) for i in range(n_components)]) plt.figure(figsize=(10, 6)) plt.scatter(scores['PC1'], scores['PC2'], alpha=0.7) for i, label in enumerate(data['MEDV']): plt.text(scores['PC1'][i], scores['PC2'][i], label) plt.xlabel('PC1') plt.ylabel('PC2') plt.title('Scores Plot') plt.grid() plt.show() # 综合评估和排序 data['PC1_score'] = X_pca[:, 0] sorted_data = data.sort_values(by='PC1_score') # 主成分回归模型 from sklearn.linear_model import LinearRegression Y = data['MEDV'].values.reshape(-1, 1) X_pca_regression = X_pca[:, 0].reshape(-1, 1) regression_model = LinearRegression() regression_model.fit(X_pca_regression, Y) # 回归方程 intercept = regression_model.intercept_[0] slope = regression_model.coef_[0][0] equation = "MEDV = {:.2f} + {:.2f} * PC1".format(intercept, slope) print("Regression Equation:", equation) # 最小二乘估计结果 from statsmodels.api import OLS X_const = np.concatenate((np.ones((506, 1)), X_pca_regression), axis=1) ols_model = OLS(Y, X_const).fit() print("OLS Regression Summary:") print(ols_model.summary())

最新推荐

recommend-type

node-v18.11.0-headers.tar.xz

Node.js,简称Node,是一个开源且跨平台的JavaScript运行时环境,它允许在浏览器外运行JavaScript代码。Node.js于2009年由Ryan Dahl创立,旨在创建高性能的Web服务器和网络应用程序。它基于Google Chrome的V8 JavaScript引擎,可以在Windows、Linux、Unix、Mac OS X等操作系统上运行。 Node.js的特点之一是事件驱动和非阻塞I/O模型,这使得它非常适合处理大量并发连接,从而在构建实时应用程序如在线游戏、聊天应用以及实时通讯服务时表现卓越。此外,Node.js使用了模块化的架构,通过npm(Node package manager,Node包管理器),社区成员可以共享和复用代码,极大地促进了Node.js生态系统的发展和扩张。 Node.js不仅用于服务器端开发。随着技术的发展,它也被用于构建工具链、开发桌面应用程序、物联网设备等。Node.js能够处理文件系统、操作数据库、处理网络请求等,因此,开发者可以用JavaScript编写全栈应用程序,这一点大大提高了开发效率和便捷性。 在实践中,许多大型企业和组织已经采用Node.js作为其Web应用程序的开发平台,如Netflix、PayPal和Walmart等。它们利用Node.js提高了应用性能,简化了开发流程,并且能更快地响应市场需求。
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

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

实现实时数据湖架构:Kafka与Hive集成

![实现实时数据湖架构:Kafka与Hive集成](https://img-blog.csdnimg.cn/img_convert/10eb2e6972b3b6086286fc64c0b3ee41.jpeg) # 1. 实时数据湖架构概述** 实时数据湖是一种现代数据管理架构,它允许企业以低延迟的方式收集、存储和处理大量数据。与传统数据仓库不同,实时数据湖不依赖于预先定义的模式,而是采用灵活的架构,可以处理各种数据类型和格式。这种架构为企业提供了以下优势: - **实时洞察:**实时数据湖允许企业访问最新的数据,从而做出更明智的决策。 - **数据民主化:**实时数据湖使各种利益相关者都可
recommend-type

可见光定位LED及其供电硬件具体型号,广角镜头和探测器,实验设计具体流程步骤,

1. 可见光定位LED型号:一般可使用5mm或3mm的普通白色LED,也可以选择专门用于定位的LED,例如OSRAM公司的SFH 4715AS或Vishay公司的VLMU3500-385-120。 2. 供电硬件型号:可以使用常见的直流电源供电,也可以选择专门的LED驱动器,例如Meanwell公司的ELG-75-C或ELG-150-C系列。 3. 广角镜头和探测器型号:一般可采用广角透镜和CMOS摄像头或光电二极管探测器,例如Omron公司的B5W-LA或Murata公司的IRS-B210ST01。 4. 实验设计流程步骤: 1)确定实验目的和研究对象,例如车辆或机器人的定位和导航。
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。
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

实现实时监控告警系统:Kafka与Grafana整合

![实现实时监控告警系统:Kafka与Grafana整合](https://imgconvert.csdnimg.cn/aHR0cHM6Ly9tbWJpei5xcGljLmNuL21tYml6X2pwZy9BVldpY3ladXVDbEZpY1pLWmw2bUVaWXFUcEdLT1VDdkxRSmQxZXB5R1lxaWNlUjA2c0hFek5Qc3FyRktudFF1VDMxQVl3QTRXV2lhSWFRMEFRc0I1cW1ZOGcvNjQw?x-oss-process=image/format,png) # 1.1 Kafka集群架构 Kafka集群由多个称为代理的服务器组成,这
recommend-type

解释这行代码 c = ((double)rand() / RAND_MAX) * (a + b - fabs(a - b)) + fabs(a - b);

这行代码的作用是随机生成一个浮点数,范围在 a 和 b 之间(包括 a 和 b)。 其中,`rand()` 函数是 C 语言标准库中的一个函数,用于生成一个伪随机整数。`RAND_MAX` 是一个常量,它表示 `rand()` 函数生成的随机数的最大值。 因此,`(double)rand() / RAND_MAX` 表示生成的随机数在 [0, 1] 之间的浮点数。 然后,将这个随机数乘上 `(a - b) - fabs(a - b)`,再加上 `fabs(a - b)`。 `fabs(a - b)` 是 C 语言标准库中的一个函数,用于计算一个数的绝对值。因此,`fabs(a - b)
recommend-type

c++校园超市商品信息管理系统课程设计说明书(含源代码) (2).pdf

校园超市商品信息管理系统课程设计旨在帮助学生深入理解程序设计的基础知识,同时锻炼他们的实际操作能力。通过设计和实现一个校园超市商品信息管理系统,学生掌握了如何利用计算机科学与技术知识解决实际问题的能力。在课程设计过程中,学生需要对超市商品和销售员的关系进行有效管理,使系统功能更全面、实用,从而提高用户体验和便利性。 学生在课程设计过程中展现了积极的学习态度和纪律,没有缺勤情况,演示过程流畅且作品具有很强的使用价值。设计报告完整详细,展现了对问题的深入思考和解决能力。在答辩环节中,学生能够自信地回答问题,展示出扎实的专业知识和逻辑思维能力。教师对学生的表现予以肯定,认为学生在课程设计中表现出色,值得称赞。 整个课程设计过程包括平时成绩、报告成绩和演示与答辩成绩三个部分,其中平时表现占比20%,报告成绩占比40%,演示与答辩成绩占比40%。通过这三个部分的综合评定,最终为学生总成绩提供参考。总评分以百分制计算,全面评估学生在课程设计中的各项表现,最终为学生提供综合评价和反馈意见。 通过校园超市商品信息管理系统课程设计,学生不仅提升了对程序设计基础知识的理解与应用能力,同时也增强了团队协作和沟通能力。这一过程旨在培养学生综合运用技术解决问题的能力,为其未来的专业发展打下坚实基础。学生在进行校园超市商品信息管理系统课程设计过程中,不仅获得了理论知识的提升,同时也锻炼了实践能力和创新思维,为其未来的职业发展奠定了坚实基础。 校园超市商品信息管理系统课程设计的目的在于促进学生对程序设计基础知识的深入理解与掌握,同时培养学生解决实际问题的能力。通过对系统功能和用户需求的全面考量,学生设计了一个实用、高效的校园超市商品信息管理系统,为用户提供了更便捷、更高效的管理和使用体验。 综上所述,校园超市商品信息管理系统课程设计是一项旨在提升学生综合能力和实践技能的重要教学活动。通过此次设计,学生不仅深化了对程序设计基础知识的理解,还培养了解决实际问题的能力和团队合作精神。这一过程将为学生未来的专业发展提供坚实基础,使其在实际工作中能够胜任更多挑战。