METHOD AND APPARATUS FOR ANALYZING PHYSICAL QUANTITIES, AND APPARATUS FOR REMOVING LINE SPECTRUM NOISE

Abstract: A physical formula the principal portion of which is determined on the basis of physical quantities is calculated by a plurality of physical formula calculation units (11), (12), ...(1m) on the basis of known information, a sigma unit (2) effects cumulative addition, an error calculator (3) calculates the difference from an actually measured value, and correction units (11a), (12a), ...(1ma) of the physical formula calculation units (11), (12), ...(1m) correct variables contained in the physical formula so that the difference described above becomes small. A series of processings described above are repeated a necessary number of times, and the corrected variables are collected by an information collection unit (4) and are outputted as an analytical result of the physical quantities. As a result, the construction can be simplified, and analysis of the physical quantities of a physical source can be accomplished with a smaller number of times of learning on the basis of the measured values.

Description: Technical Field This invention relates to method and apparatus for analyzing physical quantities and apparatus for reducing line spectrum noise. More particularly, the present invention relates to method and apparatus for calculating physical quantities of physical sources based upon measured values of physical quantities which are obtained at plural points apart from the physical sources, and observation condition at measuring timing, and relates to apparatus for reducing line spectrum noises from measured signal which includes line spectrum noises. In the following a description of background technology is given which facilitates the understanding of the invention. From past years, an apparatus is proposed which analyzes biological magnetic field sources (one species of physical sources) by disposing plural Superconducting Quantum Interference Devices (hereinafter referred to as SQUID) magnetometers in a closed condition to the living organism, the SQUID magnetometer employing SQUIDs. The apparatus performs following processings using a supercomputer. That is, a) Scattering m-numbered current elements using random numbers within a search space which is searched by the plural SQUID magnetometers. Input parameters of a current source i are position information P(x,y,z) and a current vector I(X,Y,Z), thereby 6m-numbered parameters xi,yi,zi,Xi,Yi,Zi (i=1,2,...m) are determined using random numbers. b) Calculating a total estimated error E by an estimated error calculating process which is described later. c) Repeat following procesings d) to g). d) Selecting a current element arbitrarily, and evacuating parameters of corresponding current element k and the total estimated error. That is, processings of position information Ps (xs, ys, zs) = Pk (xk, yk, zk), a current vector Is (Xs, Ys, Zs) = Ik (Xk, Yk, zk), a total estimated error Es=E are carried out. Wherein a suffix k indicates data of a current element k, and a suffix s indicates evacuated data. e) Varying parameters of the current element k by extremely small quantities using random numbers. That is, when extremely small quantities are supposed as Δx, Δy, Δz, ΔX, ΔY, ΔZ, processings of Pk (xk, y k, z k) = Pk (xk + Δx, y k + Δy, zk + Δz) Ik (Xk, Yk, Zk) = (Xk + ΔX, Yk + ΔY, Zk + ΔZ) are carried out. f) Calculating a total estimated error E by the estimated error calculating process which is described later. g) Comparing the evacuated total estimated error Es and the total estimated error E, and restoring the information which are evacuated at the processing d) when the total estimated error Es is smaller than the total estimated error E. That is, processings of position information pk (xk, yk, zk) = pk (xs, ys, zs), a current vector Ik (Xk, Yk, Zk) = Is (Xs, Ys, Zs), a total estimated error E=Es are carried out. The estimated error calculating process is as follows. I. Magnetic field at each measuring point j (j=1,2,...N) is calculated based upon parameters of each current element. That is, 1) Carrying out following processings 2) and 3) for all measuring points j. 2) Carrying out following processing 3) for all current elements i. 3) Calculating magnetic field Beij (BXeij,BYeij,BZeij) using Biot-Savart law, the magnetic field Beij being generated at a measuring point j by a current element i. 4) Calculating magnetic field Bej using the following equation, the magnetic field Bej being generated at the measuring point j by m-numbered current elements. II. Calculating an estimated error Ej based upon a measured value Bj(BXj,BYj,BZj) at each measuring point j and an estimated value Bej based upon all current elements, and calculating a total estimated error E. That is, 5) Carrying out following processing 6) for all measuring points j. 6) Calculating the estimated error Ej at each measuring point j using the following equation. Ej = (BXj - BXej) 2 + (BYj - BYej) 2 + (BZj - BZej) 2 7) Calculating the total estimated error E using the following equation. When magnetic field sources are analyzed using the apparatus mentioned above, correct analysis results seem to be obtained finally because parameters of current elements k are varied by extremely small quantities so as to decrease the total estimated error E. But, when initial condition of current elements k are determined as is illustrated in Fig. 22(A), the current elements k are varied to a condition illustrated in Fig. 22(B) when 2,400 times of processings are carried out, and the current elements k are varied to a condition illustrated in Fig. 22(C) when 3,600 times of processings are carried out. Consequently, final solutions cannot be obtained. When Figs. 22(B)(C) are compared to one another, condition of the current element k is varied very slightly, therefore a disadvantage arises in that final solution cannot be obtained even when a number of times of processings is increased. Even if 3,600 times of processings are carried out using a supercomputer, it takes about 20 minutes, therefore the apparatus cannot be utilized at all. The inventor found causes for arising the disadvantages. The causes are as follows: When the processings are carried out for once, it is not assured at all to decrease the total estimated error E, and Only processing parts of I.1),2),3) and II.5),6) of the estimated error calculating processings can be processed in parallel, while other processing parts cannot be processed in parallel, therefore speed up of calculation as a whole cannot be achieved even when parallel processor system is employed. In recent years, research of an artificial neural network is developed, and it is thought to apply the artificial neural network to the analysis of magnetic field sources mentioned above. The artificial neural network is classified into a hierarchical perceptron (refer to Fig. 23) and a Hopfield model (refer to Fig. 26). A hierarchical perceptron is constituted by an input layer for receiving an input pattern, one or more intermediate layers and an output layer for outputting an output pattern, and a plurality of neuron devices constituting each layer are interconnected to one another, as is illustrated in Fig. 23. A back propagation rule is employed as a learning rule of the hierarchical perceptron. But, when a problem to be dealt with the hierarchical perceptron becoms complicated, a large-numbered-layers construction is required, or increase of a number of neuron devices which constitute one intermediate layer is required. Consequently, a number of neuron devices as a whole becomes great, a number of weighting factors which should be determined by learning becomes extremely great, and an operation load for converging solution becomes extremely great. Specifically, when learning is carried out for about 50 times for one pattern, an error suddenly decreases so as to be settled to a first convergence value, as is illustrated in Fig. 24. When learning is continued, an error may suddenly varies for several times. Fig. 25 illustrates varying of an error for each pattern. In Fig. 25, compromising of allocation of errors continues (refer to solid lines in Fig. 25) even at a portion in which an error as a whole (refer to dashed line in Fig. 25) scarcely exists, therefore trial for decreasing an error as a whole continues. Consequently, when an objected accuracy cannot be obtained by a first convergence value, great amount of processings to the next convergence is required, and it is not assured to obtain a sufficient accuracy by the next convergence value. As is apparent from the foregoing, it is almost impossible to employ a hierarchical perceptron for analyzing magnetic field sources in practical use because of limitations concerning to learning process of physical phenomenon. A Hopfield model has an arrangement in which each neuron model is interconnected to all other neuron models, and has no classification such as input layer, intermediate layer and output layer of a hierarchical perceptron, as is illustrated in Fig. 26. And all neuron models can perform functions as input layer, intermediate layer, and output layer. When a neuron model is supposed to be a threshold device model, condition varying of a neuron model i can be modeled by one of following formulae. Ui represents an activity or output value of the neuron model i, hi represents a threshold value of the neuron model i, Wij represents a weighting factor, and i is not equal to j. The function (αi is a positive constant, and i is not equal to j in the first paragraph of the above-mentioned formula) which is defined as an evaluation function for determining weighting factors and a learning rule of a neuron model, by Hopfield, decreases the formula (2) based upon varying of each neuron model when neuron models perform varying of interior condition based upon formulae (1) in asynchronous manner. When formula (2) becomes a minimum value or a local minimum value, activities or output values of neuron models converge. When analysis is to be carried out using a Hopfield model, following works 1) to 4) for preparation are necessary. That is, 1) Converting a specific object function which is given for analysis, to a form of formula (2). 2) Converting independent variables specifically constituting the object function so as to be activities or output values Ui of neuron models, and determining a conversion rule of the activities or output values Ui, the conversion rule being similar to the formula (1). 3) Determining an input pattern to each Ui, that is initial value of each Ui, so as to be able to be converged to a pattern finally. 4) Repeating same information processing using the learning rule obtained at 1) and 2), till an output pattern converges or a function En becomes minimum. Among these works, works 1) and 2) are very difficult and it is highly possible not to be converted suitably, therefore extremely great amount of work is needed. As to the work 3), the output pattern may converge or not depending upon a given input pattern. Therefore, it is difficult to determine an input pattern which gives high convergence, and different solution may be obtained depending upon a given input pattern. As to the work 4), it takes a fairly long time from varying of inner state of one of neuron models to a new stable state settled by propagating an unbalanced state caused by the variation to the entirety of the Hopfield model. Therefore, it may take an extremely long time for obtaining convergence values by performing information processing in asynchronous manner. As is apparent from the foregoing, it is almost impossible to employ a Hopfield model for analyzing magnetic field sources in practical use. Only application for analysis of magnetic field sources is described in the foregoing, similar disadvantages arise when the system is applied to analysis of physical sources, such as pressure sources and temperature sources, a regulating rule of a system including the physical sources being able to be expressed by numerical formulae, and linear addition being realized in the system. In Superconductor Science and Technology, September 1991, pages 377 - 385 "DC SQUID systems overview" a SQUID system is described consisting of three main components: the pick-up circuit, the squid and the feedback electronics for the linearization of the transfer function. An application to biomagnetic multi-channel systems is described using DC SQUIDs with integrated pick-up loops and alternative read-out schemes. It is the object of the present invention to provide a method and an apparatus for analysing physical quantities with ease and with high speed based upon measured values obtained at plural points. To perform the object above mentioned a method for analysing physical quantities according to claim 1 is provided. As to the method for analyzing physical quantities according to claim 1, plural physical formulae are operated based upon known information, the physical formulae being determined corresponding to species of physical quantities (that is, the physical formulae have different constants from one another), and an estimated value corresponding to a measured value is obtained by cumulatively adding operation results of the physical formulae, when physical quantities of physical sources which are measurable at arbitrary points apart from individual physical sources and satisfy linear addition, are analyzed based upon physical quantities measured at plural points. Then, a difference between measured physical quantity and a value obtained by cumulatively adding the operation results of the physical formulae is obtained, variables included in each physical formula are corrected based upon the calculated difference so as to vary variables to decrease the difference. Thereafter, variables are further varied to decrease the difference to a sufficiently small value by repeating the above-mentioned series of processings based upon the varied variables. When the above-mentioned series of processings are repeated till the difference becomes sufficiently small value, the estimated value can be approximated to a measured value with extremely high accuracy. Thereby, analysis results of physical quantities can be obtained by outputting the corrected variables included in the physical formulae at the timing as analysis results of physical quantities. In this method, it is sufficient that constants included in known physical formulae which constants should be determined corresponding to physical quantities, are varied corresponding to the difference between the cumulatively added result and the measured result, therefore processings as a whole can be simplified and required time period for analysis can extremely be shortened. Analysis of physical quantities can be performed using a conventional artificial neural network and using a measured result as a teacher signal. But, operation formulae themselves should be determined by learning, and required time period for studying becomes extremely long. It is of course that the conventional artificial neural network cannot be applied when realtime learning is required. The conventional artificial neural network cannot be applied to analysis of physical quantities when physical rule to be treated is complicated. It is not assured that analysis result of physical quantities with predetermined accuracy is obtained, because operation formulae themselves are determined by learning. It can be thought that analysis of physical quantities are performed after operation formulae which are for analyze physical quantities from measued result based upon known physical formulae. A disadvantage arises in that the operation formulae become extremely complicated following increase of a number of physical sources, or that proper operation formulae may not obtained. The method employing a conventional artificial neural network can be applied to only a specific system, and is required to redetermine operation formulae when a specific system is changed to another system. A disadvantage arises in that general use is scarcely performed. On the contrary, the invention accodingly to claim 1 applies known physical formulae as they are and performs operation. Therefore, determination of physical formulae is easy and operation results are obtained with accuracy. And, required time period for analysis can extremely be shortened, because only constants which are inclused in the physical formulae, and which constants should be determined corresponding to physical quantities, are corrected. The object is also solved by an apparatus having the features of claim 2. As to the appartatus for analyzing physical quantities according to claim 2, operations based upon the physical formulae are carried out by the physical formula operating means of which number is greater than the number of cause of physical quantities which are objected for analysis, operation results output from each physical formula operating means are cumulatively added by the cumulative addition means, an error is calculated by the error calculating means by receiving the cumulatively added operation result output from the cumulative addition means and a measured value of a physical quantity, and correction at the physical formula operating means is carried out by feedbacking the calculated error to the physical formula operating means, when physical quantities of physical sources which are able to be calculated based upon predetermined operation formulae at arbitrary points apart from individual physical sources and satisfy linear addition, are calculated based upon physical quantities measured at plural points. Thereafter, results, which are corrected based upon the calculated error of causes of physical quantities in each physical formula operating means, are collected and output as analysis results of physical quantities of physical sources by the corrected result collecting means. More particularly, when physical quantities of physical sources which are able to be calculated based upon predetermined operation formulae at arbitrary points apart from individual physical sources, and when physical quantities generated by each physical source at observation points satisfy linear addition, are calculated based upon physical quantities measured at plural points, a plurality of known information exist such as mesuring time, measuring location, and the like. Therefore, the invention supplies each known information to the physical formula operating means so as to obtain corresponding operation results under the condition that variables included operation formulae of the physical formula operating means are determined arbitrary values. The invention then obtains a cumulatively added value, that is a value corresponding to a measured value of a physical quantity by supplying the operation results output from each physical formula operating means to the cumulative addition means. In this case, the obtained cumulatively added value is highly possible not to coincident with the measured value of physical quantity, because the variables are determined arbitrarily. But, the invention supplies the cumulatively added value and the actual measured value of physical quantity to the error calculating means so as to calculate an error between both values, and feedbacks the calculated error to each physical formula operating means so as to vary the variables of each physical formula operating means to decrease the error by correcting determined values of the variables corresponding to the error. Consequently, the error between both values can extremely be small by repeating the above-mentioned series of processings for a required number of times. And, quantities of physical sources can be obtained by collecting and outputting the variables by the corrected result collecting means, the variables being determined at each physical formula operating means at the timing. In this case, it may be thought that the finally obtained analysis result may be a local minimum. But this invention obtains an accurate analysis result other than a local minimum finally, because all of variables of the physical formula operating means are varied based upon the calculated error so as to vary the variables to decrease fluctuation at the entirety of a system. Furthermore, quantity of operations can extremely be decreased in comparison with a conventional method, because required operations are only operations in the physical formula operating means, cumulative addition in the cumulative addition means, error calculation in the error calculating means, and correcting operations of variables by correcting means based upon the feedbacked error. Further developments of the apparatus are given in the Subclaims. As to the appartatus for analyzing physical quantities according to claim 3, plural experimemtal formulae, constant values of each experimental formula being different from one another, are operated by supplying known information (for example, time) to plural physical formula operating means, the known information being obtained after beginning immunity measurement by introducing the exiting light in the light guide. An estimated intensity of fluorescence corresponding to measured intensity of fluorescence is obtained by cumulatively adding the operation results of experimemtal formulae. Then, a difference between the obtained estimated intensity of fluorescence and the measured intensity of fluorescence by the error calculating means. Constant values in each physical formula operating means are corrected by feedbacking the calculated error to each physical formula operating means. And, above-mentioned series of processings are repeated till the error calculated by the error calculating means becomes sufficiently small. When the error becomes sufficiently small, the corrected constants correspond at least to immunity fluorescence and non-specific adsorption fluorescence are output as immunity measurement results by the corrected results collecting means. In a conventional method, immunity measurement accuracy is lowered by being influenced of fluorescence caused by a light guide and the like, and it takes a long time till immunity measurement results are obtained. On the contrary, the invention can obtain immunity measurement results with high accuracy and with short time by using only data which are obtained at initial stage of immunity measurement. As to the appartatus for analyzing physical quantities according to claim 4, plural formulae, constant values of each formula being different from one another, are operated by supplying known information (for example, time) to plural physical formula operating means, the known information being obtained after depositing test solution onto the enzyme-immobilized membrane. An estimated electric signal corresponding to measured electric signal, which corresponds to quantity of material generated or consumed by the reaction of the material included in the test solution, is obtained by cumulatively adding the operation results of formulae. Then, a difference between the obtained estimated electric signal and the measured electric signal by the error calculating means. Constant values in each physical formula operating means are corrected by feedbacking the calculated error to each physical formula operating means. And, above-mentioned series of processings are repeated till the error calculated by the error calculating means becomes sufficiently small. When the error becomes sufficiently small, the corrected constants correspond at least to concentration of material, which reacts under the existence of enzyme, are output as measured concentration results of the material, which reacts under the existence of enzyme, by the corrected results collecting means. In a conventional method, concentration measurement accuracy is lowered by being influenced of thickness of membranes, mounting condition of membranes onto a base electrode and the like, and it takes a long time till measured concentration results are obtained when measured concentration signal is waited till the electric signal becomes stable. On the contrary, the invention can obtain measured concentration results with high accuracy and with short time by using only data which are obtained at initial stage of concentration measurement. As to the appartatus for analyzing physical quantities according to claim 5, a wave is radiated to physical sources by the wave radiating means, and reflected waves from the physical sources are received by the wave receiving means so as to obtain measurement values of the physical quantities. Operations of the physical formulae are carried out by the physical formula operating means of which number is greater than the number of causes of physical qiantities which are objected for analysis, operation results output from each physical formula operating means are cumulatively added by the cumulative addition means, an error is calculated by the error calculating means by receiving the cumulatively added result output from the cumulative addition means and a measured value of physical quantity, and correction at the physical formula operating means is carried out by the correcting means, when physical quantities of physical sources which are able to be calculated based upon predetermined operation formulae at arbitrary points apart from individual physical sources and satisfy linear addition, are calculated based upon physical quantities measured at plural points, the the physical quantity to be analyzed being physical characteristics of the physical source. Thereafter, results, which are corrected based upon the calculated error of causes of physical quantities in each physical formula calculating means, are collected and output as analysis results of physical quantities of physical sources by the corrected result collecting means. Therefore, the invention can obtain measurement results of physical quantities with high accuracy and with short time as similar as the case of claim 2. The reason is that equivalent condition, in which physical sources themselves radiate some physical quantities externally, can be generated by radiating corresponding wave to the physical sources so as to be reflected by the physical sources, when the physical quantities objected for analysis are quantities which are not actively radiated externally from the physical sources, for example, the physical quantities objected for analysis are infrared adsorption rate, ultrasonic reflection rate, and others. As to the apparatus for reducing line spectrum noise according to claim 6, operations of the physical formulae corresponding to each line spectrum noises are performed by physical formula operating means, a number of physical formula operating means corresponding to a number of species of line spectrum noises, and operation results output from each physical formula oparating means are cumulatively added by the cumulative addition means. The cumulatively added result output from the cumulative addition means and the measured spectrum signal are supplied to the error calculating means so as to calculate the error of the cumulatively added result with respect to the measured result, and the value of each physical formula operating means is corrected based upon the calculated error by the correcting means, the value regulating a cause of the line spectrum. After correction of the value regulating a cause of the line spectrum is funished, a signal reduced line spectrum noises from the spectrum signal is output from the error calculating means by interrupting only correcting processing by the correcting means. Therefore, it is no need to collect values corrected by each correcting means. The apparatus for reducing line spectrum noise can be more simplified its arrangement as a whole, and can obtain an observation signal with high quality which is reduced only line spectrum noises. Brief Description of the Drawings Figure 1 is a block diagram illustrating a physical quantity analyzing apparatus according to an embodiment of the present invention; Figure 2 is a diagram illustrating an example of an estimated error evaluating function; Figure 3 is a schematic view illustrating an example of a immunity measuring apparatus using a fluorescence; Figure 4 is a schematic diagram illustrating varying of a measured signal following passage of time, the measured signal being obtained by the immunity measuring appatarus illustrated in Fig. 3; Figure 5 is a flowchart explaining an immunity, measuring processing using the physical quantity analyzing apparatus illustrated in Fig. 1; Figure 6 is a schematic view illustrating a glucose concentration measuring apparatus; Figure 7 is a schematic diagram illustrating varying of a measured signal following passage of time, the measured signal being obtained by the glucose concentration measuring appatarus illustrated in Fig. 6; Figure 8 is a schematic block diagram illustrating an arrangement of a magnetic source analyzing apparatus using the apparatus illustrated in Fig. 1; Figure 9 is a schematic block diagram illustrating an arrangement of an eco searching apparatus using the apparatus illustrated in Fig. 1; Figure 10 is a diagram schematically illustrating relationships of wave radiating point, reflecting points and wave receiving points; Figure 11 is a block diagram illustrating an apparatus for reducing line spectrum noise; Figure 12 is a diagram illustrating a converging condition of one line spectrum noise corresponding to variation of initial values; Figure 13 is a diagram illustrating a converging condition of plural line spectrum noises; Figure 14 is a diagram illustrating a converging condition of plural line spectrum noises; Figure 15 is a block diagram illustrating an apparatus for reducing line spectrum noise which includes a direct current offset reducing function; Figure 16 is a block diagram illustrating a magnetocardiogram measuring apparatus; Figure 17 is a diagram illustrating an example of a rectangular window for extracting a Q-R-S group of an ECG; Figure 18(A) is a diagram illustrating a waveform of an MCG; Figure 18(B) is a diagram illustrating a waveform of an ECG; Figure 18(C) is a diagram illustrating a waveform of an MCG in which spectrum noises are mixed; Figures 19(A) to (F) are diagrams illustrating MCG waveforms corresponding to increase of a number of studying when a magnetocardiogram is measured using the entirety of the MCG as a teacher signal; Figures 20(A) to (E) are diagrams illustrating MCG waveforms corresponding to increase of a number of studying times when a magnetocardiogram is measured using the magnetocardiogram measuring apparatus illustrated in Fig. 10; Figure 21 is a block diagram illustrating an arrangement of a spectrum analyzer which analyzes components of desired frequencies; Figures 22(A),(B),(C) are diagrams illustraing magnetic source analysis result obtained by a conventional method; Figure 23 is a diagram schematically illustrating an arrangement of a hierarchical perceptron; Figure 24 is a diagram illustrating a relationship between an error and a number of times of learning of the hierarchical perceptron; Figure 25 is a diagram illustrating varying of an error as a whole and varying of errors corresponding to individual pattern which corresponds to a number of times of learning of the hierarchical perceptron; and Figure 26 is a diagram schematically illustrating a Hopfield model. Best Modes for Executing The Invention Referring to the attached drawings, we explain the present invention in detail. Figure 1 is a block diagram illustrating a physical quantity analyzing apparatus according to an embodiment of the present invention. The physical quantity analyzing apparatus includes plural physical formula operating units 11, 12, ...1m, a sigma unit 2 for cumulatively adding operation results g1, g2, ...gm output from the physical formula operating units 11, 12, ...1m, an error operating section 3 for receiving the cumulatively added result Oj(t) output from the sigma unit 2 and a measured value Sj(t) of a physical quantity as a teacher pattern, and for calculating a difference between the both, correcting sections 11a, 12a, ...1ma for correcting estimated variables in the physical formula operating units based upon the calculated error, and an information collecting section 4 for collecting the estimated variables in the physical formula operating units 11, 12, ...1m and for outputting the estimated variables as analysis results therefrom. The physical formula operating units 11, 12, ...1m perform operation based upon physical formulae based upon known information in response to being supplied known information such as time t, measuring locations and the like, the physical formulae being determined in each physical formula operating unit, and correct variables included in the physical formulae in response to being supplied the estimated error dj (t) {= Sj (t) - Oj (t)} output from the error operating section 3 so as to decrease the estimated error. A number of physical formula operating units is determined to be greater than a number of causes of physical quantities which are objected for analysis. The physical formula operating units may be controlled to operate in synchronous manner to one another, or may be controlled to operate in asynchronous manner. Operations of the physical quantity analyzing apparatus having the above-mentioned arrangement are as follows. It is supposed that physical quantity Oj(t) which is objected for analysis is given by the following equation. That is, it is supposed that the physical quantity Oj(t) is a linear add of m-numbered functions gi(t, ai1, ai2, ai3, ...aiL) which have time t and L-numbered unknowns ai1, ai2, ai3, ...aiL. The functions gi may include known information in their formula. Such assumption never unnatural. Such assumption is effected when pressure, temperature, sound wave, electromagnetic wave, electric field, magnetic field, light wave, gravity, corpuscular ray and the like are objected for analysis. Under the assumption, functions g1, g2, ...gm are calculated by supplying known information such as a time t, measuring locations and the like, to the m-numbered physical formula operating units 11, 12, ...1m, a cumulatively added value Oj(t) can be obtained by supplying all calculated function values to the sigma unit 2. In this case, the obtained cumulatively added value Oj(t) is different from an actually measured value Sj(t) because the variables are determined arbitrarily. Therefore, a difference between the actually measured value Sj(t) and the cumulatively added value Oj(t) is calculated by the error operating section 3, the calculated difference is feedbacked as an estimated error dj(t) to the correcting sections 11a, 12a, ...1ma of the physical formula operating units 11, 12, ...1m, unknowns in the physical formula operating units are varied so as to decrease the estimated error dj(t). The estimated error dj(t) decreases following repetition of the above-mentioned series of processings, and finally becomes nearly equal to 0. Analysis results according to physical quantities of physical sources are obtained by collecting and outputting values of unknowns in the physical formula operating units 11, 12, ...1m by the information collecting unit 4 at the timing. When an estimated error evaluating function Ej(t) is defined by the following equation, the equation (4) is obtained. Ej (t) = (1/2) {Sj (t) - Oj (t)} 2 ∂Ej (t) / ∂Oj (t) = - {Sj (t) - Oj (t)} When it is supposed that correcting of unknowns in each physical formula operating unit is performed based upon a maximum slope decreasing method, estimation of unknowns for minimizing a value of the estimated error evaluating function can be performed based upon the equation (5). In the equation (5), εk represents a learning gain (correcting gain) of an unknown ai. aik = aik - εk {∂Ej (t) / ∂aik} = aik - εk {(∂Ej (t) / ∂Oj (t)} {∂Oj (t) / ∂aik} = aik + εk {Sj (t) - Oj (t)} {∂Oj (t) / ∂aik} The following equation is obtained from the equation (3). And, the following equation is obtained by substituting the equation (6) for the equation (5). aik= aik + εk {Sj (t) - Oj (t)} [∂ {gi (t, ai1, ai2, ai3, ··· aiL)} / ∂aik] Consequently, estimation accuracy of unknowns is improved by carrying out processing of the equation (7) so as to obtain more accurate physical quantities, when a physical model can be expressed by physical formulae gi(t, ai1, a2, ai3, ...aiL) and estimated physical quantities can be calculated by the equation (3). Figure 2 is a diagram illustrating an example of an estimated error evaluating function. Estimating operation of unknowns which minimize a value of the estimated error evaluating function is desccibed, together with Fig. 2 and a table 1 which indicates conditions of the estimated error evaluating function. In the chart 1, Δaik is a correction value according to an unknown aik. As being mentioned in the foregoing, it is sufficient that the estimating of unknowns is performed to decrease the value of the estimated error evaluating function. Therefore, the correction value Δaik is determined to be negative when a slope of the estimated error evaluating function is positive, and is determined to be positive when a slope of the estimated error evaluating function is negative, by taking a sign of the slope of the estimated error evaluating function Ej(t). It may be thought that an unknown aik corresponding to a local minimum is obtained when the above-mentioned estimation of unknowns is performed, because Fig. 2 includes a local minimum point. The estimation processing is not performed for only one unknown, but is performed for all unknowns in a synchronous manner, thereby the estimated error evaluating function itself is varied following repetition of the eatimation processing. As a result, unknowns which minimize the value of the estimated error evaluating function are obtained finally. Thereafter, analysis of physical quantities can be performed by collecting and outputting the finally obtained unknowns by the information collecting unit 4. Specific Example 1 Figure 3 is a schematic diagram illustrating an example of an immunity measuring apparatus using a fluorescence. When immunity measurement is carried out using the immunity measuring apparatus, a vessel 32 for reaction is provided at one side of a light guide 31, antigen 33 for example, is fixed at a boundary face of the vessel 32 for reaction with the light guide 31. Test solution 34 is injected into the vessel 32 for reaction so as to perform antigen-antibody reaction under the condition that exiting light is guided into the light guide 31 so as to propagate interior the light guide 31 in a total reflection manner, a quantity of antigen-antibody reaction corresponding to a degree of immunity. Then, labeled antibody 35 which is obtained by labeling antibody with fluorescence material, is injected into the vessel 32 for reaction. Only labeled antibody 35 contrained close to the boundary face is exited by evanescent wave component of the exiting light. The exited fluorescence propagates interior the light guide 31 in a total reflection manner, and outgoes from an exiting light incidenting face 31a of the light guide 31. The outgone fluorescence is separated by a beam splitter 36 and the like from the exiting light, and is guided to an optical sensor 37 such as a photomultiplier so as to obtain a measured signal. Figure 4 is a schematic diagram illustrating varying of a measured signal following passage of time, the measured signal being obtained by the above-mentioned immunity measuring appatarus. An offset caused by fluorescence and the like generated by the light guide 31 itself, is obtained as the measured signal till the labeled antibody 35 is injected into the vessel 32 for reaction. After the labeled antibody 35 is injected into the vessel 32 for reaction, the labeled antibody 35 performs antigen-antibody reaction with the antigen 38 which has performed the antigen-antibody reaction so that the labeled antibody 35 is contrained close to the boundary face. And the quantity of the contarined labeled antibody 35 gradually increases, the measured signal increases like an exponential function with respect to the offset, and finally becomes a predetermined value which is determined based upon the degree of immunity reaction. Therefore, the degree of immunity cannot be measured accurately when not only the finally obtained measured signal but also the offset are not obtained accurately. By taking the point into consideration, it can be thought that the offset is estimated using a least square method (primary regression) of linear approximation. An error causes greatly in estimation of the offset depending upon the degree of immunity, so that immunity measurement cannot be performed with high accuracy. But, the quantity of immunity reaction can be estimated with high accuracy by using the physical source analyzing apparatus having the arrangement in Fig. 1. That is, it is experimentally recognized that measurement signal corresponding to outgone fluorescence can be calculated based upon the following equation. In the equation, A indicates a sum of immunity fluorescece and non-specific adsorption fluorescence, B indicates a value depending to density of immunity antibody, C indicates a beginning time of reaction, and D indicates an offset of the light guide 31. gj (t) = A {1 - e -B (t-C) } + D In this case, a number of physical formula operating units is 4 sufficiently, because a number of unknowns to be estimated is 4. Each physical formula operating unit is determined to perform processing based upon the equation (8). Each unknowns A, B, C, and D is given arbitrary initial value, respectively. The above-mentioned unknowns A, B, C and D are estimated with high accuracy by carrying out a processing based upon a flowchart illustrated in Fig. 5, when time is supplied to the physical source analyzing apparatus which is initially determined, as common measurement condition, and a measured signal which is obtained according to the supplied time is supplied as a teacher. signal. That is, degree of immunity reaction can be estimated with high accuracy. More particularly, partial differential values of the measured signal are given by the following equatuins, because it is experimentally known that the measured signal can be calculated by the equation (8). As a result, estimating rules of unknowns A, B, C and D are as follows. Therefore, in step SP1, the unknowns A, B, C and D are given initial values, in step SP2, an arbitrary time t after beginning immunity reaction, and the measured signal Sj(t) corresponding to the time t , are extracted, in step SP3, operations are carried out by physical formula operating units 11, 12, 13 and 14 based upon the estimated unknowns A, B, C and D at the present time, and cumulatively adding is carried out by the sigma unit 2, in step SP4, a difference Sj (t) - Oj (t) between the measured signal and the cumulatively added value is calculated, in step SP5, it is judged whether or not the calculated difference is equal or lesser than a desired estimated error which was previously determined. When it is judged that the calculated difference is greater than the desired estimated error, in step SP6, partial differential values of the cumulatively added value are calculated by carrying out the processing of equations (9), in step SP7, the unknowns A, B, C and D are estimated by carrying out the processing of equations (10), then processing in step SP2 is carried out again. When it is judged in step SP5 that the calculated difference is equal or lesser than the desired estimated error which was previously determined, the series of processing is finished as it is. The unknowns A, B, C and D can be estimated with high accuracy by carrying out the above-mentioned series of processing. The unknowns can be estimated with high accuracy by using initial data. A required time period can be extremely shortened because processing for estimating is simplified. Estimating unknowns in chemical reaction based upon different operating formulae, experimental formulae can be dealt with varying data processing part, that is varying software program. Of course, other parameters can be estimated. Specific Example 2 Figure 6 is a schematic view illustrating a glucose concentration measuring apparatus. When measuring of concentration of glucose is carried out using the glucose concentration measuring apparatus, for example, a hydrogen peroxide selectively penetration membrane 52, a glucoseoxidaze (hereinafter referred to as GOD) immobilized membrane 53 which is obtained by immobilizing GOD, and a diffusion limiting membrane 54 which limits diffusion of glucose to some degree, are laminated on a surface of a base electrode 51 which includes a platinum electrode 51a and a silver electrode 51b. Test solution is dropped onto the diffusion limiting membrane 54. Glucose in the dropped test solution is limited its diffusion to some degree by the diffusion limiting membrane 54, and glucose with lower concentration than the glucose concentration in the test solution, reaches to the GOD immobilized membrane 53. The following reaction is carried out under the existence of GOD. Hydrogen peroxide generated by the reaction passes through the hydrogen peroxide selectively penetration membrane 52, and is guided to the surface of the base electrode 51. A measured signal corresponding to the quantity of the hydrogen peroxide is output from the base electrode 51. Consequently, concentration of glucose can be measured based upon the measured signal. As is apparent from the reaction formula, quantity of glucose which can react, that is measurable limit of concentration of glucose, is determined by the quantity of oxigen in the test solution, therefore the measurable limit is raised by limiting diffused quantity of glucose by interposing the diffusion limiting membrane 54 instead of directly dropping the test solution onto the GOD immobilized membrane 53. Figure 7 is a schematic diagram illustrating varying of a measured signal following passage of time, the measured signal being obtained by the above-memtioned glucose concentration measuring appatarus. The measured signal gradually decreases by being influenced by the preceding measurement till the test solution is dropped. After the test solution is dropped, the measured signal increases like an exponential function because the reaction of the above-mentioned reaction formula is carried out, and the measured signal finally becomes a predetermined value which is determined based upon the concentration of glucose. Fig. 7 illustrates a condition in which quantity of oxygen in the test solution is sufficient. It is general that the concentration of glucose is measured by obtaining a maximum value of differential values of the measured signal by time, because a required time period becomes extremely long when the concentration of glucose is measured by waiting till the condition in which the measured signal doe's not vary. Furthermore, the obtained measured signal is influenced by not only the concentration of glucose, but also the thickness, laminated condition and the like of membranes laminated on the base electrode 51, it is general that a calibration is carried out using solution, glucose concentration of which is known, at every exchanging the membrane. But, the concentration of glucose can be estimated with high accuracy by using the physical source analyzing apparatus having the arrangement in Fig. 1. That is, it is recognized that the measured signal can be calculated based upon the following equation which can draw a second diffusion equation. In the following equation, A indicates a product of the concentration of glucose and a quantity determined by the quality of the membrane, B indicates a numeric value corresponding to raising speed of the reaction, C indicates a beginning time of the reaction, and D indicates a value of measured signal at the beginning time of the reaction. It is known that B becomes great when the thickness of the membranes is thin, and that B becomes great when the concentration of glucose is high, and that the reaction converges rapidly when B is great. In this case, a number of physical formula operating units is 4 sufficiently, because a number of unknowns to be estimated is 4. Each physical formula operating unit is determined to perform processing based upon the equation (11). Each unknowns A, B, C and D are given arbitrary initial values, respectively. The above-mentioned unknowns A, B, C and D are estimated with high accuracy by carrying out a processing based upon the flowchart in Fig. 5, when time is supplied to the physical source analyzing apparatus which is initially determined, as common measurement condition, and a measured signal which is obtained according to the supplied time is supplied as a teacher signal. That is, concentration of glucose can be estimated with high accuracy. More particularly, partial differential values of the measured signal are given by the following equations, because the measured signal can be calculated by the equation (11). As a result, estimating rules of unknowns A, B, C and D are as follows. A = A + εa {Sj (t) - Oj (t)} (∂Oj (t) / ∂A) B = B + εb {Sj (t) - Oj (t)} (∂Oj (t) / ∂B) C = C + εc {Sj (t) - Oj (t)} (∂Oj (t) / ∂C) D = D + εd (Sj (t) - Oj (t)} (∂Oj (t) / ∂D) Therefore, in step SP1, the unknowns A, B, C and D are given initial values, in step SP2, an arbitrary time t after beginning the reaction according to the reaction equation, and the measured signal Sj(t) corresponding to the time t, are extracted, in step SP3, operations are carried out by physical formula operating units 11, 12, 13 and 14 based upon the estimated unknowns A, B, C and D at the present time, and cumulatively adding is carried out by the sigma unit 2, in step SP4, a difference Sj(t)-Oj(t) between the measured signal and the cumulatively added value is calculated, in step SP5, it is judged whether or not the calculated difference is equal or lesser than a desired estimated error which was previously determined. When it is judged that the calculated difference is greater than the desired estimated error, in step SP6, partial differential values of the cumulatively added value are calculated by carrying out the processing of equations (12), in step SP7, the unknowns A, B, C and D are estimated by carrying out the processing of equations (13), then processing in step SP2 is carried out again. When it is judged in step SP5 that the calculated difference is equal or lesser than the desired estimated error, the series of processing is finished as it is. The unknowns A, B, C and D can be estimated with high accuracy by carrying out the above-mentioned series of processing. The unknowns can be estimated with high accuracy by using initial data. A required time period can be extremely shortened because processing for estimating is simplified. Estimating unknowns in chemical reaction based upon different operating formulae, experimental formulae can be dealt with varying data processing part, that is varying software program. Of course, other parameters can be estimated. Furthermore, it is possible that the unknown A is classified into an unknown Al which corresponds to concentration of glucose and an unknown A2 which corresponds to a quantity which is determined based upon quality of membranes. In this case, calibration becomes unnecessary because concentration of glucose can directly be estimated. In the foregoing, only immunity measurement and concentration measurement of glucose are described, it is of course possible that the apparatus can be applied to estimating various parameters in other chemical reactions. Second Embodiment Figure 8 is a schematic block diagram illustrating an arrangement of a magnetic source analyzing apparatus using the apparatus in Fig. 1. The apparatus differes from the apparatus of Fig. 1 in that N-numbered magnetic field sensors MS are disposed so as to detect magnetic field, and one measured value among N-numbered measured values and its observation condition are selected by multiplexers MX1 and MX2, respectively, so as to supply the selected measured value and the observation condition to the error operating section 3 as teacher signals, the multiplexers MX1 and MX2 being controlled by a control circuitry C. When a measured value obtained by the j-th magnetic field sensor j among the N-numbered magnetic field sensors, is determined to be Bj, and when an estimated value at the measuring point j by the magnetic field sensor j is determined to be Bej, an estimated error Ej at the measuring point j by the j-th magnetic field sensor j can be expressed by the following equation. Ej = (1/2) (Bj -Bej) 2 When Bej is formed by a magnetic field which is generated by m-numbered current elements, an estimated value Bej can be obtained based upon location information and a current vector of each current element i, and location information of the measuring point j. When a magnetic field generated at the measuring point j by the current element i is expressed as Beji, the equation (14) is effected. That is, linear addition is satisfied. When input parameters of the current element i are indicated with Wik (k=1, 2, ...L) {for example, input parameters become Wi1=xi, Wi2=yi, Wi3=zi, Wi4=Xi, Wi5=Yi, Wi6=Zi, when the input parameters are location information Oi(xi, yi, zi) and a current vector Ii(Xi, Yi, Zi)}, and when varying component of the input parameters Wik is indicated with AWik, the equation (15) should be satisfied to guide the estimated error Ej to the minimum using a maximum slope decreasing method. ΔWik∝ - ∂Ej / ∂Wik By taking these relationships into consideration, the equation (16) is obtained as an operation for correcting the input parameters Wik so as to decrease the estimated error Ej. Wik = Wik - εk · ∂Ej / ∂Wik = Wik - εk (∂Ej / ∂Bej) (∂Bej / ∂Beji) (∂Beji / ∂Wik) = Wik + εk (Bj - Bej) (∂Beji / ∂Wik) In the equation (16), εk indicates a learning gain, and i=1, 2, 3, ...m, and k=1, 2, 3, ...L. When a measured value obtained by the j-th magnetic field sensor j among the N-numbered magnetic field sensors, is determined to be Bzj, and when a magnetic field component in z-direction is determined to be Bzeji, the magnetic field component in z-direction being generated at a measuring point j by a current element i which have a current vector parallel to an x-y plane, it is known that the magnetic field component Bzeji is given by the following equation. Bzeji = (µO · Mi / 4π) {(yj - yi) cosi - (xj - xi) sini} / {(xj - xi) 2 + (yj - yi) 2 + (zj - zi) 2 } 3/2 In the equation, the unknowns xi, yi and zi respectively indicate x-, y-, and z-coordinates of the current element i, the unknown i indicates an angle of the current element i with respect to the x-axis on the x-y plane, the unknown Mi indicates a moment of the current element i, and the unknowns Xj, yj and zj respectively indicate x-, y-, and z-coordinates of the measuring point j. Consequently, analyzing of a magnetic source based upon the maximum slope decreasing method can be performed by substituting each unknown for the equation (16). That is, the equations (17) are obtained for unknowns. In the equations (17), x, y, z, T and M are very small positive values, and are learning gains. Correcting processing of unknowns based upon the equations (17) in the physical formula operating units are determined to be carried out in synchronous manner. Correcting processing of the unknowns based upon the equations (17) can be performed easily, because each partial differential term included in the equations (17) can be obtained analytically based upon the equations (18). In the equation (18), X = xj - xi, Y = yj - yi, Z = zj - zi, A = X 2 + Y 2 + Z 2 for the purpose of simplification of the equations. The unknowns with high accuracy can be obtained with relatively small number of learning times by performing magnetic source analysis using the apparatus in Fig. 8. The obtained unknowns can be collected and output as analyzed results by the information collecting unit 4. Though the learning gains for correcting the unknowns are determined to be sufficiently small value, a necessary number of learning times for convergence can be small because each unknown can be greatly corrected when the difference between the measured value Bzj and the estimated value Bzej is great. When the difference is small, each unknown can be corrected by an extremely small quantity, thereby accuracy of unknowns which are finally obtained, can be raised. When the difference between the measured value Bzj and the estimated value Bzej is small, and when unknowns are local minimum (not true solution), each unknown can be corrected so as to escape the local minimum condition because at least one partial differential term in the equation (17) has a finite value. The estimated error increases temporarily, but each unknown can be corrected so as to decrease the error again. Therefore, performing simulated annealing which is employed in other system, is not needed. The embodiment describes only a case in which the apparatus is applied to analysis for a magnetic source, but the apparatus can be applied in the similar manner for a system in which physical sources radiate some physical quantities (heat, electromagnetic wave and the like). Third Embodiment Figure 9 is a schematic block diagram illustrating an arrangement of an echo searching apparatus using the apparatus in Fig. 1. The apparatus differes from the apparatus of Fig. 1 in that sound wave is radiated from one wave radiating device 7, m-numbered wave receiving sensors i for receiving reflected waves from n-numbered reflecting points, are disposed so as to detect the reflected waves, and one measured value among n-numbered measured values and its observation condition are selected by multiplexers MX1 and MX2, respectively, so as to supply the selected measured value and the observation condition to the error operating section 3 as teacher signals, the multiplexers MX1 and MX2 being controlled by a control circuitry C. When it is supposed that the wave radiating device 7 is an extremely small point sound source, and radiates a short pulse with a time waveform of S(t), and a velocity potential is expressed by ijk, the velocity potential being generated at time j by arriving the sound wave reflected by the reflecting points k (k=1, 2, ...n) to the wave reseiving sensors i (i=1, 2, ...m), as is illustrated in Fig. 10, it is known that velocity potential ij observed by the wave receiving sensor i at time j is determined by the following equation, when sound pressure is sufficiently small. Passages of sound wave which generates velocity potential at the wave receiving sensors i, are common from the wave radiating device 7 to the reflecting point k, and are different from the reflecting point k to each wave receiving sensor i from one another. When the distance from the wave radiating device 7 to the reflecting point k is indicated with rsk, and the distance from the reflecting point k to the wave receiving sensor i is indicated with rRki, the following equation can be obtained, because velocity potential is damped in reverse proportion to distance. ijk = αk · S {t - (rsk + rRki)/c} · (1 / rsk) · (1 / rRki) In the equation, c indicates sound speed, and αk indicates a constant which proportions to the reflecting rate of the reflecting point k. Conversion between time t and time j is supposed to be performed by t=Δt·j based upon a sampling interval Δt. Furthermore, coordinates of the wave radiating device 7 is supposed to be (xs, ys, zs), coordinates of the reflecting point k is supposed to be (xk, yk, zk), and coordinates of the wave receiving sensor i is supposed to be (xi, yi, zi), for the purpose of uniting coordinates, the following equations are obtained. rsk = {(xs - xk) 2 + (ys - yk) 2 + (zs - zk) 2 } 1/2 rRki = {(xi - xk) 2 + (yi - yk) 2 + (zi - zk) 2 } 1/2 ijk can be expressed with a function {ijk =  (αk, xk, yk, zk)} which has four unknowns αk, xk, yk and zk. Therefore it is possible to apply to the apparatus in Fig. 9. The input pattern becomes j, xi, yi, zi, and the teacher pattern becomes a measured value ij of ij. As is apparent from the foregoing, echo searching with high accuracy can be performed with relatively small number of learning times, by determining studying rules of the unknowns in the physical formula operating units in Fig. 9, as the following equations. In these equations, εa, εx, εy, εz are learning gains and are positive constants. Partial differential terms in the equations (19) can be calculated by performing formula deformation arithmetically, and also can be calculated based upon numerical operation method. When the partial differential term is calculated based upon the numerical operation method, the partial differential term can easily be calculated by the following equation. ∂ijk/∂αk = { (αk + εa, xk, yk, zk) -  (αk - εa, xk, yk, z f )}/2εa The embodiment describes only a case in which the apparatus is applied to echo searching, but the apparatus can be applied in the similar manner for performing analysis (analysis for absorbance characteristic of infrared and the like) based upon physical quantities which is not radiated from physical sources. Fourth Embodiment Figure 11 is a block diagram illustrating an apparatus for reducing line spectrum noise for reducing only noises having line spectrum characteristics from a signal objected for measurement, when a frequency band of the signal and a frequency band of the noises overlap to one another. The apparatus differs from the apparatus of Fig. 1 in that the information collecting unit 4 is omitted. Such noises having line spectrum characteristics are reduced for raising analysis accuracy of data processing in a latter stage. In convention, an analoque filter or an adaptive filter is generally employed. A band elimination filter, one species of analogue filters, is simplified its arrangement and is cheap. But, the band elimination filter has disadvantages that the filter attenuates not only line spectrum noise components of a noise source but also spectrum components of a signal objected for measurement, and that the filter shifts not only amplitude but also phase in neighbouring portion to a center frequency, and that the filter inevitably lowers analysis accuracy of the signal objected for measurement. When the adaptive filter is employed, only line spectrum noise components of a noise source can be reduced, and analysis accuracy of a signal objected for measurement can be raised. But, the adaptive filter has disadvantage that the arrangement of the adaptive filter becomes complicated because a reference channel is necessary, the reference channel measuring only line spectrum noise components of a noise source. The adaptive filter cannot be applied to a case in which measuring only line spectrum noise components of a noise source is impossible. The apparatus for reducing line spectrum noise illustrated in Fig. 11 distinguishes the above-mentioned disadvantages. The apparatus supplies operation results output from physical formula operating units 11, 12, ...ln to the sigma unit 2, a number of physical formula operating units being equal to a number of line spectrum noises which should be reduced. The apparatus supplies the cumulatively added result O(t) output from the sigma unit 2, and the measured value S(t) of physical quantity as the teacher pattern, to the error operating section 3 so as to calculate the difference d(t) between the both. The apparatus feedbacks the calculated difference d(t) to the correcting sections 11a, 12a, ...1ma of the physical formula operating units 11, 12, ...1m, and outputs the difference d(t) as a signal which is reduced line spectrum noises therefrom. The number and frequency of line spectrum noises objected for reducing are previously measured. In the correcting sections 11a, 12a, ...1ma, unknowns included in physical formula are corrected so as to minimize the square value of the error d(t). Operations of the apparatus for reducing line spectrum noise having the arrangement mentioned above, are as follows. When it is supposed that frequencies fi of line spectrum noises which are causes of interfering, are measured by a measuring means which is conventionally known, corresponding line spectrum can be expressed by the following equation. g (t, Ai, i) = Ai · sin (2π fit + i) Therefore, the cumulatively added result O(t) can be expressed by the following equation. The difference d(t) becomes d (t) = S (t) - O (t). And, the partial differential values of the cumulatively added result O(t) becomes the following equations. ∂O (t) / ∂Ai = sin (2π fit + i) ∂O (t) / ∂i = Ai · cos (2π fit + i) From the foregoing equations, equations for correcting unknowns Ai and i become the following equations. Ai = Ai + εA · d (t) · sin (2π fit + i) i = i + εT · d (t) · Ai · cos (2π fit + i) In the equations, εA and εT are positive constants, respectively. Therefore, the unknowns Ai and i can be estimated with high accuracy by performing operations based upon the above-mentioned equations for correction based upon the difference d(t) in the correcting sections 11a, 12a, ...1ma of the physical formula operating units 11, 12, ...1m. When the unknowns Ai and i are accurately estimated, corresponding line spectrum noises are output from the physical formula operating units 11, 12, ...1m. The line spectrum noises are cumulatively added by the cumulative addition device 2, and are subtracted from the measured signal S(t) by the error operating section 3 after being cumulatively added by the cumulative addition device 2. Thereby, the difference d(t) is output which is obtained by reducing only line spectrum noises from the measured signal. The following software program is a software program for estimating an amplitude Ai and a phase i of one line spectrum noise. In the software program, it is defined that FNF (j, A, f, ) = A * sin (2π ·* f * j / 1000 + ), FNG (j, A, f, ) = A * cos (2π * f * j /1000 + ), εA = 0.003, εT = 0.000003 In the above-mentioned software program, when the amplitude of the line spectrum is determined to be 48.0000, the phase of the line spectrum noise is determined to be 0.9767, and the frequency of the line spectrum noise is determined to be 90.000000, and simulations are carried out by giving various initial values, the same estimated values {the amplitude Ai is 48.0006 or 47.9994, the phase i is 0.9767, the estimated error dj is ±SE-4 degree, the value of the estimated function Ej is 3E-7 degree} as the determined values are obtained as is illustrated in Fig. 12, by repeating the processing by total 10,000 times. Ai seems to have two answers, but it is caused that the calculation is carried out with 32 bit single precision floating point form so that accuracy of 5.5 digits degree at decimal number is obtained. When the amplitude of the line spectrum is determined to be 72.0000, the phase of the line spectrum noise is determined to be 1.3525, and the frequency of the line spectrum noise is determined to be 18.000000, and simulations are carried out by giving various initial values, the same estimated values {the amplitude Ai is 71.9987, the phase i is 1.3525, the estimated error dj is 1.2E-3 degree, the value of the estimated function Ej is 1.53E-6 degree} as the determined values are obtained as is illustrated in Fig. 12, by repeating the processing by total 10,000 times. That is, accuracy of 5.5 digits degree at decimal number is obtained. The following software program is a software program for estimating amplitudes Ai and phases i of plural line spectrum noises. In the software program, it is defined that FNF (j, A, f, ) = A * sin (2π * f * j / 1000 + ), FNFA (j, A, f, ) = sin (2π * f * j / 1000 + ), FNFT (j, A, f, ) = A * cos (2π * f * j / 1000 + ), εA = 0.005, εT = 0.000005, k0, k1 and that k0 and k1 are unknown numbers. In the above-mentioned software program, when it is detrermined that k0=7, and k1=7, and amplitudes, phases and frequencies of 7 line spectrum noises are determined to be 69.0000, 2.1046, 34.000000; 75.0000, 2.6461, 29.000000; 86.0000, 5.0893, 45.0000; 29.0000, 5.9306, 67.000000; 90.0000, 0.8073, 20.000000; 37.0000, 3.4211, 90.000000; 77.0000, 1.6288, 46.000000, and simulations are carried out by giving various initial values, amplitudes and phases corresponding to each line spectrum are estimated as 69.0000, 2.1046; 75.0003, 2.6461; 85.9998, 5.0894; 29.0000, 5.9308; 90.0000, 0.8073; 36.9999, 3.4211; 77.0002, 1.6288. In these cases, the estimating processing number of times is 10,000, the estimated error dj is -4.8828E-4, the estimated function Ej is 2.38419E-7. That is, each unknowns are estimated with an accuracy of 6 digits portion which is a limit of 32 bit single precision floating point form. The estimating processes of these 7 line spectrums are illustrated in Fig. 14. In this embodiment, when the cumulatively added value O(t) is defined by the following equation, direct current offset can also be cancelled. In this case, it is necessary that a physical formula operating unit 10 for performing the operation of AO = AO + εA · d (t) and for estimating unknowns, is added. Fifth Embodiment Figure 16 is a block diagram illustrating a magnetocardiogram measuring apparatus. The apparatus differs from the apparatus of Fig. 11 in that the apparatus of Fig. 16 further incluses a magnetocardiogram measuring section 61, an electrocardiogram measuring section 62, an R-wave extracting section 63 for extracting an R-wave based upon the obtained electrocardiogram (hereinagter referred to as ECG), a time extent outputting section 64 for obtaining time extents of a P-wave and T-wave which are prior and sebsequent with respect to the extracted R-wave, a triggering section for triggering the obtained magnetocardiogram (hereinafter referred to as MCG) corresponding to the R-wave of the ECG, and an magnetocardiogram limiting section 66 for inhibiting the MCG to be supplied to the error operating section 3, the inhiited MCG corresponding to the time extent of the P-wave and the T-wave. The R-wave extracting section 63 obtains cross correlations of a Q-R-S group using a rectangular window (refer to Fig. 17) which comprises a first negative wave, a positive wave and a second negative wave, which waves are determined to a predetermined time ratio (t1 : t2 : t3), and determines a time corresponding to the maximum value of the cross correlations, as a center time of the Q-R-S group, that is a time of the R-wave, by paying attention to the fact that a shape of the Q-R-S group in the ECG is negative-positive-negative pulse wave. The magnetocardiogram measuring section 61, electrocardiogram measuring section 62, and triggering section 65 are omitted their detailed descriptions because conventionally known apparatus can be employed as the sections. Operations of the magnetocardiogram measuring apparatus having the above-mentioned arrangement are as follows. The inventor has continued studying for magnetocardiogram earnestly, and has found that spectrums having the same frequencies as of frequency components of line spectrums which is caused by a commercial power source, noises of a refrigerator and the like, appear at appearing time of the R-wave and T-wave of the ECG, and that the appeared spectrums act as disturbances to an estimating processing of line spectrum noises in a magnetocardiogram measurement. The inventor has also found that a required time period for the estimating processing is extremely lengthened caused by the above-mentioned spectrums which appear at the appearing time of the R-wave and T-wave, when estimating processings of unknowns in each physical formula operating unit are performed by supplying the entirety of the magnetocardiogram to the error operating section 3 as a teacher signal. The embodiment was made based upon the above-mentioned knowledge. The embodiment improves convergence by omitting samplings of magnetocardiogram corresponding to a region from the beginning of the P-wave to the end of the T-wave, from estimating processing of unknowns, and extremely shortens the required time period for the estimating processing. More particularly, a relationship between the MCG and the ECG in ideal condition in which no spectrums as disturbances are included, is illustrated in figs. 18 (A) and (B). And a waveform which is obtained by mixing spectrum noises of 29Hz and 60 Hz to the MCG in ideal condition, is illustrated in Fig. 18 (C). When unknowns are estimated by employing the entirety of the waveform of Figs. 18 as teacher signals, using the line spectrum noise removing apparatus having the arrangement of Fig. 11, spectrum noises are greatly removed by carrying out the estimating processing by 1, 2, 3, 4 and 5 times based upon 2,500 samples of the teacher signal which is illustrated in Fig. 19 (A) (refer to Figs. 19 (B), (C), (D), (E) and (F)), and the MCG waveform (refer to Fig. 19 (F)) fairly similar to the waveform of Fig. 19 (A) is obtained. But, spectrum noises are remained to some degree in the waveform of Fig. 19 (F), and become a cause for lowering accuracy in data processings of subsequent stages. Therefore, a number of times of estimating processings should be increased when the spectrum noises are further decreased. As a result, the required time period is lengthened greatly. On the contrary, when unknowns are estimated using the magnetocardiogram measuring apparatus having the arrangement of Fig. 16, a region (refer to the region R1 in Fig. 18 (B)) which includes from the P-wave to the T-wave of the ECG illustrated in Fig. 18 (B), is extracted by the R-wave extrating section 63, and a time extent from the beginning of the P-wave to the end of the T-wave is obtained by the time extent outputting section 64 based upon the peak time of the R-wave. And the MCG obtained by the magnetocardiogram measuring section 61 is triggered corresponding to the R-wave of the ECG by the triggering section 65, the MCG corresponding to the time extent of P-wave and T-wave, is inhibited by the magnetocardiogram limiting section 66 to be supplied to the error operating section 3. When unknowns are estimated using the magnetocardiogram in the extent which is limited as such, as the teacher signal, spectrum noises are greatly removed by carrying out the estimating processing by 1, 2, 3, 4 and 5 times based upon 2,500 samples within the limited extent (refer to Figs. 20 (A), (B), (C), (D) and (E)), and the MCG waveform (refer to Fig. 20 (E)) fairly similar to the waveform of Fig. 19 (A) is obtained. When the waveforms in Figs. 19 and the waveforms in Figs. 20 are compared, it is understood that the waveform in Fig. 20 (A) is similar to an intermediate waveform (a result of estimating processing by 4.5 times) between the waveform in Fig. 19 (E) and the waveform in Fig. 19 (F). As is apparent from the result, a required time period for estimating processing can extremely be shortened because a number of times of estimating processings can extremely be decreased, when it is sufficient that spectrum noise removal as similar degree as the spectrum noise removal by the line spectrum noise removing apparatus having the arrangement of Fig. 11. When the estimating processings are carried out by a number of times which is the same as of the line spectrum noise removing apparatus having the arrangement of Fig. 11, spectrum noise removing efficiency can extremely be raised, and accuracy in data processing in subsequent stages can extremely be raised. Furthermore, operability can extremely be raised because operation for preparing a template by collecting only spectrum noises, becomes unnecessary. Sixth Embodiment Figure 21 is a block diagram illustrating an arrangement of a spectrum analyzer which analyzes components of desired frequencies. The analyzer supplies operated results output from physical formula operating units 11, 12, ...1m, to a sigma unit 2 so as to operate cumulative addition, a number of the physical formula operating units corresponding to a number of frequencies objected for analysis, supplies the cumulatively added result O(t) and a measured data S(t) objected for analysis to an error operating section 3 so as to calculate an estimated error d(t), feedbacks the estimated error d(t) as information for correction of unknowns, to correcting sections 11a, 12a, ...1ma of the physical formula operating units 11, 12, ...1m. When estimating processing for unkowns converges, converged unknowns are extracted from the physical formula operating units 11, 12, ...1m, and are output as frequency component analysis result, by an information collecting unit 4. In the correcting sectioins 11a, 12a, ...1ma, unknowns included in physical formulae are corrected so as to minimize a square value of the error d(t), for example. From past years, fast Fourier transformation (hereinafter referred to as FFT) is known as a method which can analyze frequency components with high accuracy. But, the FFT has disadvantages that a low pass filter called antialiasing filter is necessarily provided for removing unnecessary higher harmonics which are usually included in a measurement signal, and an arrangement of FFT is complicated, because the FFT is a method for performing processing based upon sampling theorem. The FFT has another disadvantage that species of applicable signals are limited, because accuracy of frequency analysis result is extremely lowered when it is not assured that data continue periodically within a sampling interval. To dissolve these disadvantages, it is proposed to use a window function such as hamming, hanning and the like. A new disadvantage arises in that a waveform is distorted after reverse filtering operation in the entirety of the sampling interval, when the window function is used in the reverse filter operation, and analysis accuracy is lowered. A disadvantage also arises in that a sampling number should inevitably be increased when analysis for a wide band is required, because only outputs at frequency intervals are obtained, which are determined based upon sampling intervals. Disadvantages further arise in that an extremely large amount of memory is needed because sampling must be performed at regular intervals even when an axis indicating frequency is expressed by logarithmic scale, and that a sampling number is constrained to be 2 N . This embodiment was made to dissolve these disadvantages. The embodiment makes the antialiasing filter unnecessary, and can use sampling data at irregular intervals, and can determine analyzing frequencies arbitrarily. More particularly, the following equation can be used as a model of a component having frequency fi in a signal objected for analysis. In the equation, Ai and Bi are constants indicating amplitudes, respectively. Relationships of the following equations are obtained from the above-mentioned equation. The following equations are obtained by substituting these equations (20) for the equation (5). d(t) = S(t) - O(t) Ai = Ai + εA · d (t) · cos (2π fit) Bi = Bi + εB · d (t) · sin (2π fit) C = C + ε C · d (t) Therefore, the operation of O(t) is carried out in the sigma unit 2, the operation for calculating the estimated error d(t) in the error operating section 3 is carried out, the operations for correcting the unknowns Ai, Bi and C are carried out based upon the estimated error d(t) and the operations of physical formulae are carried out based upon the corrected unknowns, in each physical operating unit. When the estimation of the unknowns Ai, Bi and C is converged, each converged unknown is collected by the information collecting unit 4, thereby analysis of a component of arbitrary frequency included in a signal objected for analysis, can be performed with high accuracy. Possibility of Industrial Utilization As is apparent from the foregoing, methods and apparatus for analyzing physical source according to the present invention, can calculate physical quantity of a physical source with high accuracy and with extremely short time period, in comparison with conventional analyzing apparatus which use a supercomputer and the like, based upon measured values of physical quantity, which values being obtained at plural points apart from the physical source such as a magnetic field source interior a living body and others, and observation condition at the measuring timing. The methods and apparatus can preferably be applied especially to an apparatus for measuring biological magnetic field source, and an apparatus for reducing line spectrum noise.

Claims: A method for analyzing a physical quantity of each physical source based upon measured physical quantities (Sj(t)) measured at predetermined points apart from the physical sources, in case that measurable physical quantities measurable at arbitrary points apart from an individual physical source can be calculated based upon predetermined operation formulae which include physical quantities of physical sources and observation conditions, and measurable physical quantities measurable at the arbitrary points apart from plural physical sources satisfy linear addition, characterised by the steps of; operating plural physical formulae based upon known information, the physical formulae being determined corresponding to species of physical quantities and including variables to be corrected in a subsequent step, calculating a difference (dj(t)) between a measured physical quantity (Sj(t)) and a value (Oj(t)) obtained by cumulatively adding operation results (gi) of the physical formulae, correcting plural variables (ai1, ai2, ...aiL) based upon calculated difference (dj(t)), the variables (ai1, ai2, ...aiL) being included in each physical formula, and outputting corrected variables included in each physical formula as analysis results of physical quantities, after repeating the series of processing till the difference (dj(t)) becomes sufficiently small. An apparatus for analyzing a physical quantity of each physical source based upon measured physical quantities (Sj(t)) measured at predetermined points apart from the physical sources, in case that measurable physical quantities measurable at arbitrary points apart from the individual physical source can be calculated based upon predetermined operation formulae which include physical quantities of physical sources and observation conditions, and measurable physical quantities measurable at the arbitrary points apart from plural physical sources satisfy linear addition, characterised by; physical formula operating means (11, 12, ..., 1m) for operating based upon the operation formulae, whereby a number of physical formula operating means is more than a number of causes of physical quantities which are objected for analysis, cumulative addition means (2) for cumulatively adding operation results (gi) output from each physical formula operating means (11, 12, ..., 1m), error calculating means (3) for receiving acumulatively added operation result (Oj(t)) output from the cumulative addition means (2) and a measured value (Sj(t)) of physical quantity, and for calculating an error, correcting means (11a, 12a, ...,1ma) for correcting causes of the physical quantity in each physical formula operating means (11, 12, ..., 1m) based upon the calculated error, and corrected result collecting means (4) for collecting results which are corrected by the correcting means (11a, 12a, ..., 1ma), and for outputting the collected results as analysis results of physical quantities of physical sources. An apparatus for analyzing physical quantities as set forth in claim 2, wherein the physical quantity measurable at arbitrary points apart from individual physical sources is a quantity of fluorescence labeled material (35) which are constrained to the neighbouring to a light guide (31) by antigen-antibody reaction, the predetermined formula is an experimental formula regulating intensity of fluorescence output from the light guide (31) corresponding to introducing an exiting light in the light guide (31) which exiting light propagates in a total reflection manner, the number of physical formula operating means (11, 12, ..., 1m) is a predetermined number which is not lesser than a number of unknowns included in the experimental formula, and the corrected result collecting means (4) outputs corrected unknowns as immunity measured results which corrected unknowns correspond at least to immunity fluorescence and non-specific adsorption fluorescence. An apparatus for analyzing physical quantities as set forth in claim 2, wherein the physical quantity measurable at arbitrary points apart from individual physical sources is a quantity of material which are generated or consumed under the existence of enzyme, the predetermined formula is a formula regulating strength of electric signal output from a base electrode (51), which support an enzyme-immobilized membrane (53), corresponding to depositing a test solution on the enzyme-immobilized membrane (53) directly or indirectly, the number of physical formula operating means (11, 12, ..., 1m) is a predetermined number which is not lesser than a number of unknowns included in the formula, and the corrected result collecting means (4) outputs corrected unknowns as measured concentration results of corresponding material which corrected unknowns correspond at least to concentration of material which perform reaction under the existence of the enzyme. An apparatus for analyzing physical quantities as set forth in claim 2 wherein the physical quantity to be analyzed being physical characteristics of the physical source, the apparatus further comprising; wave radiating means (7) for radiating a wave to physical sources, wave receiving means (i) for receiving reflected waves from the physical sources, and for obtaining measurement value of physical quantity, wherein Said plural physical formula operating means (11, 12, ..., 1m) are provided for operating physical formulae which correspond to the physical characteristics, said error calculating means (3) are provided for receiving the cumulatively added result output from the cumulative addition means (2) and a measured value of physical quantity which is obtained by the wave receiving means (i), and for calculating an error, and said correcting means (11a, 12a, ..., 1ma) are provided for correcting values of each physical formula operating means (11, 12, ..., 1m) based upon the calculated error, the value regulating a cause of the physical quantity. An apparatus as set forth in claim 2 for reducing line spectrum noise for reducing line spectrum noises from a measured signal which includes line spectrum noises, wherein said physical formula operating means (11, 12, ..., 1m) are provided for operating the physical formulae corresponding to line spectrum noises, a number of physical formula operating means corresponds to a number of species of line spectrum noises, and said correcting means (11a, 12a, ..., 1ma) are provided for correcting a value of each physical formula operating means based upon calculated error, the value regulating a cause of the line spectrum, and wherein said result collecting means is omitted.

Inventors: UEDA TOMOAKI, UEDA, TOMOAKI
Applicants: DAIKIN IND LTD, DAIKIN INDUSTRIES, LTD.
Application year: 1992
Main classification: G01D3
IPC codes: G01N37, G01N33, G01D3