Huber稳健估计

鲁棒估计中的Huber估计，C#源码，仅供参考

using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Windows.Forms;
using System.IO;

namespace Huber
{
    public partial class Form1 : Form
    {
        public Form1（）
        {
            InitializeComponent（）;
        }

        //输入文件，作为公共对象
        public StreamReader f1 = new StreamReader（“fin.txt“）;

        private void button1_Click（object sender EventArgs e）
        {
            int M = int.Parse（textBox1.Text）;
            //资料个数从输入框中读取
            int N = int.Parse（textBox2.Text）;
            //未知数个数

            //-----------------读数据文件
            MAT A = new MAT（M N）;
            //声明矩阵对象
            VEC b = new VEC（M）;
            VEC P = new VEC（M）;

            VEC x = new VEC（N）;
            //自变量

            //从文件中读入数据到系数矩阵
            for （int i = 0; i < M; i++）
            {
                //整行读进来
                string line = f1.ReadLine（）;
                //分割
                string[] strs = line.Split（new char[] { ‘\t‘ ‘ ‘ }
                    StringSplitOptions.RemoveEmptyEntries）;

                for （int j = 0; j < strs.Length; j++）
                    A[i j] = Double.Parse（strs[j]）;

            }

            f1.ReadLine（）;

            for （int i = 0; i < M; i++）
                b[i] = Double.Parse（f1.ReadLine（））;


            for （int i = 0; i < M; i++）
                P[i] = 1.0;
            //对所有数据赋予等权


            //------------------------------


            cHuber ob = new cHuber（）;

            //粗差系统最小二乘计算
            ob.wls（A b P out x）;
            richTextBox1.AppendText（“含粗差系统最小二乘计算结果 \n“）;
            for （int i = 0; i < N; i++）
                richTextBox1.AppendText（x[i].ToString（） + “\n“）;

            //粗差系统IGGI
            ob.Huber（A b P out x）;
            richTextBox1.AppendText（“\n“）;
            richTextBox1.AppendText（“含粗差系统Hampel计算结果 \n“）;
            for （int i = 0; i < N; i++）
                richTextBox1.AppendText（x[i].ToString（） + “\n“）;

            f1.ReadLine（）;
            for （int i = 0; i < M; i++）
                b[i] = Double.Parse（f1.ReadLine（））;

            ob.wls（A b P out x）;
            richTextBox1.AppendText（“\n“）;
            richTextBox1.AppendText（“无粗差系统最小二乘计算结果 \n“）;
            for （int i = 0; i < N; i++）
                richTextBox1.AppendText（x[i].ToString（） + “\n“）;

            f1.Close（）;
            //关闭文件

        }
    }
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       2130  2011-09-09 16:53  3 Huber稳健估计\Huber\bin\Debug\fin.txt

     文件      17408  2011-09-09 16:58  3 Huber稳健估计\Huber\bin\Debug\Huber.exe

     文件      50688  2011-09-09 16:58  3 Huber稳健估计\Huber\bin\Debug\Huber.pdb

     文件      11600  2011-11-12 15:04  3 Huber稳健估计\Huber\bin\Debug\Huber.vshost.exe

     文件        490  2010-03-17 22:39  3 Huber稳健估计\Huber\bin\Debug\Huber.vshost.exe.manifest

     文件       2840  2011-09-09 16:58  3 Huber稳健估计\Huber\Form1.cs

     文件       5337  2011-09-09 16:56  3 Huber稳健估计\Huber\Form1.Designer.cs

     文件       5817  2011-09-09 16:56  3 Huber稳健估计\Huber\Form1.resx

     文件      10550  2011-09-09 16:29  3 Huber稳健估计\Huber\Huber.cs

     文件       3742  2011-09-09 16:33  3 Huber稳健估计\Huber\Huber.csproj

     文件      17074  2011-09-09 11:10  3 Huber稳健估计\Huber\matvec.cs

     文件       2152  2011-11-12 15:03  3 Huber稳健估计\Huber\obj\x86\Debug\DesignTimeResolveAssemblyReferences.cache

     文件       6211  2011-11-12 15:03  3 Huber稳健估计\Huber\obj\x86\Debug\DesignTimeResolveAssemblyReferencesInput.cache

     文件        236  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\GenerateResource.read.1.tlog

     文件        562  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\GenerateResource.write.1.tlog

     文件        658  2011-11-12 15:04  3 Huber稳健估计\Huber\obj\x86\Debug\Huber.csproj.FileListAbsolute.txt

     文件      17408  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\Huber.exe

     文件        180  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\Huber.Form1.resources

     文件      50688  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\Huber.pdb

     文件        180  2011-09-09 16:58  3 Huber稳健估计\Huber\obj\x86\Debug\Huber.Properties.Resources.resources

     文件        486  2011-09-09 16:13  3 Huber稳健估计\Huber\Program.cs

     文件       1360  2011-09-09 16:13  3 Huber稳健估计\Huber\Properties\AssemblyInfo.cs

     文件       2860  2011-09-09 16:13  3 Huber稳健估计\Huber\Properties\Resources.Designer.cs

     文件       5612  2011-09-09 16:13  3 Huber稳健估计\Huber\Properties\Resources.resx

     文件       1090  2011-09-09 16:13  3 Huber稳健估计\Huber\Properties\Settings.Designer.cs

     文件        249  2011-09-09 16:13  3 Huber稳健估计\Huber\Properties\Settings.settings

     文件        857  2011-09-09 16:13  3 Huber稳健估计\Huber.sln

    ..A..H.     20992  2011-11-12 16:49  3 Huber稳健估计\Huber.suo

     目录          0  2012-08-20 08:47  3 Huber稳健估计\Huber\obj\x86\Debug\TempPE

     目录          0  2013-04-27 11:33  3 Huber稳健估计\Huber\obj\x86\Debug

............此处省略11个文件信息