X = np.array(X), X = np.real(X),arr_real = arr.astype('float'),这样可以吗
时间: 2024-02-26 12:58:34 浏览: 112
如果你想将一个复数类型的数组转换成实数类型的数组,同时保留实数部分,你可以使用`np.real()`函数,像这样:
```python
X = np.array(X)
X = np.real(X)
```
这将把`X`数组转换为实数类型的数组,并保留实数部分。
如果你想将一个复数类型的数组转换成实数类型的数组,并舍弃虚数部分,你可以使用`astype()`方法,像这样:
```python
arr_real = arr.astype('float')
```
这将把`arr`数组转换为实数类型的数组,并舍弃虚数部分。
需要注意的是,如果你将一个复数类型的数组转换为实数类型的数组,并舍弃虚数部分,你将会丢失一些信息。因此,如果你需要保留复数部分,你应该使用`np.real()`函数,而不是使用`astype()`方法。
相关问题
import os import numpy as np import matplotlib.pyplot as plt from PIL import Image from colorcet.plotting import arr from sklearn.cluster import SpectralClustering from sklearn.decomposition import PCA from tensorflow.keras.preprocessing import image from tensorflow.keras.applications.resnet50 import ResNet50 from tensorflow.keras.applications.resnet50 import preprocess_input # 定义加载图片函数 def load_image(img_path): img = image.load_img(img_path, target_size=(224, 224)) x = image.img_to_array(img) x = np.expand_dims(x, axis=0) x = preprocess_input(x) return x # 加载ResNet50模型 model = ResNet50(weights='imagenet', include_top=False, pooling='avg') # 加载图片并提取特征向量 img_dir = 'D:/wjd' img_names = os.listdir(img_dir) X = [] for img_name in img_names: img_path = os.path.join(img_dir, img_name) img = load_image(img_path) features = model.predict(img)[0] X.append(features) # 将特征向量转化为矩阵 X = np.array(X) X = np.real(X) arr_real = arr.astype('float') # 计算相似度矩阵 S = np.dot(X, X.T) # 归一化相似度矩阵 D = np.diag(np.sum(S, axis=1)) L = D - S L_norm = np.dot(np.dot(np.sqrt(np.linalg.inv(D)), L), np.sqrt(np.linalg.inv(D))) # 计算特征向量 eigvals, eigvecs = np.linalg.eig(L_norm) idx = eigvals.argsort()[::-1] eigvals = eigvals[idx] eigvecs = eigvecs[:, idx] Y = eigvecs[:, :2] # 使用谱聚类进行分类 n_clusters = 5 clustering = SpectralClustering(n_clusters=n_clusters, assign_labels="discretize", random_state=0).fit(Y) # 可视化聚类结果 pca = PCA(n_components=2) X_pca = pca.fit_transform(X) plt.scatter(X_pca[:, 0], X_pca[:, 1], c=clustering.labels_, cmap='rainbow') plt.show(),这行代码出现了这个numpy.ComplexWarning: Casting complex values to real discards the imaginary part The above exception was the direct cause of the following exception问题
在代码中,出现了 `X = np.real(X)` 这一行代码,它会将复数类型的特征向量转换为实数类型,导致了警告信息 `numpy.ComplexWarning: Casting complex values to real discards the imaginary part`。这个警告信息是因为转换复数类型的数据时丢失了虚部信息,可能会影响后续的计算结果。
如果你需要使用实数类型的特征向量进行后续的计算,可以使用 `numpy.absolute` 函数将复数类型的数据转换为实数类型,而不丢失虚部的信息。例如,将 `X` 转换为实数类型可以使用以下代码:
```python
X = np.absolute(X)
```
这段代码会将 `X` 中的复数数值转换为实数数值,并保留虚部的信息。如果后续计算需要使用实数类型的数据,可以使用转换后的 `X`,而不会丢失虚部的信息。
arr0 = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]) arr1 = np.array([1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24]) arr3 = np.array(input("请输入连续24个月的配件销售数据,元素之间用空格隔开:").split(), dtype=float) data_array = np.vstack((arr1, arr3)) data_matrix = data_array.T data = pd.DataFrame(data_matrix, columns=['month', 'sales']) sales = data['sales'].values.astype(np.float32) sales_mean = sales.mean() sales_std = sales.std() sales = abs(sales - sales_mean) / sales_std train_data = sales[:-1] test_data = sales[-12:] def create_model(): model = tf.keras.Sequential() model.add(layers.Input(shape=(11, 1))) model.add(layers.Conv1D(filters=32, kernel_size=2, padding='causal', activation='relu')) model.add(layers.BatchNormalization()) model.add(layers.Conv1D(filters=64, kernel_size=2, padding='causal', activation='relu')) model.add(layers.BatchNormalization()) model.add(layers.Conv1D(filters=128, kernel_size=2, padding='causal', activation='relu')) model.add(layers.BatchNormalization()) model.add(layers.Conv1D(filters=256, kernel_size=2, padding='causal', activation='relu')) model.add(layers.BatchNormalization()) model.add(layers.Conv1D(filters=512, kernel_size=2, padding='causal', activation='relu')) model.add(layers.BatchNormalization()) model.add(layers.Dense(1, activation='linear')) return model model = create_model() BATCH_SIZE = 16 BUFFER_SIZE = 100 train_dataset = tf.data.Dataset.from_tensor_slices(train_data) train_dataset = train_dataset.window(11, shift=1, drop_remainder=True) train_dataset = train_dataset.flat_map(lambda window: window.batch(11)) train_dataset = train_dataset.map(lambda window: (window[:-1], window[-1:])) train_dataset = train_dataset.shuffle(BUFFER_SIZE).batch(BATCH_SIZE).prefetch(1) model.compile(optimizer=tf.keras.optimizers.Adam(learning_rate=0.001), loss='mse') history = model.fit(train_dataset, epochs=100, verbose=0) test_input = test_data[:-1] test_input = np.reshape(test_input, (1, 11, 1)) predicted_sales = model.predict(test_input)[0][0] * sales_std + sales_mean test_prediction = model.predict(test_input) y_test=test_data[1:12] y_pred=test_prediction y_pred = test_prediction.ravel() print("预测下一个月的销量为:", predicted_sales),如何将以下代码稍作修改插入到上面的最后,def comput_acc(real,predict,level): num_error=0 for i in range(len(real)): if abs(real[i]-predict[i])/real[i]>level: num_error+=1 return 1-num_error/len(real) a=np.array(test_data[label]) real_y=a real_predict=test_predict print("置信水平:{},预测准确率:{}".format(0.2,round(comput_acc(real_y,real_predict,0.2)* 100,2)),"%")
你可以在最后添加如下代码实现 `comput_acc()` 函数的功能:
```
label = 0 # 定义标签
a = np.array(test_data[label]) # 获取测试数据的标签列
real_y = a[1:] # 实际销售数据
real_predict = y_pred[:-1] # 预测销售数据
print("置信水平:{},预测准确率:{}".format(0.2, round(comput_acc(real_y, real_predict, 0.2) * 100, 2)), "%")
```
这段代码的作用是,首先定义标签为 0,即选取测试数据中的第一列数据(即销售数据)。然后,使用 `np.array()` 函数将该列数据转换为 numpy 数组 `a`。接着,将 `a` 中的第二个元素到最后一个元素赋值给 `real_y`,这里是因为真实销售数据的第一个元素已经作为测试数据输入了模型,所以预测值中不包含该元素。然后,将 `y_pred` 中的第一个元素到倒数第二个元素赋值给 `real_predict`,这是因为预测值中的最后一个元素已经与真实值的最后一个元素相对应。最后,调用 `comput_acc()` 函数计算预测准确率,并将结果打印输出。
阅读全文