ICEEMDAN-iMPA-BiLSTM功率/风速预测 原创全新~直接替换数据即可用 适合新手小白

function [ fitness ] = fitnessFunction(tour,ACOtotalpathloss,ACOtotalenergy,ACOtotalpheromone)% function [ fitness ] = fitnessFunction(ACOedgesenergymatrix,ACOdesirebilitymatrix,tour)fitness = 0;for i = 1 : length(tour)        ACOcurrentroute = tour(i);    %ACOnextnode = tour(i+1);    %     fitness = fitness + ACOedgesenergymatrix(1, ACOcurrentNode);    fitness = fitness + ACOtotalpheromone(ACOcurrentroute, 1);% + (1/ACOtotalenergy(ACOcurrentroute, 1)); % + ACOtotalpathloss(ACOcurrentroute, 1);%     fitness = fitness + ACOdesirebilitymatrix(ACOnextnode, 1);    %fitness = fitness + Gaco.Edges(ACOcurrentNode, ACOnextnode);    end

% Import necessary toolboxesimport phased.*import comm.*% Constant DefinitionsnNodes = 10; % Number of wireless sensor nodescenterFrequency = 5e9; % Operating frequency (5 GHz - suitable for sub-6 GHz 5G)binaryDataSize = 1024; % Size of binary data to be transmitted by each node% Network and Antenna Setupnodes = cell(1,nNodes); % Create a wireless sensor network with antennasfor k = 1:nNodes    % Create an isotropic antenna element for each node, with a specified frequency range of 2 GHz to 8 GHz    nodes{k} = phased.IsotropicAntennaElement('FrequencyRange', [2e9 8e9]);end% Spatial Setup% Random 3D positions for the nodes; Z-coordinates are all 0nodePos = [100*rand(2,nNodes); zeros(1,nNodes)];% Base station is located at the origin in 3DbasePos = [0;0;0];% Define the propagation model with the same operating frequency% Assumes two-way propagation (to and from the base station)propModel = phased.FreeSpace('OperatingFrequency', centerFrequency, 'TwoWayPropagation', true);% Plot Node Positions and Antenna Directionsfigure;hold on;plot(basePos(1), basePos(2), 'ro', 'DisplayName', 'Base Station');for k = 1:nNodes    plot(nodePos(1,k), nodePos(2,k), 'bo', 'DisplayName', sprintf('Node %d', k));    quiver(nodePos(1,k), nodePos(2,k), basePos(1)-nodePos(1,k), basePos(2)-nodePos(2,k), 0.1, 'b');endxlabel('X Coordinate');ylabel('Y Coordinate');title('Node Positions and Antenna Directions');legend('Location','best');hold off;% Preallocate bit errorsbitErrors = zeros(1,nNodes);% Data Transmission and Reception% Iterate over each node, transmit data, and receive it at the base stationfor k = 1:nNodes    % Generate random binary data at each node    data = randi([0 1], binaryDataSize, 1);        % Apply BPSK modulation to the data    modData = pskmod(data, 2);    % Define the carrier signal and modulate the data onto the carrier signal    t = (0:length(modData)-1)';    carrier = cos(2*pi*centerFrequency*t);    txSig = modData .* carrier;    % Compute the angle to the base station from the node    pos_diff = basePos - nodePos(:,k);    azimuth = atan2d(pos_diff(2), pos_diff(1));    % Propagate the signal through the medium    sig = propModel(txSig, nodePos(:,k), basePos, [azimuth;0;0], [0;0;0]);    % Amplify the received signal (simple gain model for the receiving antenna)    receivedSig = sig;    % Demodulate the received signal to recover the original data    recvData = pskdemod(receivedSig, 2);    % Compare the transmitted and received data to check for transmission errors    nErrors = biterr(data, recvData);    % Output the number of errors for each node    fprintf('Node %d: %d errors\n', k, nErrors);    % Store bit errors for each node    bitErrors(k) = nErrors;end% Plot the bit errors for each nodefigure;bar(1:nNodes, bitErrors);xlabel('Node');ylabel('Bit Errors');title('Bit Errors for Each Node');