由两点坐标如何画出直线  方法1:利用直线方程 斜率加截距 方法2:数据拟合 %由两点坐标得数据拟合直线与画线 x = [,]; y = [,]; k = ((-)/(-));% 由两点坐标得到直线斜率 line = k*x+0.5;% 直线方程 xy = :;% 定义画线的 x 长度 line1 = k*xy+0.5; figure(),plot(xy,line1); % % x1=[,]; x2=[,]; X=[x1(),x2()];% 两点坐标的x值 Y=[x1(),x2()];% 两点坐

最近接触到了遗传算法以及利用遗传算法求最优解,所以就把这些相关的内容整理记录一下. 一.遗传算法简介(摘自维基百科) 遗传算法(英语:genetic algorithm (GA))是计算数学中用于解决最佳化的搜索算法,是进化算法的一种.进化算法最初是借鉴了进化生物学中的一些现象而发展起来的,这些现象包括遗传.突变.自然选择以及杂交等. 算法 选择初始生命种群 循环 评价种群中的个体适应度 以比例原则(分数高的挑中概率也较高)选择产生下一个种群. 改变该种群(交叉和变异) 直到停止循环的条件满足

MATLAB线性方程组的迭代求解法 作者:凯鲁嘎吉 - 博客园http://www.cnblogs.com/kailugaji/ 一.实验目的 1． 借助矩阵按模最大特征值,判断解方程组的Jacobi迭代法所得迭代序列的敛散性. 2． 会在Jacobi迭代法所得迭代序列收敛时,用修改后的Gauss-Seidel迭代法. 3． 会逐次超松驰迭代法. 二.实验原理 三.实验程序 四.实验内容 用上面前二种方法求解4元线性方程组的近似解,所选方程组尽可能可以用多种方法求得收敛解. 注:要注意判断迭代法

MATLAB PCHIP函数一阶求导分析 摘要:本文首先根据三次立方插值的一般表达式,得出分段三次立方插值时,每个小区间上的各次项系数.分析发现,三次项.二次项.一次项系数都与小区间端点处的一阶导数值相关,故需要求出端点处的一阶导数值,Matlab Pchip.m文件给出了方法.根据pchip.m的代码以及参考文献,归纳出求一阶导数值的公式.文章第三节着重分析内点处导数值公式的由来,以及与此相关的其它三个公式. 关键字:Matlab Pchip.m:分段三次立方插值:离散数据一阶求导:单调保形插

画decision boundary(直线) %% ============= Part 3: Optimizing using fminunc =============% In this exercise, you will use a built-in function (fminunc) to find the% optimal parameters theta. % Set options for fminuncoptions = optimset('GradObj', 'on', '

问题是这样,如果我们知道两个向量v1和v2,计算从v1转到v2的旋转矩阵和四元数,由于旋转矩阵和四元数可以互转,所以我们先计算四元数. 我们可以认为v1绕着向量u旋转θ​角度到v2,u垂直于v1-v2平面. 四元数q可以表示为cos(θ/2)​+sin(θ/2)​u,即:q0​=cos(θ/2)​,q1​=sin(θ/2)​u.x,q2=sin(θ/2)​u.y,q3=sin(θ/2)​u.z 所以我们求出u和θ/2即可,u等于v1与v2的叉积,不要忘了单位化:θ/2用向量夹角公式就能求. ma

相关Matlab函数:hist, bar, cdfplot, ksdensity (1) hist函数 n = hist(Y, x)  假设x是一个向量,返回x的长度个以x为中心的,Y的分布情况. 比如:假设x是一个5元素的向量,返回Y在以x为中心的,x长度个范围内数据直方分布. [n,xout] = hist(...)  返回n和xout.包括有数目频率和间隔位置.能够使用bar(xout, n)来绘制直方图. (2) bar函数 绘制条形图.bar(X,Y) 将Y矩阵的每一行化成一组条形图.

转自 http://www.cnblogs.com/tiandsp/archive/2012/01/07/2316006.html  感谢Dsp tian clearclc;img = imread('124.jpg');     % Read image from graphics file [m n t] = size(img);            % 获取图型大小   max = 0;min = 256;avg = 0;for i =1:1:m    for j = 1:1:n    

sigmoid.m文件 function g = sigmoid(z)%SIGMOID Compute sigmoid functoon% J = SIGMOID(z) computes the sigmoid of z. g = zeros(size(z));  初始化g ,z可以是一个数,一个向量或者一个矩阵 % ====================== YOUR CODE HERE ======================% Instructions: Compute the si

clc; clear all; close all; fx = @(x) -(0.4./sqrt(1 + x.^2) - sqrt(1+x.^2) .* (1- 0.4./(1 + x.^2))+x); [x0, f] =fminbnd(fx,0,2); % f利用负号求最小值 x = 0 :0.1: 2; y = feval(fx,x); %% ========maxvalue==== figure plot(x,-y,'b-','linewidth',2) hold on plot(x0,-

clear all; close all; clc img=imread('lena.bmp'); figure; imshow(uint8(img)); [m n]=size(img); img_max=0; img_min=256; img_avg=0; for i=1:1:m for j=1:1:n if img(i, j)<img_min img_min=img(i,j); end if img(i, j)>img_max img_max=img(i,j); end img_avg=i

clear;clc;close all % 初始化 m = [-6,-2,0,2,4,6]; sigma = [1,1,0.5,0.25,0.6,2]; h = [1,2,3,2,2.13,3.14]; x=-10:0.01:10; sumy = zeros(1,length(x)); % 绘制各分支高斯信号 figure; for i = 1:length(m) y=(h(i)/((sqrt(2*pi))*sigma(i)))*exp(-((x-m(i)).^2)/(2*sigma(i).^2

"QR_H.m" function [Q,R] = QR_tao(A) %输入矩阵A %输出正交矩阵Q和上三角矩阵R [n,n]=size(A); E = eye(n); X = zeros(n,); R = zeros(n); P1 = E; :n- s = -sign(A(k,k))*norm(A(k:n,k)); R(k,k) = -s; w = [A(,)+s,A(:n,k)']'; else w = [zeros(,k-),A(k,k)+s,A(k+:n,k)']'; R(:

至于为什么,上图.转载于MZW_BG 枚举正方形的一条边,此时有上正方形和下正方形. 最后正方形个数/4,因为每个正方形被枚举了4条边 #include <bits/stdc++.h> using namespace std; int x[509],y[509],vis[200][200]; bool ok(int x,int y) { if(x<-50||x>50||y<-50||y>50) return false; if(vis[x+50][y+50]==0) r

matlab求极限(可用来验证度量函数或者隶属度函数)可用来验证是否收敛,取值范围等等. 一.问题来源 搜集聚类资料时,又看到了隶属度函数,没错,就是下面这个,期间作者提到m趋于2是,结果趋于1,我想验证下,于是查资料. 二.不同类型的极限 2.1 基础知识 a./b表示常数a除以矩阵b中每个元素或者矩阵a除以矩阵b对应元素或者常数b:点乘方a.^b,矩阵a中每个元素按b中对应元素乘方或者b是常数. 2.2 单变量独立式子 独立式子之地的是不存在连加之类的操作. 问题:用MATLAB求(x^2+

题目链接 题意 : 求一个多边形的核的面积. 思路 : 半平面交求多边形的核,然后在求面积即可. #include <stdio.h> #include <string.h> #include <iostream> #include <math.h> using namespace std ; struct node { double x; double y ; } p[],temp[],newp[];//p是最开始的多边形的每个点,temp是中间过程中临时

题目链接 题意 : 给你一个多边形,问你该多边形中是否存在一个点使得该点与该多边形任意一点的连线都在多边形之内. 思路 : 与3335一样,不过要注意方向变化一下. #include <stdio.h> #include <string.h> #include <iostream> #include <math.h> using namespace std ; struct node { double x; double y ; } p[],temp[],n

题目链接 题意 : 给你一个多边形,问你在多边形内部是否存在这样的点,使得这个点能够看到任何在多边形边界上的点. 思路 : 半平面交求多边形内核. 半平面交资料 关于求多边形内核的算法 什么是多边形的内核? 它是平面简单多边形的核是该多边形内部的一个点集,该点集中任意一点与多边形边界上一点的连线都处于这个多边形内部.就是一个在一个房子里面放一个摄像 头,能将所有的地方监视到的放摄像头的地点的集合即为多边形的核. 如上图,第一个图是有内核的,比如那个黑点,而第二个图就不存在内核了,无论点在哪里,总

本渣想回过头来整理一下MATLAB的一些基本的知识(很多东西比较琐碎,应该系统的梳理梳理),下文中没有提到的,自己用help查即可. 此文用来存个档,便于回顾. 由于matlab各版本部分语法存在差异,可能会出现bug,用help查帮助文档即可. 如果没有装Matlab,我这里有一篇建模软件的博客:https://www.cnblogs.com/fangxiaoqi/p/10563509.html 变量名:字母数字串(第一个字符必须英文字母 | 字符间无空格 | 最多19个字符): 用%注解:

换了新公司,接触新行业,半路出家,看着别人的代码,看着api慢慢理解. 需求如下:已知两点坐标求距离. 思路,没有,站在同事的肩膀上踩路子,给的这个链接 https://developers.arcgis.com/javascript/latest/api-reference/esri-tasks-GeometryService.html#distance 综合看完,基本调用如下: new GeometryService(url).distance(new DistanceParameters(

本文主要讲解在matlab中实现Linear Regression和Logistic Regression的代码,并不涉及公式推导.具体的计算公式和推导,相关的机器学习文章和视频一大堆,推荐看Andrew NG的公开课. 一.线性回归(Linear Regression) 方法一.利用公式 : function [ theta ] = linearReg() %线性回归. X=[1 1;1 2;1 3;1 4]; %注意第一列全为1,即x0=1,第二列才为x1 Y=[1.1;2.2;2.7;3.