Matlab实现Logistic映射的递归图分析

需积分: 5 1 下载量 137 浏览量 更新于2024-12-21 收藏 481B ZIP 举报
资源摘要信息:"Logistic映射是一种用来表示生物种群动态模型的数学方程,它是由美国数学家Robert May在1976年提出的。该模型是用来描述在一定环境条件下,种群数量随时间变化的非线性动力系统。Logistic映射方程形式简单,但是具有非常复杂的动态行为,从简单的周期性到混沌状态都有可能出现。 Logistic映射的数学表达式为: x_{n+1} = a \cdot x_n \cdot (1 - x_n) 其中,x_n 表示在第n个时间步的种群数量,取值范围是[0,1],a是一个控制参数,用来表示种群增长率。当a的值变化时,Logistic映射展现出从稳定的不动点到复杂的混沌行为的演化。 在这个问题描述中,给定了a=3,初值x0=0.808,迭代200次,这些参数用于确定Logistic映射的动态行为。嵌入维数为1,延迟时间为1,这两个参数通常用于构建延迟坐标嵌入,它是递归图分析中的一个概念,有助于揭示时间序列中的非线性结构。阈值选择为0.25,这个阈值用于确定递归图中的点之间的连线规则。 递归图是一种用于分析时间序列中数据点之间关系的图示方法,它能够揭示出时间序列中隐藏的动态特性。递归图中的点代表时间序列中的数据点,而点之间的连线代表时间序列数据点之间的相似性或重复性。如果两个点的距离小于某个阈值,则这两个点在图中会用线连接起来。通过观察递归图,研究者可以判断时间序列是否具有周期性、混沌性或是其他复杂的动态特性。 对于提到的源代码文件名'logistic_mapping.m',可以推断这是一个用Matlab语言编写的脚本文件。Matlab是一种用于数值计算、可视化和编程的高级语言和交互式环境。在Matlab中创建递归图,通常需要计算每个数据点与其它数据点的距离,然后根据阈值确定是否连接这些点。问题描述中提到需要在画图时手动取消y反向的问题,这可能是指在绘图时y轴的数值方向与预期相反,需要通过编程调整坐标轴设置来修正。 最后,描述中给出了一个参考文章的链接,可能包含有关Logistic映射和递归图分析的详细说明和实例,读者可以访问该链接获取更多的理论知识和应用背景。"

把这段代码的PCA换成LDA:LR_grid = LogisticRegression(max_iter=1000, random_state=42) LR_grid_search = GridSearchCV(LR_grid, param_grid=param_grid, cv=cvx ,scoring=scoring,n_jobs=10,verbose=0) LR_grid_search.fit(pca_X_train, train_y) estimators = [ ('lr', LR_grid_search.best_estimator_), ('svc', svc_grid_search.best_estimator_), ] clf = StackingClassifier(estimators=estimators, final_estimator=LinearSVC(C=5, random_state=42),n_jobs=10,verbose=1) clf.fit(pca_X_train, train_y) estimators = [ ('lr', LR_grid_search.best_estimator_), ('svc', svc_grid_search.best_estimator_), ] param_grid = {'final_estimator':[LogisticRegression(C=0.00001),LogisticRegression(C=0.0001), LogisticRegression(C=0.001),LogisticRegression(C=0.01), LogisticRegression(C=0.1),LogisticRegression(C=1), LogisticRegression(C=10),LogisticRegression(C=100), LogisticRegression(C=1000)]} Stacking_grid =StackingClassifier(estimators=estimators,) Stacking_grid_search = GridSearchCV(Stacking_grid, param_grid=param_grid, cv=cvx, scoring=scoring,n_jobs=10,verbose=0) Stacking_grid_search.fit(pca_X_train, train_y) Stacking_grid_search.best_estimator_ train_pre_y = cross_val_predict(Stacking_grid_search.best_estimator_, pca_X_train,train_y, cv=cvx) train_res1=get_measures_gridloo(train_y,train_pre_y) test_pre_y = Stacking_grid_search.predict(pca_X_test) test_res1=get_measures_gridloo(test_y,test_pre_y) best_pca_train_aucs.append(train_res1.loc[:,"AUC"]) best_pca_test_aucs.append(test_res1.loc[:,"AUC"]) best_pca_train_scores.append(train_res1) best_pca_test_scores.append(test_res1) train_aucs.append(np.max(best_pca_train_aucs)) test_aucs.append(best_pca_test_aucs[np.argmax(best_pca_train_aucs)].item()) train_scores.append(best_pca_train_scores[np.argmax(best_pca_train_aucs)]) test_scores.append(best_pca_test_scores[np.argmax(best_pca_train_aucs)]) pca_comp.append(n_components[np.argmax(best_pca_train_aucs)]) print("n_components:") print(n_components[np.argmax(best_pca_train_aucs)])

2023-07-22 上传

报错ValueError: np.nan is an invalid document, expected byte or unicode string. 怎么修改import pandas as pd from sklearn.feature_extraction.text import CountVectorizer, TfidfVectorizer from sklearn.model_selection import train_test_split from sklearn.linear_model import LogisticRegression from sklearn.metrics import accuracy_score # 读取电影评论数据集 data = pd.read_csv(r'D:\shujukexue\review_data.csv', encoding='gbk') x = v.fit_transform(df['eview'].apply(lambda x: np.str_(x))) # 分割数据集为训练集和测试集 X_train, X_test, y_train, y_test = train_test_split(data['review'], data['sentiment'], test_size=0.2, random_state=42) # 创建CountVectorizer对象进行词频统计和向量化 count_vectorizer = CountVectorizer() X_train_count = count_vectorizer.fit_transform(X_train) X_test_count = count_vectorizer.transform(X_test) # 创建TfidfVectorizer对象进行TF-IDF计算和向量化 tfidf_vectorizer = TfidfVectorizer() X_train_tfidf = tfidf_vectorizer.fit_transform(X_train) X_test_tfidf = tfidf_vectorizer.transform(X_test) # 创建逻辑回归分类器并在CountVectorizer上进行训练和预测 classifier_count = LogisticRegression() classifier_count.fit(X_train_count, y_train) y_pred_count = classifier_count.predict(X_test_count) accuracy_count = accuracy_score(y_test, y_pred_count) print("Accuracy using CountVectorizer:", accuracy_count) # 创建逻辑回归分类器并在TfidfVectorizer上进行训练和预测 classifier_tfidf = LogisticRegression() classifier_tfidf.fit(X_train_tfidf, y_train) y_pred_tfidf = classifier_tfidf.predict(X_test_tfidf) accuracy_tfidf = accuracy_score(y_test, y_pred_tfidf) print("Accuracy using TfidfVectorizer:", accuracy_tfidf)

2023-05-27 上传

优化这段代码 for j in n_components: estimator = PCA(n_components=j,random_state=42) pca_X_train = estimator.fit_transform(X_standard) pca_X_test = estimator.transform(X_standard_test) cvx = StratifiedKFold(n_splits=5, shuffle=True, random_state=42) cost = [-5, -3, -1, 1, 3, 5, 7, 9, 11, 13, 15] gam = [3, 1, -1, -3, -5, -7, -9, -11, -13, -15] parameters =[{'kernel': ['rbf'], 'C': [2x for x in cost],'gamma':[2x for x in gam]}] svc_grid_search=GridSearchCV(estimator=SVC(random_state=42), param_grid=parameters,cv=cvx,scoring=scoring,verbose=0) svc_grid_search.fit(pca_X_train, train_y) param_grid = {'penalty':['l1', 'l2'], "C":[0.00001,0.0001,0.001, 0.01, 0.1, 1, 10, 100, 1000], "solver":["newton-cg", "lbfgs","liblinear","sag","saga"] # "algorithm":['auto', 'ball_tree', 'kd_tree', 'brute'] } LR_grid = LogisticRegression(max_iter=1000, random_state=42) LR_grid_search = GridSearchCV(LR_grid, param_grid=param_grid, cv=cvx ,scoring=scoring,n_jobs=10,verbose=0) LR_grid_search.fit(pca_X_train, train_y) estimators = [ ('lr', LR_grid_search.best_estimator_), ('svc', svc_grid_search.best_estimator_), ] clf = StackingClassifier(estimators=estimators, final_estimator=LinearSVC(C=5, random_state=42),n_jobs=10,verbose=0) clf.fit(pca_X_train, train_y) estimators = [ ('lr', LR_grid_search.best_estimator_), ('svc', svc_grid_search.best_estimator_), ] param_grid = {'final_estimator':[LogisticRegression(C=0.00001),LogisticRegression(C=0.0001), LogisticRegression(C=0.001),LogisticRegression(C=0.01), LogisticRegression(C=0.1),LogisticRegression(C=1), LogisticRegression(C=10),LogisticRegression(C=100), LogisticRegression(C=1000)]} Stacking_grid =StackingClassifier(estimators=estimators,) Stacking_grid_search = GridSearchCV(Stacking_grid, param_grid=param_grid, cv=cvx, scoring=scoring,n_jobs=10,verbose=0) Stacking_grid_search.fit(pca_X_train, train_y) var = Stacking_grid_search.best_estimator_ train_pre_y = cross_val_predict(Stacking_grid_search.best_estimator_, pca_X_train,train_y, cv=cvx) train_res1=get_measures_gridloo(train_y,train_pre_y) test_pre_y = Stacking_grid_search.predict(pca_X_test) test_res1=get_measures_gridloo(test_y,test_pre_y) best_pca_train_aucs.append(train_res1.loc[:,"AUC"]) best_pca_test_aucs.append(test_res1.loc[:,"AUC"]) best_pca_train_scores.append(train_res1) best_pca_test_scores.append(test_res1) train_aucs.append(np.max(best_pca_train_aucs)) test_aucs.append(best_pca_test_aucs[np.argmax(best_pca_train_aucs)].item()) train_scores.append(best_pca_train_scores[np.argmax(best_pca_train_aucs)]) test_scores.append(best_pca_test_scores[np.argmax(best_pca_train_aucs)]) pca_comp.append(n_components[np.argmax(best_pca_train_aucs)]) print("n_components:") print(n_components[np.argmax(best_pca_train_aucs)])

2023-07-14 上传