Лаб6-checkpoint.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      " 0.04022  0.14100  0.13815  0.07365  0.06565  0.06870  0.18999  0.14881  0.04421  0.08962 \n",
      " 0.00348  0.11226  0.13446  0.03090  0.14224  0.08289  0.07262  0.16430  0.12442  0.13243 \n",
      " 0.04838  0.16432  0.04418  0.12634  0.04943  0.03904  0.16842  0.05663  0.17604  0.12721 \n",
      " 0.13889  0.01101  0.00429  0.11578  0.05741  0.18974  0.17680  0.17248  0.01107  0.12253 \n",
      " 0.19905  0.02702  0.20071  0.03148  0.15985  0.11023  0.03834  0.12005  0.04831  0.06494 \n",
      " 0.00084  0.01783  0.10516  0.02416  0.18764  0.06411  0.06456  0.17095  0.15092  0.21383 \n",
      " 0.06504  0.10690  0.03653  0.25945  0.05396  0.25311  0.02134  0.07118  0.01996  0.11252 \n",
      " 0.14592  0.13221  0.15450  0.05521  0.17417  0.16056  0.00020  0.07417  0.09380  0.00927 \n",
      " 0.09346  0.06010  0.04621  0.15095  0.14789  0.08156  0.00791  0.16457  0.05599  0.19134 \n",
      " 0.03520  0.14376  0.03575  0.17344  0.00310  0.17769  0.07706  0.03686  0.15771  0.15943 \n",
      "\n",
      "\n",
      "   14.00    13.00     1.00    19.00     5.00     5.00    18.00    10.00     8.00     7.00 \n",
      "\n",
      "(9.58, 9.58, 9.58, 9.58, 9.58, 9.58, 9.58, 9.58, 9.58, 9.58)\n",
      "\n",
      "   30.00   -60.00    30.00    19.00     5.00   -60.00    30.00    10.00     8.00     7.00 \n",
      "\n",
      "(-0.0977, -0.0977, -0.0977, -0.0977, -0.0977, -0.0977, -0.0977, -0.0977, -0.0977, -0.0977)\n"
     ]
    }
   ],
   "source": [
    "import random\n",
    "\n",
    "def matrix_gen(n):\n",
    "    A = matrix(RR, n)\n",
    "    for i in range(n) :\n",
    "        s = 0\n",
    "        for j in range(n) :\n",
    "            x = random.uniform(0, 1)\n",
    "            A[i, j] = x\n",
    "            s += x\n",
    "        for j in range(n) :\n",
    "            A[i, j] *= s^(-1)\n",
    "    return A\n",
    "\n",
    "def opinion_gen(n) :\n",
    "    x =[]\n",
    "    for i in range(n) :\n",
    "        x.append(random.randint(1, 20))\n",
    "    y = vector(RR, n)\n",
    "    for i in range(n) :\n",
    "        y[i] = x[i]\n",
    "    return y\n",
    "\n",
    "def matrix_norm(A, n) :\n",
    "    norm = 0\n",
    "    for i in range(10) :\n",
    "        s = 0\n",
    "        for j in range(10) :\n",
    "            s += A[j, i]^2\n",
    "        norm += sqrt(s)\n",
    "    return norm\n",
    "\n",
    "def matrix_res(A, n) :\n",
    "    eps = 1\n",
    "    i = 2\n",
    "    C = A\n",
    "    B = A\n",
    "    while eps > 10^(-6) :\n",
    "        C = B\n",
    "        B = A^i\n",
    "        eps = matrix_norm(B - C, n)\n",
    "        i += 1\n",
    "    #print(\" N =\",i)\n",
    "    return B\n",
    "\n",
    "A = matrix_gen(10)\n",
    "for i in range(10):\n",
    "    for j in range(10): \n",
    "        print('{0:8.5f}'.format(A[i,j]), end = ' ') \n",
    "    print()\n",
    "print()\n",
    "print()\n",
    "\n",
    "x = opinion_gen(10)\n",
    "for i in x:\n",
    "    print('{0:8.2f}'.format(i), end = ' ') \n",
    "print()\n",
    "    \n",
    "A.numerical_approx(digits = 3)\n",
    "x.numerical_approx(digits = 3)\n",
    "\n",
    "B = matrix_res(A, 10)\n",
    "B.numerical_approx(digits = 3)\n",
    "\n",
    "x_res = B*x\n",
    "print()\n",
    "print(x_res.numerical_approx(digits = 3))\n",
    "\n",
    "x[1] = -60\n",
    "x[5] = -60\n",
    "x[2] = 30\n",
    "x[0] = 30\n",
    "x[6] = 30\n",
    "\n",
    "\n",
    "print()\n",
    "for i in x:\n",
    "    print('{0:8.2f}'.format(i), end = ' ') \n",
    "print()\n",
    "\n",
    "x_res = B*x\n",
    "x_res.numerical_approx(digits = 3)\n",
    "print()\n",
    "print(x_res.numerical_approx(digits = 3))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "SageMath 9.2",
   "language": "sage",
   "name": "sagemath"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}