Основы компьютерного моделирования физических и технических систем. Математическую модель указанной физической системы, страница 6

double[,] ci = {{2.0,1.0},   //матрица демпфирования

{1.0,2.0}};

double ki_m1_kaf = (A * lam1) / lel;

double ki_m2_kaf = (A * lam2) / lel;

double[,] ki = {{ 1.0,-1.0},

{-1.0, 1.0}};

double[,] tmp = new double[ci.GetLength(0), ci.GetLength(1)];

for (int i = 0; i < C.GetLength(0)-1; i++)

{

Array.Copy(ci,tmp,ci.Length);

matrNaConst(tmp, ci_m1_kaf);

priSummMatr(C,tmp , i);

}

for (int i = 0; i < K.GetLength(0) - 1; i++)

{

Array.Copy(ki, tmp, ki.Length);

matrNaConst(tmp, ki_m1_kaf);

priSummMatr(K, tmp, i);

}

for (int i = 0; i < tl; i++)

{

B = summMatr(K, 1, C, 2 / dT);

P = summMatr(K, -1, C, 2 / dT);

//формирование вектора F

F = new double[kolEl+1];

//for (int i = 0; i < kolEl-1; i++)

Q = 0;

q = 0;

h = double.Parse(TB_hl.Text);

Ts = double.Parse(TB_Tsl.Text);

double[] tm = fi(Q, q, h, Ts,T[i,0], lel, A);

F[0] = tm[0];

F[1] = tm[1];

Q = 0;

q = double.Parse(TB_qr.Text);

h = 0;

Ts = 0;

tm = fi(Q, q, h, Ts, T[i, 0], lel, A);

F[kolEl - 1]= tm[0];

F[kolEl] = tm[1];

//}

double[] ot;

double[] tmpT = new double[T.GetLength(1)];

//double[,] Btmp = new double[B.GetLength(0), B.GetLength(1)];

for (int p = 0; p < tmpT.GetLength(0); p++)

tmpT[p] = T[i, p];

ot = multiMatr(P, 1, tmpT, 1);

double[] Rpart = summVector(ot, 1, F, 2);

reshSLAU(B,Rpart,i+1);

}

string heder = "<?xml version=\"1.0\"?> \n<?mso-application progid=\"Excel.Sheet\"?>\n<Workbook xmlns=\"urn:schemas-microsoft-com:office:spreadsheet\"\n xmlns:o=\"urn:schemas-microsoft-com:office:office\"\n xmlns:x=\"urn:schemas-microsoft-com:office:excel\"\n xmlns:ss=\"urn:schemas-microsoft-com:office:spreadsheet\" \n xmlns:html=\"http://www.w3.org/TR/REC-html40\"> \n <DocumentProperties xmlns=\"urn:schemas-microsoft-com:office:office\">\n <Author>Igor</Author>\n  <LastAuthor>Igor</LastAuthor>  \n<Created>2012-02-17T22:38:42Z</Created>  \n<LastSaved>2012-02-17T22:50:27Z</LastSaved> \n <Version>14.00</Version> \n</DocumentProperties> \n<OfficeDocumentSettings xmlns=\"urn:schemas-microsoft-com:office:office\">  \n<AllowPNG/> \n</OfficeDocumentSettings> \n<ExcelWorkbook xmlns=\"urn:schemas-microsoft-com:office:excel\">  \n<WindowHeight>7995</WindowHeight>  \n<WindowWidth>20115</WindowWidth>  \n<WindowTopX>240</WindowTopX>  \n<WindowTopY>75</WindowTopY>  \n<ProtectStructure>False</ProtectStructure>  \n<ProtectWindows>False</ProtectWindows> \n</ExcelWorkbook> \n<Styles>  \n<Style ss:ID=\"Default\" ss:Name=\"Normal\">   \n<Alignment ss:Vertical=\"Bottom\"/>   \n<Borders/>   \n<Font ss:FontName=\"Calibri\" x:CharSet=\"204\" x:Family=\"Swiss\" ss:Size=\"11\"    \nss:Color=\"#000000\"/>   \n<Interior/>   \n<NumberFormat/>   \n<Protection/>  \n</Style> \n</Styles> \n<Worksheet ss:Name=\"Лист1\">  \n";

string footor = "</Table> \n <WorksheetOptions xmlns=\"urn:schemas-microsoft-com:office:excel\"> \n  <PageSetup>  \n  <Header x:Margin=\"0.3\"/>   \n <Footer x:Margin=\"0.3\"/>    \n<PageMargins x:Bottom=\"0.75\" x:Left=\"0.7\" x:Right=\"0.7\" x:Top=\"0.75\"/>   \n</PageSetup>   \n<Unsynced/>   \n<Selected/>   \n<Panes>    \n<Pane>     \n<Number>3</Number>     \n<ActiveRow>2</ActiveRow>     \n<ActiveCol>3</ActiveCol>    \n</Pane>   \n</Panes>   \n<ProtectObjects>False</ProtectObjects>   \n<ProtectScenarios>False</ProtectScenarios>  \n</WorksheetOptions> \n</Worksheet> \n</Workbook> ";