- 博客(6)
- 收藏
- 关注
RSS订阅原创 数值积分复华辛普森、梯形公式与误差定区间块
%JF数值积分，fun为积分函数，例如fun=@(x)3*x，a、b分别为下上限，n为划分区间，c为所选积分方法。其中n可由第二个函数通过设定误差R_n求出满足的n。function [op]=JF(fun,a,b,n,c)if (nargin<5) c='T';endx1=(b-a)/2;x2=(b+a)/2;h=(b-a)./n;Sum=0;Sum_1=0;if c=='T' op=x1*fun(x2);elseif c=='X' op=x1/3*(f...
2020-12-03 10:28:12
1399
原创 Cholesky分解（MATLAB实现）
%Cholesky分解法function [a,L]=m_skysd(A)l=length(A);L=zeros(length(A));a=L;a(1,1)=sqrt(A(1,1)-a(1,1)*L(1,1));L=a;for i=2:l L(1,i)=A(1,i)/a(1,1);endfor k=2:l for i=1:l if(i==k&&i<=k) a(k,i)=sqrt(A(k,i)-a(k,:)*L(:...
2020-08-23 13:55:26
3267
原创 时间序列分析之修正指数模型
function difn=difn(A)difn=[];A=diff(A);[r,l]=size(A);if r==1&&l>=2 for i=1:l-1 difn=[difn A(i+1)/A(i)]; endendif r>1 for i=1:r-1 difn=[difn;A(i+1,:)./A(...
2020-02-25 15:22:32
4683
原创 基于公平分席位求解各单位席位数（MATLAB）
function distribute=distribute(people,seat)%people矩阵：(n,1)n个单位人数%seat:总席位减nuseseat=ones(length(people),1);freeseat=seat;Q=@(people,useseat) people.^2./((useseat+1).*useseat);while freeseat>0...
2020-02-13 17:14:11
320
原创 删除特定后缀特定范围大小文件
该文为python代码#!/usr/bin/env python3#coding: utf-8import os,timedef sizeremove(path,a=0,c=0,last=".abcdef"):#[a,c](kb)之间的文件都被删除 oldpath=path #global total_size命名为全局变量 #path=os.path.abspath(path) ...
2020-02-13 15:05:01
115
原创 MATLAB施密特正交化
施密特正交化的矩阵写法(MATLAB)clearclca=input('所需要正交化向量构成的矩阵');[lklkl,oioi]=size(a);for i=1:oioi if i==1 da=[a(:,i)/norm(a(:,i))]; end if i>=2%将第i个的所有施密特正交化系数写于矩阵b中alpha(i)除外 b=[]...
2019-06-29 14:11:26
6723
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人