Статистические модели макроэкономики

Российский государственный гидрометеорологический университет

Экономический и социально-гуманитарный факультет

Кафедра Экономики и менеджмента

Курсовая работа

по дисциплине «Информатика»

Статистические модели макроэкономики

Санкт-Петербург

2008 г.

1. Постановка задачи.

Даны: вектор непроизводственного потребления, матрицы межотраслевого баланса и .

1. Рассчитать матрицу В.

2. Найти вектор валового выпуска х, обеспечивающий данный вектор потребления.

3. Рассчитать произведение вектора на матрицу Ах.

4. Рассчитать вектор .

Все расчеты произвести с использованием программы, написанной на алгоритмическом языке ПАСКАЛЬ.

program KURSOVOJ_PROJECT;

uses crt;

type int=0.10;

vec=array [0.10] of real;

matr=array [0.10,0.10] of real;

var x, x₁, ax, ax₁, c: vec;

e, a, a₁, b, b₁, br, br₁: matr;

m, n: int;

name:char;

procedure readvec (var x: vec; name: char; n: int);

var i: int;

begin

writeln (' Ввести вектор ', name,' размера n=', n); for i:=0 to n-1 do

begin

write (' ', name,'[', i:2,']='); readln (x[i])

end;

end;readvec

procedure writevec (var x: vec; name: char; n: int);

var i: int;

begin

writeln (' Вывести вектор ', name, ' размера n=', n); for i:=0 to n-1 do

begin

write (' ', name,'[', i:2,']='); writeln (x[i]: 10:5)

end;

end;writevec

procedure readmatr (var y: matr; name: char; n, m: int);

var i, j: int;

begin

writeln (' Ввести матрицу ', name,

' размера: n=', n,', * m=', m); for i:=0 to n-1 do

for j:=0 to m-1 do

begin

write (' ', name, '[', i:2,' ,', j:2, ']= ');

readln (y[i, j])

end;

end;readmatr

procedure writematr (var y: matr; name: char; n, m: int);

var i, j: int;

begin

writeln (' Вывести матрицу ', name,

' размера: n=', n,', * m=', m);

for i:=0 to n-1 do

for j:=0 to m-1 do

begin

write (' ', name, '[', i:2,' ,', j:2, ']= ');

writeln (y[i, j]: 8:4)

end;

end;writematr

procedure Em (var E: matr; n: int);

var i, j: int;

begin

for i:=0 to n do

for j:=0 to n do

if i=j then E[i, j]: =1 else E[i, j]: =0;

end;end E

procedure Rmatr (var a, b, c:matr; n, m: int);

var i, j: int;

begin

for i:=0 to n do

for j:=0 to m do

c[i, j]: =a[i, j]-b[i, j];

end;Rmatr

procedure obrmatr (Var AIS, AP: matr; n: int);

var i, j, k, l: integer;

a:matr;

begin

a:=AIS;

n:=n-1;

for k:=0 to n do begin

for j:=0 to n do

if (j<>k) then AP[k, j]: =-a[k, j]/a[k, k];

for i:=0 to n do

if (i<>k) then AP [i, k]: =a[i, k]/a[k, k];

for i:=0 to n do

for j:=0 to n do

if (i<>k) and (j<>k)

then AP[i, j]: = a[i, j]-a[i, k]*a[k, j]/a[k, k];

AP[k, k]:=1/a[k, k]; a:=AP;

end;

end; end obrmatr

procedure matrvec (var b: matr; var c, x: vec; n: int);

var i, j: int;

begin

for j:=0 to n-1 do begin

x[j]: =0;

for i:=0 to n-1 do

x[j]: =x[j]+b[i, j]*c[i]

end;

end; end matrvec

{ ТЕЛО ПРОГРАММЫ }

begin

clrscr;

readvec (c,'C', 2);

readmatr (a,'A', 2,2);

readmatr (a₁,'Z', 2,2);

Em (e, 2); rmatr (e, a, br, 2,2); rmatr (e, a₁, br₁, 2,2); obrmatr (br, b,2); obrmatr (br₁, b₁, 2); matrvec (B, C, x, 2); matrvec (B₁, C, x₁, 2); matrvec (A, x, ax, 2); matrvec (A₁, x₁, ax₁, 2);

if (ax[0]<=c[0]) and (ax[1]<=c[1]) then

writeln (' Экономика матрица A — эффективна') else

writeln (' Экономика матрица A — не эффективна');

if (ax₁[0]<=c[0]) and (ax₁[1]<=c[1]) then

writeln (' Экономика матрица A₁ — эффективна') else

writeln (' Экономика матрица A₁ — не эффективна');

readln

write matr (B,'B', 2,2);

write matr (B1,'Z', 2,2);

write vec (x,'x'.2);

write vec (x1,'y', 2);

end.

{Окончание программы}

1.Процедура ввода вектора

2.Процедура вывода вектора

3.Процедура ввода матрицы

4.Процедура вывода матрицы

5.Процедура вывода единичной матрицы

6.Процедура умножения вектора на матрицу

7.Процедура умножения вектора на матрицу

8.Процедура образования матрицы

Расчеты при вводе значений

При добавлении:

B[0,0]=4.6154

B[0,1]=1,0462

B[1,0]=3,0769

B[1,1]=2,0308

Вывести матрицу Z размера: n=2,m=2

Z[0,0]=2,400

Z[0,1]=0,800

Z[1,0]=0,800

Z[1,1]=1,600

Вывести вектор Х размера n=2

Х[0]=13,84 615

Х[1]=7,13 846

Вывести вектор x размера n=2

Y[0]=4,8000

Y[1]=5,6000