改写以下代码，使其能够实现“准确率 = （真阳性数 + 真阴性数）/ 总样本数”的功能，并分别对四种故障类型诊断准确率进行计算：% 对测试数据进行测试 [testErr,testY_out] = nntest(nn, testXn, testY); % 转为实际类别 [~,testYreal_out] = max(testY_out'); testYreal_out = testYreal_out'; % 计算准确率 accTest = sum(testYreal_out==testYreal)/length(testYreal);

修改后的代码如下所示：

% 对测试数据进行测试
[testErr,testY_out] = nntest(nn, testXn, testY);

% 转为实际类别
[~,testYreal_out] = max(testY_out');
testYreal_out = testYreal_out';

% 计算真阳性数、真阴性数、总样本数
TP = sum(testYreal_out == 1 & testY == 1);
TN = sum(testYreal_out == 0 & testY == 0);
totalSamples = length(testYreal_out);

% 计算准确率和四种故障类型的诊断准确率
accuracy = (TP + TN) / totalSamples;
accuracyType1 = sum(testYreal_out == 1 & testY == 1) / sum(testY == 1);
accuracyType2 = sum(testYreal_out == 0 & testY == 0) / sum(testY == 0);
accuracyType3 = sum(testYreal_out == 0 & testY == 1) / sum(testY == 1);
accuracyType4 = sum(testYreal_out == 1 & testY == 0) / sum(testY == 0);

其中，`TP` 表示真阳性数，`TN` 表示真阴性数，`totalSamples` 表示总样本数，`accuracyType1`、`accuracyType2`、`accuracyType3`、`accuracyType4` 分别表示四种故障类型的诊断准确率。