![](https://img-blog.csdnimg.cn/20201014180756916.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
Matlab && Haskell
文章平均质量分 75
Snail_Walker
Google Machine Learning Tech Lead | 行胜于言
展开
-
Matlab之贝塞尔函数
贝塞尔函数在matlab中已经有了函数,可以直接用。那么我们就没有必要重复发明轮子。以下内容摘自某文档。第一类贝塞尔函数在MatLab中用besselj(NU,Z)来表示:用MatLab的仿真代码是:clear ,clc;format longx=(0:0.01:20)';y_0=besselj(0,x);y_1=besselj(1,x);y_2原创 2014-06-30 01:14:21 · 72641 阅读 · 1 评论 -
Matlab之采样定理
当你看到我的博客的时候说明你已经做到快要吐血了,不过不用担心。可以看看我的,做的不一定对,不过可以借鉴一下! 对连续信号进行等间隔采样形成采样信号,采样信号的频谱是原连续信号的频谱以采样频率为周期进行周期性的延拓形成的。设连续信号的的最高频率为Fmax,如果采样频率Fs>2Fmax,那么采样信号可以唯一的恢复出原连续信号,否则Fs会造成采样信号中的频谱混叠现象,不可能无失真地恢复原连原创 2014-06-30 01:06:39 · 9400 阅读 · 3 评论 -
Haskell之字符串放大
--放大字符串type Picture = [String]--打印出pictureprintPic :: Picture -> IO ()printPic pic = putStr (unlines pic)--map每一个字符,把字符对应成pictureputLetter :: Char -> PictureputLetter s原创 2014-06-30 13:23:14 · 2528 阅读 · 0 评论 -
Haskell之计算毕达哥拉斯三元组
--计算毕达哥拉斯三元组triads :: Int -> [(Int,Int,Int)]triads n = [(x,y,z)|x原创 2014-06-30 13:25:17 · 1472 阅读 · 0 评论 -
Haskell之入门知识
--条件语句:--例子:--求三个数的最大值:maxThree :: Int -> Int -> Int -> IntmaxThree x y z | x>=y&&x>=z = x | y>=z = y | otherwise = z --求两个数的较原创 2014-06-30 00:57:57 · 743 阅读 · 0 评论 -
Haskell之列表的使用
Haskell中没有循环,不过可以用列表实现。--多元组和列表可以进行模式匹配,列表概括,多态函数--可以对多元组进行命名type Purchase = (String,Int)addPair :: (Int,Int) ->IntaddPair (x,y) =x+y--模式匹配name :: Purchase -> Stringname (s,i) =原创 2014-06-30 00:56:58 · 2941 阅读 · 1 评论 -
Haskell之命题与重载
import Data.List-- 定义命题的类型data Prop = Const Bool | Var Char | Not Prop | And Prop Prop | Or Prop Prop | Imply Prop Prop deriving Eq-- 定义p1、p2、p3三种命题p1 = And (Var 'A') (Not (V原创 2014-06-30 13:26:45 · 1603 阅读 · 0 评论 -
Matlab之Gauss-Newton法
这一题就是已知初始点,然后求某一点到三个圆的距离的和的最小值。使用的就是大名鼎鼎的Gauss-Newton方法。这里要求雅克比行列式,由函数 jacobian可以求出。所以这就减少了许多不必要的麻烦。% using Gauss-Newton method% initial point(0,0)% 1、with circul point (0,1),(1,1),(0,-1)% 2、w原创 2014-06-30 01:00:39 · 8632 阅读 · 0 评论 -
Matlab之求导篇
matlab求导命令diff调用格式:diff(函数) , 求的一阶导数;diff(函数, n) , 求的n阶导数(n是具体整数);diff(函数,变量名), 求对的偏导数;diff(函数, 变量名,n) ,求对的n阶偏导数;matlab求雅可比矩阵命令jacobian,调用格式:jacobian([函数;函数; 函数], [])给出矩阵:另外 解微分方程可以用d原创 2014-06-30 01:01:21 · 12598 阅读 · 0 评论 -
Matlab之使用linearization估计人口数量
这种是求解当方程数量多于未知数时,可以使用正规方程来求解。 这种linearization方法是将pt = c1*e^(c2*t),两边取ln得到线性方程组。 线性化之后呢就要进行方程求解,也就是A'Ax = A'b,那么使用x = ((A')*A)\((A')*b);就可以求出x的值,也就是lnc1和c2的值。记得一定要像上面那样写,A'*A在前面,A'*b在原创 2014-06-30 01:02:10 · 1426 阅读 · 0 评论 -
Matlab之求随机生成数字信号单极性不归零码的时域图并求功率谱
》这一题的重点就是首先把码元表示出来,这里我是使用了randi([0,1],1,N);随机生成了0或者1,生产了N个随机数。》接着表示基本码元,这里的是单极性不归零NRZ码元,所以gt1= ones(1,N_sample);%不归零。之后就应用基本码元函数来生成所需要的最后的序列。》再根据所生成的每一个码元把基本码元添加到里面。yt1((i-1)*N_sample+原创 2014-06-30 01:03:39 · 14967 阅读 · 1 评论 -
Matlab之声音处理:对wav音频信号量化
对于matlab,大家应该比较熟悉。今天小哥就带大家来做均匀量化。量化是将模拟信号转化为数字信号必不可少的一步。均匀量化相对比较简单,就是比如将信号在[-1,1]分成相同的很多段,要是某个值在其中某一段,就将其的值归为那一类。每一类都可以用一个唯一的二进制编码表示。现在我们要做的就是对wav的音频信号进行量化。原创 2014-06-30 01:11:09 · 14566 阅读 · 0 评论 -
Matlab之rand函数的使用
rand函数的意义:rand是0-1的均匀分布,randn是均值为0方差为1的正态分布 用法自己help rand或者help randn一下就知道了 rand(n)或randn(n)生成n*n的随机数矩阵 rand(n,m)或randn(m,n)生成m*n的随机数矩阵。 Y = rand(n)Y = rand(m,n)Y = rand([m n])Y = ra原创 2014-06-30 01:04:42 · 8745 阅读 · 0 评论 -
Matlab之画出你的签名
我们来做一些有意思的事情,通过外界的扫描仪或者照相机之类的,获取你写在纸上的签名,图片一般保存为jpg格式。然后用matlab把图片中的签名通过某种采样方法来采取尽可能少的点来用贝塞尔曲线把你的签名勾勒出来并保存为pdf格式。这个东西感觉起来好像很高大上,我们不妨把这个任务分而治之。现在我们需要做的是采集图片->采集样本点->用贝塞尔曲线画出->保存为pdf格式。 高级做法:应用图原创 2014-06-30 01:08:06 · 3899 阅读 · 0 评论 -
Matlab之贝塞尔曲线的逆应用
这一题是已知两个端点以及控制点与端点的斜率和经过的一点。其实也就是要求出两个控制点。这一题很显然已知端点与控制点有垂直切线,所以可得到x2 = -1, x3 = 1。这样就可以把x(t)表示出来了。对于求y(t)的表达式我用到了对称性,所以y2 = y3。然后再用已知点求出x(t)中t的值代入y(t)便可求出y2与y3。这一题便解决了。这一题要注意的是画图的时候不能有变量,必须是确定的。所以可以看原创 2014-06-30 01:08:51 · 1858 阅读 · 0 评论 -
Matlab之用牛顿均差求插值多项式
% page 151 3% 写出一个插值多项式% input :vector(x,y),interpolation point x0% output:the result at x = x0function y0 = page_157_1(x ,y,x0)format longn = length(x);for j = 1:n %fill in y colum原创 2014-06-30 01:12:03 · 9320 阅读 · 0 评论 -
Matlab之与图像有关的函数
函数一:I = imread('path');%读图像函数二:imshow(I);%显示图像函数三:imcontour(I,3);%扣轮廓函数四:figure,imshow(im,[]);title('Raw'); %显示原图im=im2bw(im); %转二值图像figure,imsho原创 2014-06-30 01:07:15 · 1147 阅读 · 0 评论 -
Matlab之将非严格占优矩阵化为严格占优矩阵
这个是只是用行变换将非严格占优矩阵通过行变换转换为严格占有矩阵。伪代码如下:Input matrix:AOutput: If A can transform into a dominance with our method,we can get the dominance. Else output’Cannot transform the原创 2014-06-30 01:12:57 · 2086 阅读 · 0 评论 -
Matlab之三次样条画图和表达式
这一题是得到数据点(0,3),(1,5),(2,4),(3,1)并得到它的三次样条表达式和画出三次样条后的图图形。 以及对数据点(-1,3),(0,5),(3,1),(4,1),(5,1)并得到它的三次样条表达式和画出三次样条后的图图形。 用函数spline可以直接得到,都是如果是要求自然三次样条呢?那就可以在数组y的左右两侧添0。如:csa = spline(原创 2014-06-30 01:13:46 · 9002 阅读 · 0 评论 -
Matlab之用最小二乘建立模型预测值以下例子使用1960,1970,1990和2000的人口估计1980的人口。分别用了直线估计和抛物线估计
以下例子使用1960,1970,1990和2000的人口估计1980的人口。分别用了直线估计和抛物线估计。% page 199 3% this problem is to estimate the populatipn in 1980 and% compare the error using different method to estimate% one use line原创 2014-06-30 01:05:44 · 4057 阅读 · 3 评论 -
Matlab之定位技术DV-HOP的实现
该程序在matlab环境下完成无线传感器中的DV-Hop算法的实现。可算出最短路经算法计算节点间跳数、每个信标节点的校正值并可用用跳数估计距离。%~~~~~~~~~~ DV-Hop算法 ~~~~~~~~~~ ~~~~~~~~~~ ~~~~~~~~~~ % BorderLength-----正方形区域的边长,单位:m% NodeAmount-------网络节点的个数% BeaconAmo转载 2014-11-22 22:10:50 · 7498 阅读 · 18 评论 -
Matlab之DPSK
% communication DPSK demodulation% name:Chen Yu% student ID:12353032function Lab9%% 参数设置%Ts=1/100; %码元周期Rs=1/Ts; %码元频率n1=10;原创 2014-06-30 23:31:45 · 5528 阅读 · 0 评论