CN109992805A - Method and device for simulating wind speed - Google Patents
Method and device for simulating wind speed Download PDFInfo
- Publication number
- CN109992805A CN109992805A CN201711487726.2A CN201711487726A CN109992805A CN 109992805 A CN109992805 A CN 109992805A CN 201711487726 A CN201711487726 A CN 201711487726A CN 109992805 A CN109992805 A CN 109992805A
- Authority
- CN
- China
- Prior art keywords
- groups
- sequence
- time
- random
- wind speed
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000001914 filtration Methods 0.000 claims abstract description 21
- 238000001228 spectrum Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims abstract description 10
- 230000006870 function Effects 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 230000005055 memory storage Effects 0.000 claims description 2
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims 2
- 238000004088 simulation Methods 0.000 abstract description 14
- 230000003595 spectral effect Effects 0.000 abstract description 6
- 238000013211 curve analysis Methods 0.000 abstract description 4
- 238000010248 power generation Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 20
- 238000004590 computer program Methods 0.000 description 10
- 230000005611 electricity Effects 0.000 description 6
- 230000006854 communication Effects 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 239000004744 fabric Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 239000004065 semiconductor Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000011109 contamination Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000011017 operating method Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
- G06F17/156—Correlation function computation including computation of convolution operations using a domain transform, e.g. Fourier transform, polynomial transform, number theoretic transform
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Abstract
The invention discloses a method and a device for simulating wind speed. The method comprises the following steps: randomly generating two groups of random sequences which accord with normal distribution; respectively carrying out filtering processing on the two groups of random sequences, and correspondingly generating two groups of first time sequences meeting the requirement of a specified power spectrum; respectively converting the two groups of first time sequences into second time sequences based on the specified standard deviation; and generating a wind speed sequence satisfying a Weibull distribution based on the two groups of second time sequences. Therefore, the random wind speed sequence which not only meets the requirement of certain power spectral density, but also meets the requirement of Weibull distribution can be obtained very conveniently, so that the characteristic of wind speed of a certain time scale (such as a middle time scale) can be met more closely, and the simulation requirements of equipment control, power curve analysis and power grid control of wind power generation can be met.
Description
Technical field
The present invention relates to the technical fields of wind-power electricity generation more particularly to a kind of method and apparatus for simulating wind speed and storage to be situated between
Matter.
Background technique
As its people are to the pay attention to day by day of environmental protection, because wind-power electricity generation has cleaning, free of contamination characteristic, wind-force hair
Electricity receives public will be widely welcomed.However, the randomness and fluctuation of wind speed, equipment life, fortune to wind power generating set
Row performance and the stability of power grid have larger impact, thus in research air-blower control strategy, power curve analysis and blower
When accessing the power swing and power quality of power grid, need to establish Wind speed model adaptable therewith.
Currently, the Wind speed model for wind-powered electricity generation emulation specifically includes that basic wind speed, gradual change wind speed, fitful wind and RANDOM WIND are built
Vertical model.
Applicant it has been investigated that, the prior art can only be simulated according to certain basic characteristics of wind speed, when selection
When essential characteristic is 1 (such as basic wind speed feature), although operation is very easy, the precision of simulation is lower；When trial increases
When essential characteristic simulates wind speed, operation is extremely complex, and and is unsatisfactory for the power spectral density of RANDOM WIND and wanting for Weibull distribution
It asks.
How wind speed to be simulated in the way of convenient and fast, it is made to meet the power spectral density of RANDOM WIND and Weibull point
The requirement of cloth becomes technical problem urgently to be resolved.
Summary of the invention
It is in order to solve operation complexity and simulation precision low, it is unable to satisfy the power spectral density and Weibull distribution of RANDOM WIND
It is required that the problem of, the embodiment of the invention provides a kind of method, apparatus for simulating wind speed.
In a first aspect, providing a kind of method for simulating wind speed.Method includes the following steps:
It is random to generate two groups of random sequences for meeting normal distribution；
Two groups of random sequences are filtered respectively, when corresponding two groups of satisfactions of generation specify the first of power spectrum requirement
Between sequence；
It is poor based on specified value, respectively by two groups at the first time it is Sequence Transformed be the second time series；
Based on two group of second time series, the wind series for meeting Weibull distribution are generated.
Second aspect provides a kind of device for simulating wind speed.The device includes:
Random number generation unit, for generating two groups of random sequences for meeting normal distribution at random；
Sequential filtering unit, for being filtered respectively to two groups of random sequences, corresponding two groups of satisfactions of generation are specified
The first time sequence that power spectrum requires；
Scale conversion unit, for based on specified value it is poor, respectively by two groups at the first time Sequence Transformed be second when
Between sequence；
Sequence generating unit, for being based on two group of second time series, generation meets Weibull distribution wind series.
The third aspect provides a kind of device for simulating wind speed.The device includes:
Memory, for storing program；
Processor, for executing the program of the memory storage, it is above-mentioned each that described program executes the processor
Method described in aspect.
Fourth aspect provides a kind of computer readable storage medium.Finger is stored in the computer readable storage medium
It enables, when run on a computer, so that computer executes method described in above-mentioned various aspects.
5th aspect, provides a kind of computer program product comprising instruction.When the product is run on computers,
So that computer executes method described in above-mentioned various aspects.
6th aspect, provides a kind of computer program.When the computer program is run on computers, so that calculating
Machine executes method described in above-mentioned various aspects.
On the one hand, foregoing invention embodiment can be by generating two groups of random sequences for meeting normal distribution, respectively at random
Frequency domain processing is carried out to two groups of random sequences, two groups of first time time sequences for meeting specified power spectrum requirement can be generated；
On the other hand, foregoing invention embodiment can be poor based on specified value, can be Sequence Transformed at the first time by two groups of times respectively
For the time series for meeting specific criteria difference.
It is based on above-mentioned two groups of time serieses as a result, easily can both have been met very much certain power spectral density requirement,
Meet the random wind speed of Weibull distribution requirement again, meets certain time scale (such as middle time scale) so as to more proper
The characteristics of wind speed, and then can satisfy and the equipment of wind-power electricity generation is controlled, the simulation demand of power curve analysis and power grid control.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, will make below to required in the embodiment of the present invention
Attached drawing is briefly described, it should be apparent that, drawings described below is only some embodiments of the present invention, for
For those of ordinary skill in the art, without creative efforts, it can also be obtained according to these attached drawings other
Attached drawing.
Fig. 1 is the flow diagram of the method for the simulation wind speed of one embodiment of the invention；
Fig. 2 is the map schematic diagram of the random sequence S1 of one embodiment of the invention；
Fig. 3 is the map schematic diagram of the random sequence S2 of one embodiment of the invention；
Fig. 4 is the power spectrum schematic diagram of one embodiment of the invention；
Fig. 5 is the map schematic diagram of the random sequence S3 of one embodiment of the invention；
Fig. 6 is the map schematic diagram of the random sequence S4 of one embodiment of the invention；
Fig. 7 is the map schematic diagram of the random sequence S5 of one embodiment of the invention；
Fig. 8 is the map schematic diagram of the random sequence S6 of one embodiment of the invention；
Fig. 9 is the map schematic diagram of the wind series S7 of one embodiment of the invention；
Figure 10 is the structural schematic diagram of the device of the simulation wind speed of one embodiment of the invention；
Figure 11 is the block schematic illustration of the device of the simulation wind speed of one embodiment of the invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
Every other embodiment obtained without creative efforts, shall fall within the protection scope of the present invention.
It should be noted that in the absence of conflict, the features in the embodiments and the embodiments of the present application can phase
Mutually combination.The application is described in detail below with reference to the accompanying drawings and in conjunction with the embodiments.
Fig. 1 is the flow diagram of the method for the simulation wind speed of one embodiment of the invention.
As shown in Figure 1, the method for simulation wind speed may comprise steps of: S110, random two groups of generation meet normal state point
The random sequence (being indicated with S1 and S2) of cloth；S120 is respectively filtered two groups of random sequences, generates and meets specified function
Two groups of first time sequences (being indicated with S3 and S4) that rate spectrum requires；S130, based on specified value it is poor, respectively by two group first when
Between it is Sequence Transformed be the second time series (being indicated with S5 and S6)；S140, is based on two groups of time serieses, and generation meets Weibull point
The wind series (being indicated with S7) of cloth.
In the present embodiment, wind series can be the wind series of middle time scale.Middle time scale week is indicated with Δ t
Phase in the random wind speed period of generation that is, to be simulated, specifically can be 600s (second).Furthermore it is possible to which total week of random wind speed is arranged
Phase, can be counted with data indicates.
In step s 110, can use Generating Random Number obtain two groups of standardized normal distribution random sequence S1 and
S2。
In some embodiments, S1 and S2 can indicate white Gaussian noise.
Fig. 2 is the map schematic diagram of the random sequence S1 of one embodiment of the invention.
As shown in Fig. 2, the lateral numerical value of the S1 map can indicate sequence number, and such as: 1,2 ... N.Wherein, N=
4096.The Vertical Numerical of the map can indicate the value size of random number in the map, and middle time scale period Δ t can be
600s.It can be with one data point of generation in every 10 minutes.Total cycle time can be 4096*600s (i.e. about January).
Fig. 3 is the map schematic diagram of the random sequence S2 of one embodiment of the invention.
As shown in figure 3, the lateral numerical value of the S2 map can indicate sequence number, and such as: 1,2 ... N.Wherein, N=
4096.The Vertical Numerical of the map can indicate the value size of random number.In the map, middle time scale period Δ t can be with
For 600s.It can be with one data point of generation in every 10 minutes.Total periodicity can be 4096*600s (i.e. January).
Referring to figs. 2 and 3, S1 and S2 can be respectively as follows: S1 (a1, a2…aN)~N (0,1), S2 (b1, b2..., bN)~N
(0,1).Wherein ,~N (0,1) indicates to obey (0,1) normal distribution.
This example is filtered by the way of Fourier inversion is at time-domain signal again after to frequency domain filtering processing.For convenience
The implementation of fast Flourier (fft) transformation, N can take 2 integer power, for example, being intended to obtain one month simulation wind speed, can use N
=4096.When one data point of acquisition in 10 minutes, total time can be T=N* Δ t=4096*10*60s.Δ t is the middle time
The scale period is 10*60s (second) here.
In some embodiments, the chronomere of middle time scale includes following one or more: minute, hour and
Month.Such as 10 minutes, 5 minutes and 1 month.
In the step s 120, be filtered respectively to two groups of random sequences includes: to be divided by the way of bandpass filtering
It is other that two groups of random sequences are filtered.
In some embodiments, the mode of bandpass filtering includes one or more of following manner: analog filtering side
Formula (physical filter), digital filtering mode (Fourier inversion method, transfer function method, mobile Return Law etc.).
In some embodiments, it can be taken in signal processing by the filtering method that white Gaussian noise generates gauss heat source model
The method of filtering.
Fig. 4 is the power spectrum schematic diagram of one embodiment of the invention.
As shown in figure 4, lateral numerical value can indicate cycle period, frequency representation can be specifically used.Vertical Numerical can be with table
Show numerical values recited.It may include following 3 wave crests in the figure: atmosphere peak (day synoptic peak), day peak (Diurnal
) and turbulent flow peak (Turbulent peak) peak.fsAnd fhIt can indicate the range of ordinate value on the figure.Vertical line can be with table
Show it is equidistant take a position, take the value at a position to can be the value of response curve (as shown in circled).
In some embodiments, the implementation of Fourier inversion mode may include following S1201-S1204 step:
S1201, the bandwidth f on power spectrums~fhIn range, equidistantly takes N/2 point and P (c is normalized1, c2...,
cN/2), make its variance and is
Wherein,
Wherein,
In expression formula 1, fsIt can be bandwidth on power spectrum, Δ t is time ruler in can indicating in the middle time scale period
Spend the period.
In expression formula 1, fhIt can indicate to count for bandwidth on power spectrum, N, Δ t is that the middle time scale period can
With the time scale period in indicating.
S1202 constructs following first frequency domain sequence: F1 [d0；d1；…；dN-1]=[0；c1(a1+iaN-1)；c2(a1+
iaN-2)；…；cN/2-1aN/2-1+iaN/2+1；cN/2aN/2；cN/2-1(aN/2-1-iaN/2+1)；…c1(a1-iaN-1)]。
In above-mentioned formula, a and c are above-mentioned S1The respective value of sequence and P sequence, i are imaginary symbols.
S1203 constructs the second frequency domain sequence F2 [e0；e；…；eN-1]=[0；c1(b1+ibN-1)；c2(b1+ibN-2)；…；
cN/2-1(bN/2-1+ibN/2+1)；cN/2bN/2；cN/2-1(bN/2-1-ibN/2+1)；…c1(b1-ibN-1)].B is S2 sequence respective value in formula.
S1204, to frequency domain sequence F1 (d0…..dN-1) and frequency domain sequence F2 (e0……eN-1) Fourier's contravariant is carried out respectively
(IFFT) is changed, Fig. 5 and two time domain sequences S3 (f shown in fig. 6 are obtained1…..fn) and S4 (g1……gn).S3 and S4 can be with table
Show gauss heat source model.
With reference to Fig. 5 and Fig. 6, S3 is by F1 by obtaining after Fourier inversion method (IFFT).S4 is by F2 by obtaining after IFFT
It arrives.
Fig. 5 is the map schematic diagram of the random sequence S3 of one embodiment of the invention.
As shown in figure 5, the lateral numerical value of the S3 map can indicate sequence number, and such as: 1,2 ... N.Wherein, N=
4096.Vertical Numerical can indicate that S3 corresponds to numerical value.
Fig. 6 is the map schematic diagram of the random sequence S4 of one embodiment of the invention.
As shown in fig. 6, the lateral numerical value of the S4 map can indicate sequence number, and such as: 1,2 ... N.Wherein, N=
4096.Vertical Numerical can indicate that S4 corresponds to numerical value.
In step s 130, which may include following S1301-S1304 sub-step:
S1301 obtains the Weibull Function of historical wind speed；
S1302 obtains Weibull distribution parameters according to Weibull Function；
S1303 is based on Weibull distribution parameters, and it is poor (such as σ) to obtain specified value；
S1304, it is poor based on specified value, respectively convert two groups of time domain sequences to the time sequence for meeting specific criteria difference
Column.
In the present embodiment, Weibull Function can be as follows:
In above-mentioned expression 3, f (x) can indicate Weibull distribution probability density function；X can be independent variable, i.e. wind
Speed；K and λ is model parameter, and e is natural constant.
In S1303, Weibull distribution parameters are based on, obtain that specified value is poor, may include: obtain the first parameter the
The values of powers of two parameter power；Using the half of values of powers as the square value of specified value difference.
For example, σ2=λk/ 2 (expression formulas 4)
In above-mentioned expression formula 4, σ can be poor for specified value.λ can be the first parameter, and k can be the second parameter.
Fig. 7 is the map schematic diagram of the random sequence S5 of one embodiment of the invention.
As shown in fig. 7, the lateral numerical value of the S5 map can indicate sequence number, and such as: 1,2 ... N.Wherein, N=
4096.Vertical Numerical can indicate that S5 corresponds to numerical value.Wherein, S5 corresponds to numerical value and corresponds to numerical value and the product of σ equal to S3.Wherein, σ
It can be poor for specified value.
Fig. 8 is the map schematic diagram of the random sequence S6 of one embodiment of the invention.As shown in figure 8, the transverse direction of the S6 map
Numerical value can indicate sequence number, such as: 1,2 ... N.Wherein, N=4096.Vertical Numerical can indicate that S6 corresponds to numerical value.Its
In, S6 corresponds to numerical value and corresponds to numerical value and the product of σ equal to S4.Wherein, σ can be poor for specified value.
With reference to Fig. 7 and Fig. 8, can be converted to obtain N (0, σ by S3 and S42) distribution two sequence S5 (h1……hN)、S6
(j1……jN)。
Fig. 9 is the map schematic diagram of the wind series S7 of one embodiment of the invention.
In step S140, which may include following S1401-S1403 sub-step:
S1401, obtain two groups of time serieses square and value；
S1402 handles with value evolution, obtains root value；
S1403 is based on root value, generates wind series.
Specifically, can be to S5 (hn)、S6(jn) sequence carries out that S7 (x is calculated by serial numbern)
In expression formula 5, sequence xnThe random air speed value of simulation output when can be for n (n=1,2 ..., N) digital point；hn
The correspondence numerical value of S5 when can be n digital point；jnThe correspondence numerical value of S6 when can be n digital point.
On the one hand, foregoing invention embodiment can be by generating two groups of random sequences for meeting normal distribution, respectively at random
Frequency domain processing is carried out to two groups of random sequences, two groups of first time time sequences for meeting specified power spectrum requirement can be generated；
On the other hand, foregoing invention embodiment can be poor based on specified value, can be Sequence Transformed at the first time by two groups of times respectively
For the time series for meeting specific criteria difference.
It is based on above-mentioned two groups of time serieses as a result, easily can both have been met very much certain power spectral density requirement,
Meet the random wind speed of Weibull distribution requirement again, meets certain time scale (such as middle time scale) so as to more proper
The characteristics of wind speed, and then can satisfy and the equipment of wind-power electricity generation is controlled, the simulation demand of power curve analysis and power grid control.
It should be noted that in the absence of conflict, those skilled in the art can according to actual needs will be above-mentioned
The sequence of operating procedure is adjusted flexibly, or above-mentioned steps are carried out the operation such as flexible combination.For simplicity, repeating no more
Various implementations.In addition, the content of each embodiment can mutual reference.
Figure 10 is the structural schematic diagram of the device of the simulation wind speed of one embodiment of the invention；
As shown in Figure 10, the device 100 for simulating wind speed may include: random number generation unit 101, sequential filtering unit
102, scale conversion unit 103 and sequence generating unit 104.Wherein, random number generation unit 101 can be used for generating two at random
Group meets the random sequence of normal distribution；Sequential filtering unit 102 can be used for being filtered place to two groups of random sequences respectively
Reason, the corresponding first time sequence for generating two groups of satisfactions and power spectrum being specified to require；Scale conversion unit 103 can be used for based on finger
Standard deviation is determined, respectively by two groups of first time Sequence Transformed second time serieses for being；Sequence generating unit 104 can be used for base
In two group of second time series, generation meets Weibull distribution wind series.
In some embodiments, the method for filtering processing includes: analog filtering method and/or digital filtering.
In some embodiments, which can also include: standard deviation acquiring unit.Standard deviation acquiring unit can be used for
Obtain the Weibull Function of historical wind speed；According to Weibull Function, Weibull distribution parameters are obtained；Based on Weibull
It is poor to obtain specified value for distribution parameter.
In some embodiments, historical wind speed can be the data provided by meteorological observatory, or the number voluntarily measured
According to.For example, historical wind speed can be the wind speed of a period of time (such as 1 year, January, the long time scale of some months, 100 days).One
In a little embodiments, sequence generating unit 104 can be also used for obtaining two groups of time serieses square and value；To with value evolution at
Reason, obtains evolution value；Based on evolution value, wind series are generated.
In some embodiments, random sequence can be white Gaussian noise.
In some embodiments, first time sequence is gauss heat source model.
In some embodiments, wind series can be with are as follows: the wind series of middle time scale.
In some embodiments, the chronomere of middle time scale includes following one or more: minute, hour and
Month.
It should be noted that the device of the various embodiments described above can be used as the method for each embodiment of the various embodiments described above
In executing subject, the corresponding process in each method may be implemented, realize identical technical effect, for sake of simplicity, in this respect
Content repeats no more.
In the above-described embodiments, can come wholly or partly by software, hardware, firmware or any combination thereof real
It is existing.For example, encryption/decryption element is integrated in one unit, two individual units can also be divided into.In another example will request
Receiving unit and request transmitting unit are substituted with a coffret.When implemented in software, can entirely or partly with
The form of computer program product is realized.The computer program product includes one or more computer instructions, when it is being counted
When being run on calculation machine, so that computer executes method described in above-mentioned each embodiment.Load and execute on computers institute
When stating computer program instructions, entirely or partly generate according to process or function described in the embodiment of the present invention.The calculating
Machine can be general purpose computer, special purpose computer, computer network or other programmable devices.The computer instruction can
To store in a computer-readable storage medium, or computer-readable deposit from a computer readable storage medium to another
Storage media transmission, for example, the computer instruction can pass through from a web-site, computer, server or data center
Wired (such as coaxial cable, optical fiber, Digital Subscriber Line (DSL)) or wireless (such as infrared, wireless, microwave etc.) mode are to another
A web-site, computer, server or data center are transmitted.The computer readable storage medium can be computer
Any usable medium that can be accessed either includes the data such as one or more usable mediums integrated server, data center
Store equipment.The usable medium can be magnetic medium, (for example, floppy disk, hard disk, tape), optical medium (for example, DVD) or
Person's semiconductor medium (such as solid state hard disk Solid State Disk, SSD) etc..
Figure 11 is the block schematic illustration of the device of the simulation wind speed of one embodiment of the invention.
As shown in figure 11, which may include central processing unit (CPU) 1101, can be according to being stored in read-only deposit
Program in reservoir (ROM) 1102 is loaded into the program in random access storage device (RAM) 1103 from storage section 1108
And execute the various operations that each embodiment in Fig. 1 is done.In RAM1103, also it is stored with each needed for system architecture operation
Kind program and data.CPU 1101, ROM 1102 and RAM 1103 are connected with each other by bus 1104.Input/output (I/O)
Interface 1105 is also connected to bus 1104.
I/O interface 1105 is connected to lower component: the importation 1106 including keyboard, mouse etc.；Including such as cathode
The output par, c 1107 of ray tube (CRT), liquid crystal display (LCD) etc. and loudspeaker etc.；Storage section including hard disk etc.
1108；And the communications portion 1109 of the network interface card including LAN card, modem etc..Communications portion 1109 passes through
Communication process is executed by the network of such as internet.Driver 1110 is also connected to I/O interface 1105 as needed.It is detachable to be situated between
Matter 1111, such as disk, CD, magneto-optic disk, semiconductor memory etc. are mounted on as needed on driver 1110, so as to
In being mounted into storage section 1108 as needed from the computer program read thereon.
Particularly, according to an embodiment of the invention, may be implemented as computer above with reference to the process of flow chart description
Software program.For example, the embodiment of the present invention includes a kind of computer program product comprising be tangibly embodied in machine readable
Computer program on medium, the computer program include the program code for method shown in execution flow chart.At this
In the embodiment of sample, which can be downloaded and installed from network by communications portion 1109, and/or from removable
Medium 1111 is unloaded to be mounted.
The apparatus embodiments described above are merely exemplary, wherein described, unit can as illustrated by the separation member
It is physically separated with being or may not be, component shown as a unit may or may not be physics list
Member, it can it is in one place, or may be distributed over multiple network units.It can be selected according to the actual needs
In some or all of the modules achieve the purpose of the solution of this embodiment.Those of ordinary skill in the art are not paying creativeness
Labour in the case where, it can understand and implement.
Through the above description of the embodiments, those skilled in the art can be understood that each embodiment can
It realizes by means of software and necessary general hardware platform, naturally it is also possible to pass through hardware.Based on this understanding, on
Stating technical solution, substantially the part that contributes to existing technology can be embodied in the form of software products in other words, should
Computer software product may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, CD, including several fingers
It enables and using so that a computer equipment (can be personal computer, server or the network equipment etc.) executes each implementation
Method described in certain parts of example or embodiment.
Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations；Although
Present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that: it still may be used
To modify the technical solutions described in the foregoing embodiments or equivalent replacement of some of the technical features；
And these are modified or replaceed, technical solution of various embodiments of the present invention that it does not separate the essence of the corresponding technical solution spirit and
Range.
Claims (17)
1. a kind of method for simulating wind speed, which comprises the following steps:
It is random to generate two groups of random sequences for meeting normal distribution；
The random sequence described in two groups is filtered respectively, when corresponding two groups of satisfactions of generation specify the first of power spectrum requirement
Between sequence；
It is poor based on specified value, it is respectively the second time series by first time is Sequence Transformed described in two groups；
Based on the second time series described in two groups, the wind series for meeting Weibull distribution are generated.
2. the method according to claim 1, wherein the method for the filtering processing include: analog filtering method and/
Or digital filtering.
3. the method according to claim 1, wherein based on specified value it is poor, respectively by first described in two groups when
Between it is Sequence Transformed be the second time series before, further includes:
Obtain the Weibull Function of historical wind speed；
According to the Weibull Function, Weibull distribution parameters are obtained；
Based on the Weibull distribution parameters, it is poor to obtain the specified value.
4. the method according to claim 1, wherein it is described be based on two groups described in the second time series, generate full
The wind series of sufficient Weibull distribution, comprising:
Time series described in obtaining two groups square and value；
To described and value evolution processing, root value is obtained；
Based on the evolution value, the wind series are generated.
5. the method according to claim 1, wherein wherein；
The random sequence is white Gaussian noise.
6. the method according to claim 1, wherein wherein；
The first time sequence is gauss heat source model.
7. method according to claim 1 to 6, which is characterized in that the wind series are as follows: middle time scale
Wind series.
8. the method according to the description of claim 7 is characterized in that wherein:
The chronomere of the middle time scale includes following one or more: minute, hour, day and the moon.
9. a kind of device for simulating wind speed characterized by comprising
Random number generation unit, for generating two groups of random sequences for meeting normal distribution at random；
Sequential filtering unit is filtered for the random sequence described in two groups respectively, and corresponding two groups of satisfactions of generation are specified
The first time sequence that power spectrum requires；
Scale conversion unit was respectively the second time by first time is Sequence Transformed described in two groups for poor based on specified value
Sequence；
Sequence generating unit, for generating the wind series for meeting Weibull distribution based on the second time series described in two groups.
10. device according to claim 9, which is characterized in that the method for the filtering processing includes: analog filtering method
And/or digital filtering.
11. device according to claim 9, which is characterized in that further include:
Standard deviation acquiring unit, for obtaining the Weibull Function of historical wind speed；According to the Weibull Function, obtain
Take Weibull distribution parameters；Based on the Weibull distribution parameters, it is poor to obtain the specified value.
12. device according to claim 9, which is characterized in that wherein:
The sequence generating unit is also used to: obtain two groups described in time series square and value；At described and value evolution
Reason, obtains root value；Based on the root value, the wind series are generated.
13. device according to claim 9, which is characterized in that wherein；
The random sequence is white Gaussian noise.
14. device according to claim 9, which is characterized in that wherein；
The first time sequence is gauss heat source model.
15. the device according to any one of claim 9-14, which is characterized in that the wind series are as follows: middle time ruler
The wind series of degree.
16. according to the method for claim 15, which is characterized in that wherein:
The chronomere of the middle time scale includes following one or more: minute, hour, day and the moon.
17. a kind of device for simulating wind speed characterized by comprising
Memory, for storing program；
Processor, for executing the program of the memory storage, described program makes the processor execute such as claim
Method described in any one of 1-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711487726.2A CN109992805B (en) | 2017-12-29 | 2017-12-29 | Method and device for simulating wind speed |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711487726.2A CN109992805B (en) | 2017-12-29 | 2017-12-29 | Method and device for simulating wind speed |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109992805A true CN109992805A (en) | 2019-07-09 |
| CN109992805B CN109992805B (en) | 2022-10-18 |
Family
ID=67111442
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711487726.2A Active CN109992805B (en) | 2017-12-29 | 2017-12-29 | Method and device for simulating wind speed |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109992805B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10154200C1 (en) * | 2001-11-07 | 2003-03-06 | Infineon Technologies Ag | Random number sequence generation method for simulation of noise using computer program product and computer program |
| CN105119320A (en) * | 2015-09-15 | 2015-12-02 | 东北大学 | Distributed wind power plant fan optimized arrangement system and method |
| CN106500648A (en) * | 2016-12-08 | 2017-03-15 | 北京国网富达科技发展有限责任公司 | Power transmission circuit caused by windage monitoring method and device based on dynamic wind |
-
2017
- 2017-12-29 CN CN201711487726.2A patent/CN109992805B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE10154200C1 (en) * | 2001-11-07 | 2003-03-06 | Infineon Technologies Ag | Random number sequence generation method for simulation of noise using computer program product and computer program |
| CN105119320A (en) * | 2015-09-15 | 2015-12-02 | 东北大学 | Distributed wind power plant fan optimized arrangement system and method |
| CN106500648A (en) * | 2016-12-08 | 2017-03-15 | 北京国网富达科技发展有限责任公司 | Power transmission circuit caused by windage monitoring method and device based on dynamic wind |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109992805B (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cui et al. | Characterizing diffused stellar light in simulated galaxy clusters | |
| Ismael et al. | Probabilistic hosting capacity enhancement in non-sinusoidal power distribution systems using a hybrid PSOGSA optimization algorithm | |
| Zhao et al. | Accurate universal models for the mass accretion histories and concentrations of dark matter halos | |
| Larik et al. | An improved algorithm for optimal load shedding in power systems | |
| Edelman et al. | Cover Your Basis: Comprehensive Data-driven Characterization of the Binary Black Hole Population | |
| Hong et al. | Hour-ahead wind speed and power forecasting using empirical mode decomposition | |
| CN111028893A (en) | Crystal growth prediction method and device | |
| Dovgun et al. | Distributed power quality conditioning system for three-phase four-wire low voltage networks | |
| Basu et al. | The MLP distribution: a modified lognormal power-law model for the stellar initial mass function | |
| Andreadou et al. | Smart meter traffic in a real LV distribution network | |
| Kitiashvili | Data assimilation approach for forecast of solar activity cycles | |
| Chen et al. | An energy storage assessment: Using frequency modulation approach to capture optimal coordination | |
| Bashir et al. | Development of frequency weighted model reduction algorithm with error bound: application to doubly fed induction generator based wind turbines for power system | |
| Das et al. | Performance of hybrid filter in a microgrid integrated power system network using wavelet techniques | |
| Zhang et al. | A novel routing algorithm for power line communication over a low-voltage distribution network in a smart grid | |
| Frontiere et al. | Simulating hydrodynamics in cosmology with crk-hacc | |
| Liaquat et al. | Modeling and formulation of optimization problems for optimal scheduling of multi-generation and hybrid energy systems: Review and recommendations | |
| Liu et al. | A comprehensive survey of accurate and efficient aggregation modeling for high penetration of large-scale wind farms in smart grid | |
| Cieślik | Mathematical modeling of the dynamics of linear electrical systems with parallel calculations | |
| Capligins et al. | FPGA-Based antipodal chaotic shift keying communication system | |
| CN108229676A (en) | The training method and device of neural network model, electronic equipment and computer storage media | |
| Priyadarshini et al. | Significance of harmonic filters by computation of short-time Fourier transform-based time–frequency representation of supply voltage | |
| Conti et al. | A multi-port hardware energy meter system for data centers and server farms monitoring | |
| CN109992805A (en) | Method and device for simulating wind speed | |
| Kołek et al. | Analysis of the practical implementation of flicker measurement coprocessor for AMI meters |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |