• 群智能优化算法测试函数matlab源代码 global M; creatematrix(2); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %画ackley图。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % ackley x from[-5 5] % x=-5:0.01:5; % [x,y]=...
群智能优化算法测试函数matlab源代码

global M;
creatematrix(2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画ackley图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% ackley x from[-5 5]
% x=-5:0.01:5;
% [x,y]=meshgrid(x);
% temp1=x.^2+y.^2;
% temp2=cos(2*pi*x)+cos(2*pi*y);
% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);
% axis([-5,5,-5,5]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的ackley图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%Rotated ackley x from[-5 5
% x=-5:0.01:5;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% temp1=x.^2+y.^2;
% temp2=cos(2*pi*x)+cos(2*pi*y);
% z=20+exp(1)-20*exp(-0.2*sqrt(temp1/2))-exp(temp2/2);
% axis([-5,5,-5,5]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画cigar图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%cigar
x=-5:0.01:5;
[x,y]=meshgrid(x);
z=x.^2+(10^4)*y.^2;
axis([-5,5,-5,5]);
meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的cigar图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%rotated cigar x from [-5 5]
% x=-7:0.01:7;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% z=x.^2+(10^4)*y.^2;
% axis([-7,7,-7,7]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画ellipse图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%ellipse x from [-5 -5]
% x=-5:0.01:5;
% [x,y]=meshgrid(x);
% z=x.^2+10.^4*(y.^2);
% axis([-5,5,-5,5]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画griewangk图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%griewangk
% x=-600:1:600;
% [x,y]=meshgrid(x);
% temp1=(x.^2+y.^2)/4000;
% temp2=cos(x)*cos(y/sqrt(2));
% z=temp1-temp2+1;
% axis([-600,600,-600,600]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的griewangk图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rotated griewangk
% x=-600:1:600;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% temp1=(x.^2+y.^2)/4000;
% temp2=cos(x)*cos(y/sqrt(2));
% z=temp1-temp2+1;
% axis([-600,600,-600,600]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画nocon_rastrigin图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%nocon_rastrigin
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% temp1=0;
% temp2=0;
% z=0;
% if abs(x)<1/2
% temp1=x;
% else
% temp1=round(2.*x)/2;
% end
% if abs(y)<1/2
% temp2=y;
% else
% temp2=round(2.*y)/2;
% end
% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;
% axis([-10,10,-10,10]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的nocon_rastrigin图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rotated nocon_rastrigin
% x=-15:0.05:15;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% temp1=0;
% temp2=0;
% z=0;
% if abs(x)<1/2
% temp1=x;
% else
% temp1=round(2.*x)/2;
% end
% if abs(y)<1/2
% temp2=y;
% else
% temp2=round(2.*y)/2;
% end
% z=z+temp1.^2-10*cos(2*pi.*temp1)+temp2.^2-10*cos(2*pi.*temp2)+20;
% axis([-15,15,-15,15]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画rastrigin图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%rastrigin
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% temp1=x.^2-10*cos(2*pi.*x);
% temp2=y.^2-10*cos(2*pi.*y);
% z=temp1+temp2+20;
% axis([-10,10,-10,10]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的rastrigin图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rotated rastrigin
% x=-15:0.05:15;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% temp1=x.^2-10*cos(2*pi.*x);
% temp2=y.^2-10*cos(2*pi.*y);
% z=temp1+temp2+20;
% axis([-15,15,-15,15]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画rosenbrock图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%rosenbrock
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% z=100*(y-x.^2).^2+(x-1).^2;
% axis([-10,10,-10,10]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画schwefel图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%schwefel
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;
% axis([-10,10,-10,10]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的schwefel图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rotated schwefel
% x=-15:0.05:15;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% z=(x-x.^2).^2+(x-1).^2+(x-y.^2).^2+(y-1).^2;
% axis([-15,15,-15,15]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画sphere图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%sphere
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% z=x.^2+y.^2;
% axis([-10,10,-10,10]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画sumcan图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%sumcan
% x=-0.16:0.0005:0.16;
% [x,y]=meshgrid(x);
% z=1./(10.^(-5)+abs(x)+abs(x+y));
% z=log10(z);
% axis([-0.16,0.16,-0.16,0.16]);
% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画tablet图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%tablet
% x=-10:0.05:10;
% [x,y]=meshgrid(x);
% z=10.^4*x.^2+y.^2;
% axis([-10,10,-10,10]);
% mesh(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画Weierstrass图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Weierstrass
% x=-1:0.005:1;
% [x,y]=meshgrid(x);
% z=0;
% for k=0:20
% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));
% end
% temp=0;
% for k=0:20
% temp=temp+0.5^k*cos(2*pi*3^k*0.5);
% end
% z=z-2*temp;
% axis([-1,1,-1,1]);
% meshc(x,y,z);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%画旋转的Weierstrass图。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%Rotated Weierstrass
% x=-1.5:0.005:1.5;
% [x,y]=meshgrid(x);
% for i=1:size(x,1)
% for j=1:size(y,1)
% p=[x(i,j),y(i,j)]';
% x(i,j)=M(1,:)*p;
% y(i,j)=M(2,:)*p;
% end
% end
% z=0;
% for k=0:20
% z=z+0.5^k*cos(2*pi*3^k*(x+0.5))+0.5^k*cos(2*pi*3^k*(y+0.5));
% end
% temp=0;
% for k=0:20
% temp=temp+0.5^k*cos(2*pi*3^k*0.5);
% end
% z=z-2*temp;
% axis([-1.5,1.5,-1.5,1.5]);
% meshc(x,y,z);
%
xlabel('x');
ylabel('y');
zlabel('f(x,y)');
colormap jet;
% clear;
%saveas(gcf,'temp.jpg')


展开全文
• function [mt1,mt2,upl_mt1,lor_mt1,upl_mt2,lor_mt2,match_weight,err]=temp_temp_matching(t1,t2,upl_1,lor_1,upl_2,lor_2,th,cal)
• 十一、基本文件操作函数 1）文件创建函数 ...tempdir ：返回系统暂存地址名 。 tempname ：返回系统暂存文件名 。 fullfile :创建文件名 2）文件打开、读取、存储函数 open:打开文件 ---用法：op...
十一、基本文件操作函数

1）文件创建函数

filemaker ：把文件名与文件中函数名分开 。

filesep ：文件目录分隔。

fileparts ：把目标文件名拆分成字符串形式输出 。

tempdir ：返回系统暂存地址名 。

tempname ：返回系统暂存文件名 。

fullfile :创建文件名

2）文件打开、读取、存储函数

open:打开文件 ---用法：open('name')

load :将文件导入到MATLABworkspace中。loadfilename load filename X Y Z,

load-asciifilename ,load-mat filename,load('arg1', 'arg2', 'arg3', ...)

save:存储文件.save filename, save filename content, savefilename options

save('filename', 'var1', 'var2', ...)

daqread ：读数据采集工具箱(.daq)中的文件 。

[data, time, abstime, events, daqinfo]= daqread('filename')

uiimport :调用文件导入助手

3)简单文件输入输出函数

fopen ：打开文件。

fid = fopen(filename)

fid = fopen(filename,permission)

Permission值

打开模式

‘r’

只读模式

‘w’

可写入模式，覆盖已有内容

‘a’

可写入模式，写入数据到文件末尾

‘r+’

可读写模式，

‘w+’

可写入模式，覆盖已有内容

‘a+’

可写入模式，写入数据到文件末尾

‘A’

无自动删除添加

‘W’

无自动删除写入

fclose  关闭文件。status=fclose(fid) ,status=fclose('all')

fprintf:以指定格式将数据写入文件 。count = fprintf(fid, format, A, ...)

实例：

 

x=1:0.1:2;

y=[x,sin(x)];

fid=fopen('sin.txt','wt');

fprintf(fid,'%5.2f %5.6f\n',y);

fclose(fid)

fscanf ：以指定格式读取文件中的内容 。A =fscanf(fid, format) , [A,count] = fscanf(fid, format, size)

fread ：读取文件中的二进制数据 。A =fread(fid, count ,precision, skip, machineformat)

fwrite ：将二进制数据写入到文件中 。

fseek ：设置文件位置指针 。status = fseek(fid,offset, origin)

ftell ：获取文件位置指针 。position = ftell(fid)

fwind ：把文件位置指针移动到文件起始位置 。fwind(fid)

fgetl ：按行读取文件内容，并删去换行符 。tline = fgetl(fid)

fgets ：按行读取文件内容，并保留换行符 。tline = fgets(fid, nchar)

frewind ：将文件指针移到文件的起始位置。 frewind(fid)

feof ：测试是否为文件的末尾 。message = ferror(fid)

ferro ：测试文件输入输出是否有错误 。message = ferror(fid,'clear') , [message,errnum] = ferror(...)

disp ：将目标文件所包含的信息以不同形式显示，如文本、数组等形式。

disp(obj)

4)电子表格操作函数

xlsfinfo ：检测文件中是否包含Microsoft Excel电子表格文件 。

typ = xlsfinfo(filename) [typ,desc, fmt] = xlsfinfo(filename)

xlsread ：读Microsoft Excel电子表格文件中的内容 。

num = xlsread(filename)，num = xlsread(filename,sheet)，num =xlsread(filename, range)，num = xlsread(filename, sheet, range, 'basic')

[num,txt,raw,X]=xlsread(filename,sheet, range, 'basic')

xlswrite :将数据写入Microsoft Excel电子表格文件中 。

xlswrite(filename, M)，xlswrite(filename, M, sheet, range)，

[status, message] = xlswrite(filename, ...)

wk1finfo ：检测文件中是否包含Lotus 1-2-3 WK1电子表格。

[extens, typ] = wk1finfo(filename)

wk1read ：读取Lotus 1-2-3 WK1电子表格文件 。

M = wk1read(filename) ，M = wk1read(filename,r,c,range)

wk1write ：将数据写入Lotus1-2-3 WK1电子表格文件中 。

wk1write(filename,M,r,c)

展开全文
• 算法步骤如下： 1、mat2gray是将输入数据F归一... 4、对于数据F中介于最大值和最小值之间的数据是这样处理的：假设temp1为介于最大值和最小值中的某一值，经过mat2gray 变换后变为：temp_last=(temp1-min(F))*det;
• 今天研究C# 调用 matlab 编译的动态库，在将简单加法运算在C# 中调用成功后，满怀开心的心情进一步的调用matlab中regress 函数，结果灾难的报错了，经过半天的奋斗各种百度，谷歌，未能找到解决方案。 .m文件如下：...
• Canny function [ canny ] = canny...temp=rgb2gray(rgb); canny=edge(temp,'canny'); end 灰度 temp=rgb2gray(rgb); 播放视频 clear all source = VideoReader('d:\v\\bus\1.avi'); myObj = VideoWrite...
Canny
function [ canny ] = canny( rgb )
temp=rgb2gray(rgb);
canny=edge(temp,'canny');
end

灰度
temp=rgb2gray(rgb);
播放视频
clear all
source = VideoReader('d:\v\\bus\1.avi');
myObj = VideoWriter('bus1.avi');
open(myObj);
c= source.numberOfFrames;
for i=1:c
f= read(source,i);
f=rgb2gray(f);
f=edge(f,'canny');
writeVideo(myObj,uint8(f));
imshow(f);
if(i>20)
break;
end
end
close(myObj);

类型转换
r1=int16(r1);

二值化
R=im2bw(rs,0.5);

显示图像
figure(2);
imshow(R);

反相，绝对值差
rs=255-rs;
rs=imabsdiff(rs,r2);

读入图片
Imread

旋转
A=imrotate（A，angle,' 旋转实现的方法'，'BBox'）

二值图面积计算
I=imread('bw.bmp');   % 读入二值图像[L,num]=bwlabel(I,8);  %区域标记，       STATS=regionprops(L, 'all');for i=1:num    area(i)=STATS(i).Area;  %计算各区域的面积。end
二值图掩码运算
f1=cat(3,f,f,f);  %前景1维转3维
f1=uint16(f);
r1=uint16(r);
immultiply(r1,f1); %得到运算后的图像

sprintf
for i=1:11
filename=sprintf('f (%d).jpg',i);
t=imread(filename);
t=im2bw(t,0.5);
[L,num]=bwlabel(t,8);
ff(i)=num;
end

灰度图直方图
a=rgb2gray(x33);imhist(a);
转载于:https://www.cnblogs.com/ahuo/p/4089067.html
展开全文
• 代价函数： function [jVal,gradient]=costFunc(theta) %load data fid=fopen('D:\graduate design\data\track2\training.txt'); temp=textscan(fid,'%f%f%f%f%f%f%f%f%f%f%f%f',1000000); temp=cell2mat(temp)...
• 最进看大牛程序有一段 eval(['temp(:,:,ori) = bg_r' num...eval函数是以字符串的形式运行代码下面是matlab官方的注释：eval(expression) evaluates the MATLAB code in the string expression. If you use eval wit...
 最进看大牛程序有一段  eval(['temp(:,:,ori) = bg_r' num2str(r) '{' num2str(ori) '};']);  eval函数是以字符串的形式运行代码下面是matlab官方的注释：eval(expression) evaluates the MATLAB code in the string expression. If you use eval within an anonymous function, nested function, or function that contains a nested function, the evaluated expression cannot create a variable. 其中一个好处就是假如我要对a1,a2,a3,a4,……,a100分别赋予1,2,3，……，100。这时eval就发挥作用了。 for i=1:100     eval(['a' num2str(i) '=' num2str(i)]); end当然还有其他例子，比如批量存数据或图片文件等等。那么开始提到的例子也就好解释了。注意eval中的中括号在两个以上字符串出现时一定要有，起连接作用。如： input：[‘hello’ ‘world’]output：helloworld
转载于:https://www.cnblogs.com/penguins/p/3519361.html
展开全文
•  fs=8000;%采样频率 Fp=1000; Fs=1500; wp=2*pi*Fp/fs;...%line12,窗函数 ...%line13,若引入此行的窗...不调用matlab本身自带的fir函数（因为我必须用c语言设计）</p>
• 原创Matlab产生三角波的函数任意循环-tri_wave.m 自己写的三角波函数，可以方便调用。 参数接口： function y=tri_wave starting_value:三角波起始值 ending_value:三角波终点值 sub_interval:三角波区间...
• matlab三维数组转置函数 我需要读入一个ps0.mat文件，我使用 load ps0.mat ,读到工作空间的temp中，temp的value值是<3 * 640 * 480 unit8> ,我用imshow(temp)却报错，不知道有没有三维数组转置的函数或者别...
• 插值法实现sin函数： %calculate and print the sine function %input: x %output: sin(x) similar function y = sin2(x) %save a copy of x x_temp = x; %calculate the interpolation polynomial %save the ...
• matlab 画图函数，文件的读写   %读取excel数据并画图 A= xlsread('excel路径\excel名','sheet1','a3:b35'); [temp, col] = size(A); plot(A(:,1),A(:,2:col)); xlabel('时间(s)'); legend('销量') F=...
• 关于matlab中bitget函数使用规则的一次小探索1、问题提出2、查阅资料3、硬着头皮的探索4、归纳5、结论 在用matlab进行数字图像处理编程仿真时同学问到了我bitget函数，不巧，对出现的问题我也十分不清楚（问题应该很...
• Matlab数据输入: A = input(“请输入一个大于2的数”);函数： Floor(m) –对m取整数 Tic - toc :测量时间间隔,用于了解程序执行时间。求余： mod 求余函数，mod(X，Y)的符号与Y相同 rem 求余函数，rem(X，Y)的...
• 1. 求实部和虚部 t1=real(temp.');t2=imag(temp.'); 2。 转置 .' 3. circshift 4. fi
• public Matrix<byte> padarray(Matrix<byte> src, int RowPad, int ColPad) { int n = src.Rows;... Mat temp1 = new Mat(n, m + ColPad * 2, DepthType.Cv32F, 1); Matrix<byte>
• temp=157.5; a=temp ./xk; plot(xk,a); %hold on %plot(xk,a); %hold off grid minor; axis([250 1650 0 0.6]); xlabel('5级暗纹的距离'); ylabel('a的距离'); grid minor; 添加细密的网格线 grid on；添加...
• 在学习了线性回归与逻辑...（一） 运用Matlab函数进行回归 %linear regression with matlab inner function x=[1;2;3;4;5;6;7]; y=[2.1;5;5.8;8.2;10.5;11;15]; temp=ones(7,1); X=[temp x]; b=regress(y,X) z=b
• 利用 importdata 函数 data = importdata("123.txt",'',1000) %读取文件，从第1行一直到第1000行 实际应用 data = importdata("123.txt",'',10000) %读取 [m n]= size(data) %寻找含有%的地方 for i =1:...
• ; margin-right:0pt">%绘制调制波形00101010 ... margin-right:0pt">Akk=temp(1:length(Akk));            %码元扩展 ; margin-right:0pt"> </p>
• 实验室的项目中很多要与txt交互，研究所比较喜欢搞这个，有的时候还和一些比较怪的文件格式如.Jin交互，与这些乱七八糟的东西交互与和txt交互差不多，主要的函数有以下几种： temp=dlmread(get(handles.path1,'...
• 返回的则是灰度范围在每份灰度级中的像素个数f = imread('temp.jpg'); g = rgb2gray(f); h = imhist(g, 25);bar（horz， z, width） 绘制条形图horz 表示 水平增量，需要与z的行数相同。 z 表示要显示的条形图的数
• CDF示例代码：cdf.mfunction[xTime,yPercentage]=cdf(initValue,step,endValue,sample); xTime=[]; yPercentage=[]; totalNum=length(sample); fori=initValue:step:endValue temp=length(find(sample&...
• 参考链接：https://www.computationalimaging.cn/2020/01/matlab-train-test-val-inedx-number.html 功能：将[1,N]之间的整数随机抽1%写到val中; 1%写到test中; 其余写到train中 N = 10000; A = [1:N]; A_temp ...
• 我的环境是MATLAB 2016a+vs2015+GPUwindows，想运行一段行人再识别的train源代码，代码是从github上下载的，但是总：显示未定义函数或变量 'net'。如下图所示、 ![图片说明]...
• ; margin-right:0pt">%绘制调制波形00101010 ... margin-right:0pt">Akk=temp(1:length(Akk));            %码元扩展 ; margin-right:0pt"> </p>
• 问题：找出y=x4-4x3+3x+5 （xÎ[0,6]）在[0,6]之间的最小值。 思路：用33位0，1变量来编码x,3位编码整数，30位编码小数。理论上30位编码小数可以将最小值对应的x精确到小数点后九位. ... temp=x; x=y; y=tem...
• 如何读取数据程序进行优化？ 1、预先为变量分配内存； （1）原始程序如下： 程序通过fgetl()逐行读取数据，并对数据进行处理，代码如下： ...temp=d2r/3600; t=[]; freq=100; j=0; fid=fopen(...

matlab 订阅