- 博客(15)
- 收藏
- 关注
RSS订阅原创 利用偏微分方程学习偏微分方程
clear;N = 7; M = 7; T = 10;c = 2; sigma = 2;dat = DataGen(N,M,T,c,sigma);[a,f] = learner(dat);%%%% we check the solution a going with sigma %%%%for sigma = 0.5:0.1:3 dat = DataGen(N,M,T,c,sigma); [a,f] = learner(dat); disp(sum(sum(a)));endclea
2022-04-02 22:31:33
181
原创 图像处理的傅里叶级数方法matlab代码
clear;P = imread('pic.jpg');I = rgb2gray(P);I = double(I);s = size(I);N_ = s(1); M_ = s(2);imwrite(uint8(I),'I0.jpg');N = 50;Ccc = zeros(N+1,N+1);Ccs = zeros(N+1,N+1);Csc = zeros(N+1,N+1);Css = zeros(N+1,N+1);for E = 0:N for n = 0:E m = E
2022-03-04 19:47:07
2681
原创 方程数值解不稳定的算例
解方程{di+xi+λϕ(xi)Φ(xi)=0,i=1,⋯ ,p∏i=1pΦ(xi)=r\left\{\begin{aligned}&d_i+x_i+\lambda\frac{\phi(x_i)}{\Phi(x_i)}=0,i=1,\cdots,p\\&\prod_{i=1}^p\Phi(x_i)=r\end{aligned}\right.⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧di+xi+λΦ(xi)ϕ(xi)=0,i=1,⋯,pi=1∏pΦ(xi)=r其中rrr,did_
2020-08-13 11:29:27
370
原创 GLM的MLE的渐进正态性
(GLM β\betaβ 的MLE的渐近正态性的证明）设 y∽exp(yTθ−b(θ))y\backsim \exp(y^T\theta-b(\theta))y∽exp(yTθ−b(θ)), θ=u(xTβ)\theta = u(x^T\beta)θ=u(xTβ)Score nSn=1n∑j=1nxju.(xjTβ)(yj−h(xjTβ))\sqrt{n}S_n = \frac{1}{\sqrt{n}} \sum\limits_{j=1}^n x_j \overset{.}{u}(x_j^T\bet
2020-06-02 17:30:26
1188
原创 重对数律的R语言模拟
Law of Iterated Logarithm for Weiner Processlim supt→∞∣Bt∣2tloglogt=1\limsup_{t\to\infty}\frac{|B_t|}{\sqrt{2tloglogt}}=1t→∞limsup2tloglogt∣Bt∣=1n = 40 #模拟次数N = 100000 #一次模拟中的长度func = functi...
2020-02-20 22:01:05
294
原创 GTK+Linux连连看游戏
init.c#include"LTA.h"#include<stdio.h>void table_init(){ int i, tmp1, tmp2, x1, y1, x2, y2; for(i=0; i<TABLE_COL+2; i++) _table[0][i] = _table[TABLE_ROW+1][i] = 1; for(i=0; i<T...
2020-02-16 21:10:54
357
原创 Backward nonlinear diffusion
pic0 = imread('pic0.jpg');f = imresize(pic0, 0.4);f = rgb2gray(f);s = size(f); N = s(1); M = s(2);for i = 1:10 f = convolute(f, N, M, 1);endimwrite(uint8(f),'pic.jpg');pic = imread('pic.jpg')...
2020-02-16 20:56:04
169
原创 Nonlinear isotropic diffusion filtering
∂tu=div(D(∇uσ)∇u)\partial_{t} u=\operatorname{div}\left(D\left(\nabla u_{\sigma}\right) \nabla u\right)∂tu=div(D(∇uσ)∇u) D(∇uσ)=g(∣∇uσ∣)ID\left(\nabla u_{\sigma}\right)=g\left(\left|\nabla u_{\si...
2020-02-12 14:16:48
298
原创 Nonlinear diffusion by PDE in image processing
∂u∂t=div(c(∣∇u∣2)∇u)\frac{\partial u}{\partial t}=div(c(|\nabla u|^2)\nabla u)∂t∂u=div(c(∣∇u∣2)∇u)当c=1c=1c=1时是热方程。Images are piecewise smooth functions. We want to smooth the image while preserving...
2020-02-09 16:21:39
162
原创 attenuating high frequency by Fourier transform
Fourier平面距(0,0)(0,0)(0,0)较远的点是高噪声，这Fourier变换的形式可以看出F(k1,k2)=∫f(x,y)exp(−ik1x−ik2y)F(k_1,k_2) = \int f(x,y)exp(-ik_1x-ik_2y)F(k1,k2)=∫f(x,y)exp(−ik1x−ik2y)较大的(k1,k2)(k_1,k_2)(k1,k2)是fff向较大的sin...
2020-02-09 07:47:40
61
原创 热方程与光滑处理图片
最简单的smooth一张图片的方法是热方程ut−a2Δu=0u(0,x,y)=f(x,y)u_t - a^2\Delta u=0\\u(0,x,y) = f(x,y)ut−a2Δu=0u(0,x,y)=f(x,y)To explain this phenomenon, on the one hand, 数值解u(x,y)=u(x,y)+(a^2)*(u(x+1,y)+u(x-1,y)+u(...
2020-02-07 14:38:53
312
1原创 摆三游戏成功率的随机模拟
# cmd: C:\python.lnk D:\simulator.py > D:\1.txt# cmd Ctrl+C 退出当前程序class node(object): def _init(self, val = 0, prev = 0, next = 0): self.data = val self.prev = prev self.next = nextdef i...
2020-02-03 23:02:09
133
原创 Finding edges in images 和 decomposing images into geometry and texture
Finding edges 是 image processing 成熟的理论MATLAB自带edge函数% MATLAB自带的edge函数pic = imread('pic.jpg');tmp = size(pic);N = tmp(1); M = tmp(2);pic_ = rgb2gray(pic);pic2 = edge(pic_,'canny');for x = 1:N ...
2020-02-02 21:25:28
159
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人