{ "cells": [ { "cell_type": "code", "execution_count": 307, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import numpy as np\n", "\n", "from keras.layers import Dense\n", "from keras.models import Sequential\n", "from keras import optimizers" ] }, { "cell_type": "code", "execution_count": 916, "metadata": { "collapsed": false }, "outputs": [], "source": [ "X = np.random.uniform(low=0, high=1, size=(300, 30))\n", "y = np.array([0, 1] * (X.shape[0] // 2))" ] }, { "cell_type": "code", "execution_count": 917, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# n_hidden = X.shape[0]\n", "# W0 = np.ones_like(X.T)\n", "# b0 = -X.ravel() + 0.001\n", "# weights_0 = [W0, b0]" ] }, { "cell_type": "code", "execution_count": 918, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# hidden_activations = np.maximum(np.dot(X, W0) + b0, 0)\n", "# np.linalg.matrix_rank(hidden_activations)" ] }, { "cell_type": "code", "execution_count": 919, "metadata": { "collapsed": true }, "outputs": [], "source": [ "n_hidden = X.shape[0]\n", "W0 = X / (np.linalg.norm(X, axis=1, keepdims=True) ** 2 + 1e-8)\n", "W0 = W0.T\n", "b0 = -0.98\n", "weights_0 = [W0, b0]" ] }, { "cell_type": "code", "execution_count": 920, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "300" ] }, "execution_count": 920, "metadata": {}, "output_type": "execute_result" } ], "source": [ "hidden_activations = np.maximum(np.dot(X, W0) + b0, 0)\n", "np.linalg.matrix_rank(hidden_activations)" ] }, { "cell_type": "code", "execution_count": 921, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "0.0" ] }, "execution_count": 921, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.linalg.det(hidden_activations)" ] }, { "cell_type": "code", "execution_count": 931, "metadata": { "collapsed": true }, "outputs": [], "source": [ "n_hidden = X.shape[0]\n", "W0 = X\n", "W0 = W0.T\n", "b0 = -0.98 * np.linalg.norm(X, axis=1) ** 2\n", "weights_0 = [W0, b0]" ] }, { "cell_type": "code", "execution_count": 932, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "300" ] }, "execution_count": 932, "metadata": {}, "output_type": "execute_result" } ], "source": [ "hidden_activations = np.maximum(np.dot(X, W0) + b0, 0)\n", "np.linalg.matrix_rank(hidden_activations)" ] }, { "cell_type": "code", "execution_count": 933, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "5.278958284816851e-214" ] }, "execution_count": 933, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.linalg.det(hidden_activations)" ] }, { "cell_type": "code", "execution_count": 934, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([[ 0.20354264, 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0.05563168, 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0.01919508, 0.06842773, 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0.2154047 , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0.09565226, 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ],\n", " [ 0. , 0.18351664, 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0.28104634, 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ],\n", " [ 0. , 0. , 0.16006979, 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ,\n", " 0. , 0. , 0. , 0. , 0. ]])" ] }, "execution_count": 934, "metadata": {}, "output_type": "execute_result" } ], "source": [ "hidden_activations[:3]" ] }, { "cell_type": "code", "execution_count": 935, "metadata": { "collapsed": false }, "outputs": [], "source": [ "# W1 = np.linalg.solve(hidden_activations, y)[:, np.newaxis]\n", "W1 = np.dot(np.linalg.pinv(hidden_activations), y)[:, np.newaxis]\n", "b1 = np.zeros(1)\n", "weights_1 = [W1, b1]" ] }, { "cell_type": "code", "execution_count": 936, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([[ 8.12960810e-14],\n", " [ 1.00000000e+00],\n", " [ 1.53071833e-13],\n", " [ 1.00000000e+00],\n", " [ 5.55215829e-15]])" ] }, "execution_count": 936, "metadata": {}, "output_type": "execute_result" } ], "source": [ "np.dot(hidden_activations, W1)[:5]" ] }, { "cell_type": "code", "execution_count": 937, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Epoch 1/5\n", "1s - loss: 2.0071e-08 - acc: 1.0000\n", "Epoch 2/5\n", "0s - loss: 2.0071e-08 - acc: 1.0000\n", "Epoch 3/5\n", "0s - loss: 2.0071e-08 - acc: 1.0000\n", "Epoch 4/5\n", "0s - loss: 2.0071e-08 - acc: 1.0000\n", "Epoch 5/5\n", "0s - loss: 2.0071e-08 - acc: 1.0000\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 937, "metadata": {}, "output_type": "execute_result" } ], "source": [ "model = Sequential()\n", "first_layer = Dense(output_dim=n_hidden, input_dim=X.shape[1], activation='relu',\n", " weights=weights_0)\n", "first_layer.trainable = False\n", "model.add(first_layer)\n", "second_layer = Dense(output_dim=1, activation='linear', weights=weights_1)\n", "second_layer.trainable = False\n", "model.add(second_layer)\n", "\n", "optimizer = optimizers.Adam(lr=0.001)\n", "model.compile(optimizer=optimizer, loss='mse',\n", " metrics=['accuracy'])\n", "\n", "model.fit(X, y, nb_epoch=5, verbose=2)" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "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.5.2" } }, "nbformat": 4, "nbformat_minor": 1 }