代码搜索:动力学模型
找到约 5,890 项符合「动力学模型」的源代码
代码结果 5,890
www.eeworm.com/read/491628/6431598
m globaladaptive.m
%文件名:globaladaptive.m
%程序员:郭迟
%编写时间:2004.4.7
%函数功能:本函数将完成Watson模型下全局自适应水印
%输入格式举例:[result,alf]=globaladaptive('c:\lenna.jpg',0.1,1983,'c:\wm.jpg')
%参数说明:
%image为输入原始图像
%equal为期望的感知质量
%seed为随机数
www.eeworm.com/read/491628/6431618
m watsondistorsion.m
%文件名:watsondistorsion.m
%程序员:郭迟
%编写时间:2004.4.2
%函数功能:本函数将完成Watson模型下图像感知质量度量.
%输入格式举例:result=watsondistorsion('c:\lenna.jpg','c:\test.jpg')
%参数说明:
%original为原始图像
%test加有水印的图像
%result为两者误差估计
f
www.eeworm.com/read/484469/6577401
txt prony1.txt
x=0:0.05:10;%原理想信号
y=220*exp(-0.25*x).cos(2*pi*1.5*x+pi)+110*exp(0.25*x).cos(2*pi*10*x+pi/2);
p=6% p为模型阶数
sf=1/0.001 %sf为采样频率
[F,D,A,theta]=exprony_m a(y,p,sf) %调用prony算法函数
function [F,D,A,theta
www.eeworm.com/read/484393/6587390
m ar2.m
function [a, e] = AR2(x, p)
% AR自相关算法谱估计,输入采样数据x,预测阶数p,输出模型传输函数系数a
% 返回H(z)=e/(1+A(z))中的A,e是sigma^2
N = length(x);
A = zeros(p, N+p);
for n=1:p
for m=1:N+p
if m=N+1
www.eeworm.com/read/484393/6587394
m ar.m
function [a, e] = AR(x, p)
% AR自协方差算法谱估计,输入采样数据x,预测阶数p,输出模型传输函数系数a
% 返回H(z)=e/(1+A(z))中的A,e是sigma^2
N = length(x);
A = zeros(p, N-p);
for n=1:p
for m=1:N-p
A(n, m) = x(p-n+m); %A'
www.eeworm.com/read/484398/6587415
m ar2.m
function [a, e] = AR2(x, p)
% AR自相关算法谱估计,输入采样数据x,预测阶数p,输出模型传输函数系数a
% 返回H(z)=e/(1+A(z))中的A,e是sigma^2
N = length(x);
A = zeros(p, N+p);
for n=1:p
for m=1:N+p
if m=N+1
www.eeworm.com/read/484398/6587419
m ar.m
function [a, e] = AR(x, p)
% AR自协方差算法谱估计,输入采样数据x,预测阶数p,输出模型传输函数系数a
% 返回H(z)=e/(1+A(z))中的A,e是sigma^2
N = length(x);
A = zeros(p, N-p);
for n=1:p
for m=1:N-p
A(n, m) = x(p-n+m); %A'
www.eeworm.com/read/263879/11338078
m jm99asmfun.m
function f=test(data,n,m)
%年中国大学生数学建模竞赛a题:自动化车床管理模型一
%参见《数学的实践与认识》2000.1.p36-40
%费用函数
%n--examing interval
%m--Max number of products
%x--Living
%y--average loss
x=data;
num_prod=min((n*ceil
www.eeworm.com/read/263879/11338098
m jm99a1.m
%99年中国大学生数学建模竞赛a题:自动化车床管理模型一
%参见《数学的实践与认识》2000.1.p36-40
clear;opt=inf;
for n=[10 20 30]
for m=[200 300 400]
[x0,op]=fmins('jm99afun',[n m]);
if op(8)
www.eeworm.com/read/263879/11338270
m quadeg5.m
%微积分例5:微分方程(rk4,ode45)
%需用模型函数quadeg5fun.m
clear;close;
[t,y]=rk4('quadeg5fun',[0,7],1,0.03);
plot(t,y);
hold on;
[t,y]=ode45('quadeg5fun',[0,7],1);
plot(t,y,'r:');
hold off;