A secure aggregation method and system for power time series data communication

Abstract

The present invention relates to the technical field of secure data transmission, and discloses a secure aggregation method and system for power time series data communication. The method comprises: preprocessing the collected power time series data, decomposing the preprocessed power time series data, fitting similar power time series data by Gaussian fitting, ensuring secure transmission of power time series data by digital modulation, and finally realizing aggregation of power time series data by carrier aggregation, thereby reducing redundancy in the process of power time series data transmission and improving the security and efficiency of power time series data communication.

Description

Technical Field

[0001] The present invention relates to the technical field of secure data transmission, and in particular to a secure aggregation method and system for power time series data communication. Background Art

[0002] In the past decade, the total electricity consumption in the society has been increasing day by day. With the construction of smart grids, smart meters have entered thousands of households, from mechanical meters in the last century to today's IoT smart meters. The popularity of these smart meters has provided massive load data for the power system, but at the same time, it has also brought huge pressure to the safe transmission of load data.

[0003] Chinese patent CN115834259B discloses a method for secure transmission of data over a high-risk communication line. The method only obtains the data to be transmitted as the main body data; uses the data characteristics of the data to be transmitted to copy the data to obtain the clone data; and uses the main body data and the clone data to transmit data simultaneously to complete the secure transmission process, which increases the cost of secure transmission and has great limitations. Summary of the invention

[0004] 1. Technical issues to be resolved

[0005] In view of the deficiencies in the prior art, the present invention provides a secure aggregation method and system for power timing data communication, which has the advantages of accuracy and security, and solves the problem of secure transmission of load data.

[0006] (II) Technical solution

[0007] In order to solve the above technical problem of secure transmission of load data, the present invention provides the following technical solution:

[0008] The present invention discloses a secure aggregation method for power time series data communication, which specifically comprises the following steps:

[0009] S1. Collect historical power time series data to establish a power time series data sample set, and preprocess the power time series data in the power time series data sample set to obtain preprocessed power time series data, wherein the power time series data in the sample set is a collection of load data;

[0010] S2. Decomposing the preprocessed power time series data based on a time series decomposition algorithm to obtain power time series data after time series decomposition;

[0011] S3, fitting the power time series data after time series decomposition to obtain a fitted and transmittable power time series signal, and modulating the fitted power time series signal by digital modulation to obtain a modulated power time series signal;

[0012] S4, aggregating the modulated power timing signal through carrier aggregation technology;

[0013] The present invention realizes the secure aggregation of power time series data communication by preprocessing the collected power time series data and performing decomposition, fitting modulation and aggregation operations, thereby reducing redundancy in the power time series data transmission process and improving the security and efficiency of power time series data communication.

[0014] Preferably, the preprocessing of the power time series data in the power time series data sample set to obtain the preprocessed power time series data comprises the following steps:

[0015] S11. Missing and supplementing load data;

[0016] Setting a load data missing interval threshold and supplementing different load data according to the load data missing interval;

[0017] When the load data missing interval is 1 day, the average load data value of that month in the sample set is calculated, and the average load data value of that month is used as the supplementary data for the current missing load data;

[0018] When the load data missing interval is 1 hour, the average load data value of that day in the sample set is calculated, and the average load data value of that day is used as the supplementary data for the current missing load data;

[0019] When the load data missing interval is 1 minute, the average load data value of that hour in the sample set is calculated, and the average load data value of that hour is used as the supplementary data for the currently missing load data;

[0020] The calculation formula for the average load data value of the sample is as follows:

[0021]

[0022] Where n is the number of samples, y k is the kth missing load data, is the average load data value, y i is the i-th load data in the sample set;

[0023] S12, normalization processing of load data;

[0024] The min-max normalization method is used to normalize all load data in the sample set to the range of [0,1];

[0025] The calculation formula of the min-max normalization method is as follows:

[0026]

[0027] in, is the new value of the i-th load data in the normalized sample set, y i is the i-th load data value in the sample set before normalization, y min is the minimum value of the load data in the sample set, y max is the maximum value of the load data in the sample set.

[0028] The present invention sets a load data missing interval threshold and supplements different load data according to the load data missing interval, and at the same time implements preprocessing of the power time series data by normalization processing, thereby improving the accuracy of the power time series data.

[0029] Preferably, the decomposing the preprocessed power time series data based on the time series decomposition algorithm to obtain the power time series data after time series decomposition comprises the following steps:

[0030] Decomposing the preprocessed power time series data by Fourier decomposition;

[0031] The Fourier series expansion is as follows:

[0032]

[0033] Where N is the total number of load data, α0 is the DC component, Y(t) is a sequence of load data about time t, ω j is the angular frequency of the jth load data, a j is the coefficient of the jth load data, b j is the coefficient of the jth load data;

[0034] Based on the angular frequency components and the periodic variation law, Y(t) is decomposed and reconstructed. The decomposed and reconstructed Y(t) is as follows:

[0035] Y(t)=α0+D(t)+W(t)+F a (t)+L0(t);

[0036] Where D(t) is the daily component of a series of load data at time t; W(t) is the weekly component of a series of load data at time t, and F a L0(t) is the high-frequency component of a segment of the load data sequence at time t, and L0(t) is the low-frequency component of a segment of the load data sequence at time t.

[0037] The present invention decomposes the preprocessed power time series data by using Fourier decomposition, thereby removing redundant data in the power time series data and improving the efficiency of power time series data transmission.

[0038] Preferably, fitting the power time series data after time series decomposition to obtain a fitted transmittable power time series signal comprises the following steps:

[0039] Based on the most recent ten sets of power time series data collected, the selected data are fitted by Gaussian fitting;

[0040] The Gaussian fitting expression is as follows:

[0041]

[0042] Where m represents the order of Gaussian fitting, d l Indicates the height of the power time series data curve l, b l represents the coordinate of the center position of the power time series data curve l on the x-axis, c l represents the width of the curve l, F(x l ) is the Gaussian fitting function, x l is the coordinate value of the power timing data curve l in the x-axis.

[0043] The present invention fits multiple groups of power time series data into a group of transmittable power time series signals by using Gaussian fitting, thereby reducing redundant data existing in the multiple groups of power time series data and ensuring the efficiency of the power time series data transmission process.

[0044] Preferably, the modulating the fitted power timing signal by digital modulation comprises the following steps:

[0045] S31, amplitude keying modulation;

[0046] S32, phase shift keying modulation;

[0047] S33, frequency shift keying modulation.

[0048] Preferably, the amplitude shift keying modulation comprises the following steps:

[0049] Modulation is performed by adjusting the amplitude value of the fitted transmittable power timing signal;

[0050] The expression of setting power timing signal is:

[0051]

[0052] Among them, g represents the matrix pulse, R s Indicates the pulse width of the power timing signal s, a η is the level value of the nth load, the level value is set to be represented as 0 or 1 in the binary signal, t represents time, and s(n) represents the power timing signal of n loads;

[0053] The time domain expression of amplitude keying modulation is:

[0054]

[0055] Among them, the pulse width of the power timing signal s is R s , S MASK (t) is the time domain of amplitude keying modulation, g represents the matrix pulse, a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, and t represents the time;

[0056] The amplitude keying modulation time domain is summed to obtain the amplitude modulated power timing signal

[0057] Preferably, the phase shift keying modulation comprises the following steps:

[0058] Phase modulating the amplitude modulated power time series signal through the fitted transmittable power time series signal;

[0059] The time domain expression of phase shift keying modulation is set as:

[0060]

[0061] Among them, S MPSK (t) is the time domain of phase shift keying modulation, g represents the matrix pulse, and the pulse width of the power timing signal s is R s , a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, t represents the time, It is represented as the phase corresponding to the nth load;

[0062] in:

[0063] The phase-shift keying modulated time domain is summed to obtain the phase modulated power timing signal

[0064] Preferably, the frequency shift keying modulation comprises the following steps:

[0065] The modulated power timing signal is modulated to a different frequency by frequency shift keying modulation, keeping the amplitude and phase unchanged;

[0066] The time domain expression of frequency shift keying modulation is set to:

[0067]

[0068] Among them, S MFSK (t) is the time domain of frequency shift keying modulation, g represents the matrix pulse, and the pulse width of the power timing signal s is Rs , a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power time series signal, t represents time, Δf c is the carrier frequency offset;

[0069] The frequency shift keying modulated time domain is summed to obtain the power timing signal after frequency shift keying modulation

[0070]

[0071] The power timing signal after frequency shift keying modulation is set to be the modulated power timing signal.

[0072] The present invention modulates the fitted power timing signal by using digital modulation, modulates the power timing signal to be transmitted into a safe, reliable and transmittable power timing signal, and ensures the security and accuracy of power timing data transmission.

[0073] Preferably, aggregating the modulated power timing signal by using carrier aggregation technology comprises the following steps:

[0074] The modulated power timing signal is set to be transmitted through the carrier unit;

[0075] Multiple narrow-bandwidth carrier units are aggregated in the frequency domain through the DFT / IDFT algorithm and modulated onto an optical carrier for transmission.

[0076] According to the principle of code division multiple access, each bit is set to be sent as a set of code chip sequences;

[0077] The chip is a component unit of a data signal. A data signal is composed of multiple coded signals, one of which is a chip.

[0078] When representing a chip sequence, 1 is used to represent binary 1, and -1 is used to represent binary 0;

[0079] and limiting the order of the chip sequence based on the number of carrier components in the transmission process;

[0080] Assume that a Z-order chip sequence is allocated to M carrier units, then Z must meet the following conditions:

[0081] Z ≥ 2M;

[0082] When the chip sequence is of order 4, it means that there are four sub-functions in the chip sequence, namely:

[0083] {wa o (x);x∈(0,T),o=0,1,2,3};

[0084] Among them, wa o (x) is the relationship function of the value of the chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, and o represents the number of sub-functions;

[0085] The following relationship is satisfied between each sequence of the chip sequence:

[0086]

[0087] Among them, wa o (x) is the relationship function of the value of the code chip sequence of the modulated power timing signal on the horizontal axis x, wa p (x) is a function that limits the value of the chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, o represents the number of sub-functions, and p represents the limiting relationship about o;

[0088] After allocating different chip sequences, the spectrum of each signal is spread to keep different signals independent and orthogonal in the code domain;

[0089] The formula for spreading the spectrum of each signal is as follows:

[0090]

[0091] w (h) =[wa h,1 ,wa h,2 ,...,wa h,n ];

[0092] Among them, q h represents the spread spectrum signal for the hth channel, represents the nth feature of the hth signal, wa (h) is the chip sequence assigned to the hth signal, T is the period, wa h,n represents the nth chip in the chip sequence;

[0093] The spread spectrum formula is as follows:

[0094] [x n (k) w k,1 ,x n (k) w k,2 ,...,x n (k) w k,n ];

[0095] After the signals from different channels are spread n times, the aggregate signal Q generated by adding them is:

[0096]

[0097] Wherein, Q is the aggregate signal, and M represents the number of carrier units.

[0098] The present invention aggregates the modulated power timing signal by using carrier aggregation technology, and completes spectrum aggregation of multiple narrow-bandwidth carrier units in the frequency domain through the DFT / IDFT algorithm, and modulates them to the optical carrier for transmission, thereby reducing the number of data transmission rounds and improving the efficiency of power timing data transmission.

[0099] The present invention also discloses a safety aggregation system for power time series data communication, which specifically includes: a data acquisition module, a database, a data processing module, a data communication decomposition module, a data communication fitting module, a data communication modulation module and a data communication aggregation module;

[0100] The data acquisition module is used to collect historical power time series data to establish a power time series data sample set and to collect power time series data in real time;

[0101] The database is used to store the power time series data collected by the data acquisition module;

[0102] The data processing module is used to pre-process the collected power time series data;

[0103] The data communication decomposition module is used to decompose the preprocessed power time series data;

[0104] The data communication fitting module is used to fit the decomposed power time series data;

[0105] The data communication modulation module is used to modulate the fitted power time series data;

[0106] The data communication aggregation module is used to aggregate the modulated power time series data.

[0107] (III) Beneficial effects

[0108] Compared with the prior art, the present invention provides a secure aggregation method and system for power time series data communication, which has the following beneficial effects:

[0109] 1. The invention realizes the secure aggregation of power time series data communication by preprocessing the collected power time series data and performing decomposition, fitting modulation and aggregation operations, thereby reducing the redundancy in the power time series data transmission process and improving the security and efficiency of power time series data communication.

[0110] 2. The present invention decomposes the pre-processed power time series data by using Fourier decomposition, thereby removing redundant data in the power time series data and improving the efficiency of power time series data transmission.

[0111] 3. The present invention uses Gaussian fitting to fit multiple sets of power time series data into a set of transmittable power time series signals, thereby reducing redundant data in the multiple sets of power time series data and ensuring the efficiency of the power time series data transmission process.

[0112] 4. The invention uses digital modulation to modulate the fitted power timing signal, modulating the power timing signal to be transmitted into a safe, reliable and transmittable power timing signal, thereby ensuring the security and accuracy of power timing data transmission.

[0113] 5. This invention aggregates the modulated power timing signal by using carrier aggregation technology, and completes spectrum aggregation of multiple narrow bandwidth carrier units in the frequency domain through the DFT / IDFT algorithm, and modulates them to the optical carrier for transmission, thereby reducing the number of data transmission rounds and improving the efficiency of power timing data transmission. BRIEF DESCRIPTION OF THE DRAWINGS

[0114] Figure 1 It is a schematic diagram of the structure of the process of the secure aggregation method of power timing data communication of the present invention. DETAILED DESCRIPTION

[0115] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0116] The present invention discloses a secure aggregation method for power time series data communication, which specifically comprises the following steps:

[0117] S1. Collect historical power time series data to establish a power time series data sample set, and preprocess the power time series data in the power time series data sample set to obtain preprocessed power time series data;

[0118] The power time series data in the sample set is a collection of load data;

[0119] S11. Missing and supplementing load data;

[0120] Setting a load data missing interval threshold and supplementing different load data according to the load data missing interval;

[0121] When the load data missing interval is 1 day, the average load data value of that month in the sample set is calculated, and the average load data value of that month is used as the supplementary data for the current missing load data;

[0122] When the load data missing interval is 1 hour, the average load data value of that day in the sample set is calculated, and the average load data value of that day is used as the supplementary data for the current missing load data;

[0123] When the load data missing interval is 1 minute, the average load data value of that hour in the sample set is calculated, and the average load data value of that hour is used as the supplementary data for the currently missing load data;

[0124] The calculation formula for the average load data value of the sample is as follows:

[0125]

[0126] Where n is the number of samples, y k is the kth missing load data, is the average load data value, y i is the i-th load data in the sample set;

[0127] S12, normalization processing of load data;

[0128] The min-max normalization method is used to normalize all load data in the sample set to the range of [0,1];

[0129] The calculation formula of the min-max normalization method is as follows:

[0130]

[0131] in, is the new value of the i-th load data in the normalized sample set, y i is the i-th load data value in the sample set before normalization, y min is the minimum value of the load data in the sample set, y max is the maximum value of the load data in the sample set;

[0132] S2. Decomposing the preprocessed power time series data based on a time series decomposition algorithm to obtain power time series data after time series decomposition;

[0133] Decomposing the preprocessed power time series data by Fourier decomposition;

[0134] The Fourier series expansion is as follows:

[0135]

[0136] Where N is the total number of load data, α0 is the DC component, Y(t) is a sequence of load data about time t, ω j is the angular frequency of the jth load data, a jis the coefficient of the jth load data, b j is the coefficient of the jth load data;

[0137] Furthermore, based on the angular frequency component and the periodic variation law, Y(t) is decomposed and reconstructed. The decomposed and reconstructed Y(t) is as follows:

[0138] Y(t)=α0+D(t)+W(t)+F a (t)+L0(t);

[0139] Where D(t) is the daily component of a series of load data at time t; W(t) is the weekly component of a series of load data at time t, and F a L0(t) is the high frequency component of a sequence of load data at time t, and L0(t) is the low frequency component of a sequence of load data at time t;

[0140] S3, fitting the power time series data after time series decomposition to obtain a fitted and transmittable power time series signal, and modulating the fitted power time series signal by digital modulation to obtain a modulated power time series signal;

[0141] Based on the most recent ten sets of power time series data collected, the selected data are fitted by Gaussian fitting;

[0142] The Gaussian fitting expression is as follows:

[0143]

[0144] Where m represents the order of Gaussian fitting, d l Indicates the height of the power time series data curve l, b l represents the coordinate of the center position of the power time series data curve l on the x-axis, c l represents the width of the curve l, F(x l ) is the Gaussian fitting function, x l is the coordinate value of the power time series data curve l in the x-axis;

[0145] Modulating the power timing signal by digital modulation includes the following steps:

[0146] S31, amplitude keying modulation;

[0147] Modulation is performed by adjusting the amplitude value of the fitted transmittable power timing signal;

[0148] The expression of setting power timing signal is:

[0149]

[0150] Where g represents the matrix pulse, R s Indicates the pulse width of the power timing signal s, a η is the level value of the nth load, the level value is set to be represented as 0 or 1 in the binary signal, t represents time, and s(n) represents the power timing signal of n loads;

[0151] The time domain expression of amplitude keying modulation is:

[0152]

[0153] Among them, the pulse width of the power timing signal s is R s , S MASK is the time domain of amplitude keying modulation, g represents the matrix pulse, a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, and t represents the time;

[0154] The amplitude keying modulation time domain is summed to obtain the amplitude modulated power timing signal

[0155] S32, phase shift keying modulation;

[0156] Phase modulating the amplitude modulated power time series signal through the fitted transmittable power time series signal;

[0157] The time domain expression of phase shift keying modulation is set as:

[0158]

[0159] Among them, S MPSK (t) is the time domain of phase shift keying modulation, g represents the matrix pulse, and the pulse width of the power timing signal s is R s , a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, t represents the time, It is represented as the phase corresponding to the nth load;

[0160] in:

[0161] The phase-shift keying modulated time domain is summed to obtain the phase modulated power timing signal

[0162] S33, frequency shift keying modulation;

[0163] The modulated power timing signal is modulated to a different frequency by frequency shift keying modulation, keeping the amplitude and phase unchanged;

[0164] The time domain expression of frequency shift keying modulation is set to:

[0165]

[0166] Among them, S MFSK (t) is the time domain of frequency shift keying modulation, g represents the matrix pulse, and the pulse width of the power timing signal s is R s , a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power time series signal, t represents time, Δf c is the carrier frequency offset;

[0167] The frequency shift keying modulated time domain is summed to obtain the power timing signal after frequency shift keying modulation

[0168]

[0169] The power timing signal after frequency shift keying modulation is set to be the modulated power timing signal.

[0170] S4, aggregating the modulated power timing signal through carrier aggregation technology;

[0171] The modulated power timing signal is set to be transmitted through the carrier unit;

[0172] Multiple narrow-bandwidth carrier units are aggregated in the frequency domain through the DFT / IDFT algorithm and modulated onto an optical carrier for transmission.

[0173] According to the principle of code division multiple access, each bit is set to be sent as a set of code chip sequences;

[0174] The chip is a component unit of a data signal. A data signal is composed of multiple coded signals, one of which is a chip.

[0175] Furthermore, when representing the chip sequence, 1 is used to represent binary 1, and -1 is used to represent binary 0;

[0176] Further, the order of the chip sequence is limited based on the number of carrier components in the transmission process;

[0177] Assume that a Z-order chip sequence is allocated to M carrier units, then Z must meet the following conditions:

[0178] Z ≥ 2M;

[0179] Furthermore, when the chip sequence is of order 4, it means that there are four sub-functions in the chip sequence, namely:

[0180] {wa o(x);x∈(0,T),o=0,1,2,3};

[0181] Among them, wa o (x) is the relationship function of the value of the chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, and o represents the number of sub-functions;

[0182] Furthermore, the following relationship is satisfied between the various sequences of the chip sequence:

[0183]

[0184] Among them, wa o (x) is the relationship function of the value of the code chip sequence of the modulated power timing signal on the horizontal axis x, wa p (x) is a function that limits the value of the chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, o represents the number of sub-functions, p represents the limit relationship about o, and C represents a constant;

[0185] Furthermore, after allocating different chip sequences, the spectrum of each signal is spread to keep different signals independent and orthogonal in the code domain;

[0186] The formula for spreading the spectrum of each signal is as follows:

[0187]

[0188] w (h) =[wa h,1 ,wa h,2 ,...,wa h,n ];

[0189] Among them, q h represents the spread spectrum signal for the hth channel, represents the nth feature of the hth signal, wa (h) is the chip sequence assigned to the hth signal, T is the period, wa h,n represents the nth chip in the chip sequence;

[0190] The spread spectrum formula is as follows:

[0191] [x n (k) w k,1 ,x n (k) w k,2 ,...,x n (k) w k,n ];

[0192] Furthermore, after the signals of different paths are spread n times, the aggregate signal Q generated by adding them is:

[0193]

[0194] Wherein, Q is the aggregate signal, and M represents the number of carrier units;

[0195] The present invention also discloses a safety aggregation system for power time series data communication, which specifically includes: a data acquisition module, a database, a data processing module, a data communication decomposition module, a data communication fitting module, a data communication modulation module and a data communication aggregation module;

[0196] The data acquisition module is used to collect historical power time series data to establish a power time series data sample set and to collect power time series data in real time;

[0197] The database is used to store the power time series data collected by the data acquisition module;

[0198] The data processing module is used to pre-process the collected power time series data;

[0199] The data communication decomposition module is used to decompose the preprocessed power time series data;

[0200] The data communication fitting module is used to fit the decomposed power time series data;

[0201] The data communication modulation module is used to modulate the fitted power time series data;

[0202] The data communication aggregation module is used to aggregate the modulated power time series data.

[0203] Although embodiments of the present invention have been shown and described, it will be appreciated by those skilled in the art that various changes, modifications, substitutions and variations may be made to the embodiments without departing from the principles and spirit of the present invention, and that the scope of the present invention is defined by the appended claims and their equivalents.

Claims

1. A secure aggregation method for power time series data communication, characterized in that: The following steps are involved: S1. Collect historical power time series data to establish a power time series data sample set, and preprocess the power time series data in the power time series data sample set to obtain preprocessed power time series data, wherein the power time series data in the sample set is a collection of load data; The preprocessing of the power time series data in the power time series data sample set to obtain the preprocessed power time series data comprises the following steps: S11. Missing and supplementing load data; Setting a load data missing interval threshold and supplementing different load data according to the load data missing interval; When the load data missing interval is 1 day, the average load data value of that month in the sample set is calculated, and the average load data value of that month is used as the supplementary data for the current missing load data; When the load data missing interval is 1 hour, the average load data value of that day in the sample set is calculated, and the average load data value of that day is used as the supplementary data for the current missing load data; When the load data missing interval is 1 minute, the average load data value of that hour in the sample set is calculated, and the average load data value of that hour is used as the supplementary data for the currently missing load data; The calculation formula for the average load data value of the sample is as follows: Where n is the number of samples, y k is the kth missing load data, is the average load data value, y i is the i-th load data in the sample set; S12, normalization processing of load data; The min-max normalization method is used to normalize all load data in the sample set to the range of [0,1]; The calculation formula of the min-max normalization method is as follows: in, is the new value of the i-th load data in the normalized sample set, y i is the i-th load data value in the sample set before normalization, y min is the minimum value of the load data in the sample set, y max is the maximum value of the load data in the sample set; S2. Decomposing the preprocessed power time series data based on a time series decomposition algorithm to obtain power time series data after time series decomposition; Decomposing the preprocessed power time series data based on the time series decomposition algorithm to obtain the power time series data after time series decomposition includes the following steps: Decomposing the preprocessed power time series data by Fourier decomposition; The Fourier series expansion is as follows: Where N is the total number of load data, α0 is the DC component, Y(t) is a sequence of load data about time t, ω j is the angular frequency of the jth load data, a j is the coefficient of the jth load data, b j is the coefficient of the jth load data; Based on the angular frequency components and the periodic variation law, Y(t) is decomposed and reconstructed. The decomposed and reconstructed Y(t) is as follows: Y(t)=α0+D(t)+W(t)+F a (t)+L0(t); Where D(t) is the daily component of a series of load data at time t; W(t) is the periodic component of a series of load data at time t, and F a L0(t) is the high frequency component of a sequence of load data at time t, and L0(t) is the low frequency component of a sequence of load data at time t; S3, fitting the power time series data after time series decomposition to obtain a fitted and transmittable power time series signal, and modulating the fitted power time series signal by digital modulation to obtain a modulated power time series signal; The step of fitting the power time series data after time series decomposition to obtain a fitted transmittable power time series signal comprises the following steps: Based on the most recent ten sets of power time series data collected, the selected data are fitted by Gaussian fitting; The Gaussian fitting expression is as follows: Where m represents the order of Gaussian fitting, d l Indicates the height of the power time series data curve l, b l represents the coordinate of the center position of the power time series data curve l on the x-axis, c l represents the width of the curve l, F(x l ) is the Gaussian fitting function, x l is the coordinate value of the power time series data curve l in the x-axis; S4, aggregating the modulated power timing signal through carrier aggregation technology; The aggregating the modulated power timing signal by using the carrier aggregation technology includes the following steps: The modulated power timing signal is set to be transmitted through the carrier unit; Multiple narrow-bandwidth carrier units are aggregated in the frequency domain through the DFT / IDFT algorithm and modulated onto an optical carrier for transmission. According to the principle of code division multiple access, each bit is set to be sent as a set of code chip sequences; The chip is a component unit of a data signal. A data signal is composed of multiple coded signals, one of which is a chip. When representing a chip sequence, 1 is used to represent binary 1, and -1 is used to represent binary 0; and limiting the order of the chip sequence based on the number of carrier components during transmission; Assume that a Z-order chip sequence is allocated to M carrier units, then Z must meet the following conditions: Z ≥ 2M; When the chip sequence is of order 4, it means that there are four sub-functions in the chip sequence, namely: {wa o (x);x∈(0,T),o=0,1,2,3}; Among them, wa o (x) is the relationship function of the value of the chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, and o represents the number of sub-functions; The following relationship is satisfied between each sequence of the chip sequence: Among them, wa o (x) is the relationship function of the code chip sequence of the modulated power timing signal on the horizontal axis x, T represents the period, wa p (x) is a function that limits the value of the chip sequence of the modulated power time series signal on the horizontal axis x, o represents the number of sub-functions, p represents the limiting relationship about o, and C represents a constant; After allocating different chip sequences, the spectrum of each signal is spread to keep different signals independent and orthogonal in the code domain; The formula for spreading the spectrum of each signal is as follows: wa (h) =[wa h,1 ,wa h,2 ,...,wa h,n ]; Among them, q h represents the spread spectrum signal for the hth channel, represents the nth feature of the hth signal, wa (h) is the chip sequence assigned to the hth signal, T is the period, wa h,n represents the nth chip in the chip sequence; The spread spectrum formula is as follows: [x n (k) wa k,1 ,x n (k) wa k,2 ,...,x n (k) wa k,n ]; After the signals from different channels are spread n times, the aggregate signal Q generated by adding them is: Wherein, Q is the aggregate signal, and M represents the number of carrier units.

2. A secure aggregation method for power time series data communication according to claim 1, characterized in that: The method of modulating the fitted power timing signal by digital modulation includes the following steps: S31, amplitude keying modulation; S32, phase shift keying modulation; S33, frequency shift keying modulation.

3. A secure aggregation method for power time series data communication according to claim 2, characterized in that: The amplitude keying modulation comprises the following steps: The amplitude value of the fitted transmittable power time series signal is adjusted to perform modulation, thereby obtaining an amplitude modulated power time series signal; The expression of setting power timing signal is: Where, g represents the matrix pulse, Rs represents the pulse width of the power timing signal s, and a η is the level value of the nth load, the level value is set to be represented as 0 or 1 in the binary signal, t represents time, and s(n) represents the power timing signal of n loads; The time domain expression of amplitude keying modulation is: Among them, the pulse width of the power timing signal s is Rs, S MASK (t) is the time domain of amplitude keying modulation, g represents the matrix pulse, a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, and t represents the time; The amplitude keying modulation time domain is summed to obtain the amplitude modulated power timing signal 4. A secure aggregation method for power time series data communication according to claim 3, characterized in that: The phase shift keying modulation comprises the following steps: Phase modulating the amplitude modulated power time series signal through the fitted transmittable power time series signal to obtain a phase modulated power time series signal; The time domain expression of phase shift keying modulation is set as: Among them, S MPSK (t) is the time domain of phase shift keying modulation, g represents the matrix pulse, the pulse width of the power timing signal s is Rs, a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, t represents the time, It is represented as the phase corresponding to the nth load; in: The phase-shift keying modulated time domain is summed to obtain the phase modulated power timing signal 5. A secure aggregation method for power time series data communication according to claim 4, characterized in that: The frequency shift keying modulation comprises the following steps: The phase-modulated power timing signal is modulated to a different frequency by frequency shift keying modulation, keeping the amplitude and phase unchanged; The time domain expression of frequency shift keying modulation is set to: Among them, S MFSK (t) is the time domain of frequency shift keying modulation, g represents the matrix pulse, the pulse width of the power timing signal s is Rs, a η is the level value of the nth load, f c It is represented as the carrier frequency of the cth power timing signal, t represents time, Δf c is the carrier frequency offset; The frequency shift keying modulated time domain is summed to obtain the power timing signal after frequency shift keying modulation The power timing signal after frequency shift keying modulation is set to be the modulated power timing signal.

6. A system for implementing the secure aggregation method for power time series data communication according to any one of claims 1 to 5, characterized in that: It includes a data acquisition module, a database, a data processing module, a data communication decomposition module, a data communication fitting module, a data communication modulation module and a data communication aggregation module; The data acquisition module is used to collect historical power time series data to establish a power time series data sample set and to collect power time series data in real time; The database is used to store the power time series data collected by the data acquisition module; The data processing module is used to pre-process the collected power time series data; The data communication decomposition module is used to decompose the preprocessed power time series data; The data communication fitting module is used to fit the decomposed power time series data; The data communication modulation module is used to modulate the fitted power time series data; The data communication aggregation module is used to aggregate the modulated power time series data.

Images