qq_41985559的博客

原创 【代码】《C程序设计》谭浩强.清华大学出版社

例4.1在例3.5的基础上对程序进行改造。题目要求解得ax^ 2+bx+c=0方程的根。由键盘输入a,b,c，假设a,b,c的值任意，并不保证b ^ 2-4ac≥0。要求在程序中进行判别，如果b^ 2-4ac≥0，就计算并输出方程的两个实根，如果b^2-4ac<0，就输出“方程无实根”的信息。#include<stdio.h>#include<math.h>in...

原创 【代码】《C程序设计》谭浩强.清华大学出版社

例1.1要求在屏幕上输出一以下一行信息。This is a C program.#include<stdio.h>int main(){ printf("This is a C program.\n"); return 0;}例1.2求两个整数之和#include<stdio.h>int main(){ int a, b, su...

原创 【Matlab代码】分形树木

程序设计的思路：1，我们先设定主干；2，我们将主干设定的比例缩小；3，将第二步得到的缩小的主干旋转一定的角度，然后接到主干上面预定的分叉的位置。

原创 【Matlab代码】Sierpinsk地毯

1，将一个大正方形作为最基本的图形；2，将第一步得到的数组缩小到三分之一；3，将第二部得到的数组复制八份，平移到相应的位置。

原创 【Matlab代码】Koch曲线

步骤：1，将基本的图形即线段用数组表示出来。由于是线段，我们只需要列出端点坐标即可；2，将第一步得到的数组做如下处理：首先将线段缩小到原来的三分之一。然后将得到的线段复制四份：将第一段放在最左边，将第二段逆时针旋转60°，左端点与第一段的右端点相接；将第三段顺时针旋转60°，左端点与第二段的右端点相接；最后将第四段左端点与第三段右端点相接。这样我们得到了最初屏幕上的第五个图。...

原创 【Matlab代码】斑点猫头鹰的生存问题的演化过程

A=[0 0 0.33;0.18 0 0;0 0.71 0.94];[pc,lambda]=eig(A);%求A的特征值和对应的特征向量[Y,I]=sort(diag(abs(lambda)),'descend')%对特征值的绝对值降序排列temp=diag(lambda);lambda=temp(I)%输出按特征值的绝对值降序排列的特征值Lambda_norm=[norm(lambda(1);norm(lambda(2));norm(lambda(3))]%对三个特征值的模pc=pc(:,I)

原创 【Python代码】数据排序

1…sort_index()方法在指定轴上根据索引进行排序，默认升序.sort_index(axis=0,ascending=True)import pandas as pdimport numpy as npb = pd.DataFrame(np.arange(20).reshape(4,5), index=['c','a', 'd', 'b'])b2…sort_index(axis=0,ascending=True)import pandas as pdimport numpy as

原创 【Python代码】散点图的绘制

import numpy as npimport matplotlib.pyplot as pltfig, ax = plt.subplots()ax.plot(10*np.random.randn(100),10*np.random.randn(100),'o')ax.set_title('Simple Scatter')plt.show()扫码关注“图像处理与模式识别研究所”解锁更多技能哟。...

原创 【Python代码】极坐标图的绘制

import numpy as npimport matplotlib.pyplot as pltN = 20theta = np.linspace(0.0, 2*np.pi, N, endpoint = False)radii = 10*np.random.rand(N)width = np.pi / 4*np.random.rand(N)ax = plt.subplot(111, projection = 'polar')bars = ax.bar(theta, radii, width=

原创 【Python代码】直方图的绘制

import numpy as npimport matplotlib.pyplot as pltnp.random.seed(0)mu, sigma = 100, 20#均值和标准差a = np.random.normal(mu, sigma, size = 100)plt.hist(a, 20, normed = 1, histtype = 'stepfilled', facecolor = 'b', alpha = 0.75)plt.title('Histogram')plt.show(

原创 【Python代码】饼图的绘制

import matplotlib.pyplot as pltlabels = 'Frogs', 'Hogs', 'Dogs', 'Logs'sizes = [15, 30, 45, 10]explode = (0, 0.1, 0, 0)plt.pie(sizes, explode=explode, labels=labels, autopct = '%1.1f%%', shadow= False, startangle =90)plt.show()import matplo

原创 【Python代码】图像的手绘效果

from PIL import Imageimport numpy as npa = np.asarray(Image.open('C:/Users/xpp/Desktop/Lena.png').convert('L'))depth = 10.#(0-100)grad = np.gradient(a)#取图像灰度的梯度值grad_x,grad_y=grad#分别取横纵图像梯度值grad_x=grad_x*depth/100.grad_y=grad_y*depth/100.A = np.sqr

原创 【Python代码】图像的变换

from PIL import Imageimport numpy as npa = np.array(Image.open("C:/Users/xpp/Desktop/Lena.png"))print(a.shape, a.dtype)b = [255, 255, 255]-aim = Image.fromarray(b.astype('uint8'))im.save('C:/Users/xpp/Desktop/Lena2.png')注意：b = [255, 255, 255]-a

原创 【Python代码】图像的数组表示

from PIL import Imageimport numpy as npim = np.array(Image.open("C:/Users/xpp/Desktop/Lena.png"))print(im.shape,im.dtype)

原创 【Python代码】计算A^2+B^3，其中，A和B是一维数组

import numpy as npdef npSum(): a = np.array([0, 1, 2, 3, 4]) b = np.array([9, 8, 7, 6, 5]) c = a**2 + b**3 return cprint(npSum())

原创 【Matlab代码】不同参数的线性变换效果

原创 【Matlab代码】多幅图像D，存储了一组索引图像，MAP为颜色索引表。分别以上述的3种方式显示它们

load mri%载入Matlab自带的核磁共振图像imshow(D(:,:,7),map);%显示多幅中的一幅figure,montage(D,map);%同一窗口显示%转化为电影figuremove=(immovie(D,map));%设定颜色表movie(move);%播放电影

原创 【Matlab代码】显示多幅图像

%ex2_7.mI=imread('pout.tif');%读取图像%在不同窗口显示figure;%创建一个新的窗口imshow(I);figure;imshow(I,[]);%在相同窗口显示figure;subplot(1,2,1);imshow(I);subplot(1,2,2);imshow(I,[]);

原创 【Matlab代码】图像文件的读取、显示、回写

%ex2_6.m%读取图像文件I=imread('gantrycrane.png');%显示图像imshow(I);%写回到文件imwrite(I,'gantrycrane.tif','TIFF');

原创 【Matlab代码】查看图像文件信息

imfinfo('pout.tif');%查看图像文件信息

原创 【Matlab代码】读入一幅tif图像文件，并在写入磁盘时将tif图像转换为bmp图像

I=imread('pout.tif');%读入图像whos I%查看图像变量信息%通过whos命令可以看到读入的高为291、宽为240的灰度图像I就是一个291*240的二维矩阵imwrite(I,'pout.bmp');%将图像写如文件pout.bmp,同时起到了转换文件类型的作用...

原创 【Matlab代码】产生一幅亮度按对角线方向的余弦规律变化的灰度图，比较一维方法和二维方法所需的时间

A=rand(3000,3000);f=zeros(3000,3000);u0=100;v0=100;tic;%开始计时%一维方法%外层for循环开始for r=3000 u0x=u0*(r-1) %内层for循环开始 for c=1:3000 v0y=v0*(c-1); f(r,c)=A(r,c)*cos(u0x+v0y); %内层for循环结束 end%外层for循环结束endt1=toc%停止计时并记录时

原创 【Matlab代码】使用switch分支和while循环

%ex2_2.marg=input('Input argument:');total=0;detail=0;%switch语句开始switch arg %分支1 case 1 i=1; %外层while语句开始 while(i<=5) total=total+1; i=i+1; j=1; %内层while语句开始

原创 【Matlab代码】if语句和for循环及其嵌套

%ex2_1.marg=input('input argument:');%提示输入arg变量total=0;detail=0;%if语句开始if(arg==1) %外层for语句开始 for i=1:1:5 total=total+1; %内层for语句开始 for j=1:0.1:2 detail=detail+total; %内层for语句结束 end %外层

原创 【Matlab代码】捕食者与被捕食者问题的演化过程

求A的特征值与特征向量。A=[0.5 0.4;-0.104 1.1];[pc,lambda]=eig(A);%求A的特征值和对应的特征向量[Y,I]=sort(diag(abs(lambda)),'descend');%对特征值的绝对值降序排列temp=diag(lambda);lambda=temp(I)%输出按特征值的绝对值降序排列pc=pc(:,I)%与特征值对应的特征向量...

原创 【Matlab代码】特征值与特征向量基本概念

函数：[V,D]=eig(A)或[V,D]=eig(A,‘nobalance’)功能：求矩阵A的特征值和特征向量说明：（1）生成特征矩阵D和特征向量构成的矩阵V，使得AV=VD。矩阵D由A的特征值在主对角线构成的对角矩阵。V是由A的特征向量按列构成的矩阵。（2）[V,D]=eig(A）中，先对A作相似变换再求A的特征值和特征向量；而[V,D]=eig(A,‘nobalance’)中，直接求矩阵A的特征值和特征向量。...

原创 【Matlab代码】有持续外力下单自由度阻尼系统的数值解

function dx=odefun4(t,x)global c;%全局变量dx=zeros(2,1);%存储x1,x2的导数dx(1)=x(2);%第一个方程dx(2)=20*sin(5*t)-20*c*x(2)-100*x(1);%第二个方程注意：修改为：dx(2)=20sin(5t)-20cx(2)-100*****x(1);%第二个方程

原创 【Matlab代码】外力单自由度阻尼系统的数值解

function dx=odefun2(t,x)dx=zeros(2,1);%存储x1,x2的导数dx(1)=x(2);%第一个方程dx(2)=-20*x(2)-100*x(1);%第二个方程[t,x]=ode45('odefun2',[0,4],[1;0])

原创 【Matlab代码】无外力下单自由度阻尼系统的解析解

x=dsolve('D2x+2*c*w*Dx+w^2*x=0')

原创 【Matlab代码】微分方程的解析解和数值解

Matlab求微分方程的解析解用法：dsolve(‘equation’,‘condition;,‘v’)说明：（1）equation是方程式，condition是条件，v是自变量（确省为t）（2）若不带条件，则解中带积分常数。（3）如果没有显示解，则系统尝试给出隐氏解。（4）如果无隐氏解，则返回空符号。格式：（1）y’表示为Dy，y’‘表示为D2y，依次类推。（2）有多个方程或多个条件时，写多个对应的参数即可。求微分方程y’’-5y’+6y=e^axy=dsolve('D2y-5*Dy+6

原创 【Matlab代码】生成物浓度散点图

化学反应生成物浓度

原创 【Matlab代码】斐波拉契数列

实验：显示数列的函数类型function fib1(n)%显示数列前n项fn=[1,1];%将前两项放到数组fn中for i=3:n fn=[fn,fn(i-2)+fn(i-1)];%将第i项添加到数组fn中endplot(fn)%前n项显示成曲线

原创 【Matllab代码】数据拟合简介

做数据拟合的方法：在matlab软件中，做数据拟合的函数是polyfit。函数：polyfit(x,y,n)功能：用n阶多项式拟合数据(x,y)，使得在数据点处误差的平方和最小。说明：参数x和y都是数组，分别存放数据的横坐标和竖坐标；参数n是多项式的阶，是人为指定的。函数polyfit的返回值，只是多项式的系数。如果需要计算多项式的值，调用函数polyva即可。函数：polyval(p,x)功能：根据多项式的系数，计算多项式的值。说明：参数p，是函数polyfit的返回值，参数x，是指定的横

原创 【Matlab代码】程序设计

(1)for循环：for i=范围表达式可执行语句end(2)while循环：while 逻辑表达式可执行语句end(3)if-else-end分支：if 逻辑表达式可执行语句else if 逻辑表达式可执行语句else可执行语句end说明：上述结构中，可以没有else if段和else段。（4）switch-case-end分支：switch 表达式case 常量表达式1可执行语句块1case 常量表达式2可执行语句块2case 常量表达式n可执行语句块n

原创 【Matlab代码】M文件

将要执行的命令，写在1个文件中，就是M文件。如果需要执行的命令较多，或者需要重复执行命令，就应该采用M文件方式。MATLAB自带M文件编辑器Notebook，M文件的扩展名是.m。M文件应保存在current directory（当前工作目录）下。...

原创 【Matlab代码】二维特殊图形

常见的统计图，都有对应的函数。

原创 【Matlab代码】hold指令

将多次绘制的曲线，显示在同一图形窗口中。函数：hold on，hold off功能：hold on，开启图形保持功能；hold off，关闭图形保持功能。说明：MATLAB默认的配置是hold off。

原创 【Matlab代码】用ezplot函数绘制

函数：ezplot(F,[xmin,xmax])功能：绘制符号函数F在区间[xmin,xmax]内的图像。参数说明：参数F是函数表达式，只含有一个变量，参数区间可缺省，默认为[-2pi,2pi]。注释：可以用ezplot函数，绘制参数方程与隐函数的图形，不能设置图形的格式。ezplot('cos(4*x)*sin(10*x)*exp(-abs(x))',[-3,3,-1,1])说明：参数[-3,3,-1,1]是指定x坐标和y坐标的范围，这是为了得到更好的效果。...

原创 【Matlab代码】plot函数

绘制二维图形，最基本的函数是plot。函数：plot(X,Y,‘s’)功能：连接参数所给点列，形成一条或多条折线。参数说明：参数(X,Y)是点列的坐标，二者是同长度的数组或矩阵。参数s是图形格式字符串，可以设定曲线的线型和颜色。注释：横坐标紧密时，效果就是曲线图。格式字符串可以缺省，由系统自动设置格式。可以有多组参数。增加一组参数，就多画一条曲线。X=-4:0.1:4;Y=(exp(X)+exp(-X))*0.5;plot(X,Y)X=-4:0.1:4;Y1=(exp(X)+

原创 【Matlab代码】数组及其运算

数组就是，只有一行或一列的矩阵，数组具有矩阵的所有特性，同时还有创建数组的特殊指令，还有异于矩阵运算的数组运算。特殊数组的创建：linespace(a,b,n)给出区间[a,b]内的n个等分点构成的数组。数组的运算：...