Update pca_audio_features.ipynb

2D subplots for 3 features
view angle

Author: fs446

pca_audio_features.ipynb

@@ -144,6 +144,60 @@
     "## Plot Directions of the PCs and Original Data"
    ]
   },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "3d0d6363",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "fig = plt.figure()\n",
+    "fig.set_size_inches(9,3)\n",
+    "\n",
+    "ax = fig.add_subplot(131)\n",
+    "ax.plot(X[:, 0], X[:, 1], \".\", color=\"gray\", ms=1)\n",
+    "ax.axis(\"square\")\n",
+    "ax.set_xlabel(\"original feature 1: true_peak_lin\")\n",
+    "ax.set_ylabel(\"original feature 2: rms_db\")\n",
+    "ax.grid(True)\n",
+    "\n",
+    "ax = fig.add_subplot(132)\n",
+    "ax.plot(X[:, 0], X[:, 2], \".\", color=\"gray\", ms=1)\n",
+    "ax.axis(\"square\")\n",
+    "ax.set_xlabel(\"original feature 1: true_peak_lin\")\n",
+    "ax.set_ylabel(\"original feature 3: lufs_db\")\n",
+    "ax.grid(True)\n",
+    "\n",
+    "ax = fig.add_subplot(133)\n",
+    "ax.plot(X[:, 1], X[:, 2], \".\", color=\"gray\", ms=1)\n",
+    "ax.axis(\"square\")\n",
+    "ax.set_xlabel(\"original feature 2: rms_db\")\n",
+    "ax.set_ylabel(\"original feature 3: lufs_db\")\n",
+    "ax.grid(True)\n",
+    "\n",
+    "plt.tight_layout()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "c5dd0957",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "fig = plt.figure()\n",
+    "ax = fig.add_subplot(111, projection=\"3d\")\n",
+    "ax.plot(X[:, 0], X[:, 1], X[:, 2], \".\", color=\"gray\", ms=1)\n",
+    "ax.axis(\"square\")\n",
+    "ax.set_xlabel(\"original feature 1: true_peak_lin\")\n",
+    "ax.set_ylabel(\"original feature 2: rms_db\")\n",
+    "ax.set_zlabel(\"original feature 3: lufs_db\")\n",
+    "ax.set_title(\"data cloud in original coordinate system\")\n",
+    "ax.grid(True)\n",
+    "ax.azim = -44\n",
+    "ax.elev = 28"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -174,7 +228,9 @@
     "ax.set_zlabel(\"original feature 3: lufs_db\")\n",
     "ax.set_title(\"data cloud in original coordinate system\")\n",
     "ax.legend()\n",
-    "ax.grid(True)"
+    "ax.grid(True)\n",
+    "ax.azim = -44\n",
+    "ax.elev = 28"
    ]
   },
   {
@@ -235,7 +291,9 @@
     "ax.set_zlabel(\"PC 3\")\n",
     "ax.set_title(\"data cloud in PC coordinate system\")\n",
     "ax.legend()\n",
-    "ax.grid(True)"
+    "ax.grid(True)\n",
+    "ax.azim = -37\n",
+    "ax.elev = 28"
    ]
   },
   {
@@ -306,7 +364,9 @@
     "ax.set_zlabel(\"new feature 3\")\n",
     "ax.set_title(\"rank-{0:d} approximation of data\".format(r_des))\n",
     "ax.legend()\n",
-    "ax.grid(True)"
+    "ax.grid(True)\n",
+    "ax.azim = -44\n",
+    "ax.elev = 28"
    ]
   },
   {
@@ -380,6 +440,8 @@
     "plt.title(\"PCA data dimensionality reduction to {0:d}D\".format(dim_des))\n",
     "ax.legend()\n",
     "ax.grid(True)\n",
+    "ax.azim = -37\n",
+    "ax.elev = 28\n",
     "\n",
     "print(\"reduced to dimension\", X_dim_red.shape)"
    ]