Noisy Data Synthesizer
Generates simulated experimental observation pairs (x_i, y_i) corrupted by controlled Gaussian white noise ε_i ~ N(0, σ²).
Formulate the linear least-squares regression problem for noisy discrete data, construct the Vandermonde design matrix, assemble and solve the normal equations (Xᵀ X) β = Xᵀ y, and quantitatively determine the optimal regression coefficient β. Explore interactive numerical methods, evaluate error attenuation and verify convergence orders.
Ready to execute Step = 1.
Adjust Step, then execute computation.
Mathematical operations are computed deterministically using 64-bit IEEE 754 precision. Residuals reflect ordinary least-squares regression across successive observations.
Ordinary least squares is sensitive to rogue outliers because squaring residuals gives disproportionate weight to extreme deviations. The simulation executes linear least-squares regression over standardized data series across settings k = 1 to 6. The linear regression fit of the response y = β · k confirms the optimal regression parameter β = 2.500 with R² = 1.000000.
Generates simulated experimental observation pairs (x_i, y_i) corrupted by controlled Gaussian white noise ε_i ~ N(0, σ²).
Constructs the M × p design matrix X with column vectors corresponding to basis functions 1, x, x².
Evaluates the symmetric positive semi-definite Gram matrix and right-hand moment vector Xᵀ y.
Factors Xᵀ X = L Lᵀ and executes forward/backward substitution to solve for optimal regression weights β.
Computes the residual vector r = y - Xβ, sum of squared errors SSE, and coefficient of determination R².
Monitors the condition number κ(Xᵀ X) = (κ(X))² to assess roundoff amplification as polynomial degree increases.
Formulate the linear least-squares regression problem for noisy discrete data, construct the Vandermonde design matrix, assemble and solve the normal equations (Xᵀ X) β = Xᵀ y, and quantitatively determine the optimal regression coefficient β. Standardized input settings: 1.0, 2.0, 3.0, 4.0, 5.0, 6.0 step.
Load M = 100 observation points generated by the underlying physical law y = 2.500 x + 1.200 + ε.
Assemble the M × 2 Vandermonde matrix with row entries [1, x_i] for each observation point.
Compute the matrix product Xᵀ X and right-hand projection vector Xᵀ y using double precision accumulators.
Solve the symmetric positive definite 2 × 2 normal system via Cholesky factorization or direct Cramer elimination.
Verify that the residual vector r = y - Xβ is strictly orthogonal to the column space of X (Xᵀ r = 0).
Log the fitted parameter response across observation settings k = 1 to 6 and perform linear regression to confirm the target β = 2.500.
Points are your recorded observations. The solid regression line validates Least-Squares Curve Fitting mathematical scaling.
| # | Observation Setting (k) (step) | Fitted Model Output Reading (reading) | Observation Setting Index k | Model Response y = β · k | Remove |
|---|
No readings yet. Execute a computation and record it.