注册 登录
编程论坛 Matlab

MATLAB新手快来 新手问题交流帖

wyy022 发布于 2006-11-20 17:45, 12555 次点击
我也是个MATLAB新手,我也是个乐于助人的青年,


你是不是和我一样呢?


进来吧,发表你不懂的问题,讲讲你的看法



这里是你我的舞台,再幼稚的问题也要有勇气说出来



来吧,说出你的问题。这里的每个新手 都会用心的回答你的每一个问题。

[此贴子已经被作者于2006-11-20 17:49:53编辑过]

87 回复
#52
wgm_cwang2007-04-12 16:38
回复:(wyy022)MATLAB新手快来 新手问题交流帖
想要选出方波(矩形波)形状的图形,已有很多规一化的图形(如一万个),每个图形由一组数组(如200个元素)勾画而成,并且数值中每一个值都已知,那么怎么样规定条件使得选出的图形比较符合方波要求,即通带比较平坦,陡度越大越好。
我选阈值0.9和0.1,当数组中的具体值大于0.9或者小于0.1的个数够大,就选那些图形。可是选出的图陡度是够好,通带不平(在大于0.9和小于0.1时候),不太合乎方波要求。是不是可以比较每个数组中相邻两个元素值相差的大小来选?不是很大就好?但是这样得出的图会不会是正弦波图啊?
请问各位,规定什么样的条件能够筛选出比较符合方波要求的图。
急!
先谢谢各位大侠啦!
#53
hujiang2007-04-14 12:03
我有一系列的Q1.dat,Q2.DAT,......数据文件,格式如下
1 0.5 1.2
2 1.0 1.05




N 5000 1.25
请问如何把文件连续读入matlab,并且以第二列和第三列作图!不知哪位可否告知!
#54
yaohaiyan2007-04-19 09:37
请问有谁会编用matlab
#55
yaohaiyan2007-04-19 09:38
请问有谁会编用matlab实现混沌同步的程序吗?
#56
tommy07232007-04-22 00:09
微分方程求解并做图,恳请高手帮忙,小弟不胜感激~~~
dN1/dt=[r1*(1-N1/K1)-m*N2]*N1;
dN2/dt=[r2*(1-N2/K2)-n*N1]*N2;
其中r1=0.01,r2=0.05,K1=1000,K2=1000,m=9e-6,n=4e-5;
N1(0)=0,N2(0)=0.
#57
wyy0222007-04-24 15:17
回复:(yaohaiyan)请问有谁会编用matlab实现混沌同步...

我这里可以下载

信息隐藏技术实验教程

[推荐]信息隐藏技术实验教程

作者: 王丽娜 郭迟 李鹏编著


需用超星阅读器

这本书第二章有这方面内容,你可以看看.

我的网盘 http://tel01.ibmdisk.com/?wyy0313

欢迎来数字水印论坛一看.

#58
wyy0222007-04-24 15:20
以下是引用yaohaiyan在2007-4-19 9:37:23的发言:
请问有谁会编用matlab

对了,这本书里是实现混沌随机序列的,不知和你说的是一个东东吗?

#59
du2192007-04-24 19:25
有没有关于周期性脉冲函数构造的程序?
#60
wyy0222007-04-25 13:11

里面有伪随机发生器的构造.

#61
f1111aaaa2007-04-26 10:52
fsolve具体怎么用啊
#62
skyfeng7092007-04-27 16:35

怎样做多维g=f(x,y,z)的三次样条插值呢?

#63
bigspring2007-05-02 09:28

帮忙介绍一下matlab中的ode45或ode23函数,
谢谢

请发到我的信箱bigsping@163.com
十分感谢!!!

[此贴子已经被作者于2007-5-2 9:30:00编辑过]

#64
leonelll2007-05-05 22:08
大牛请帮助解决以下这道题目:
允许输入1~2个参数
如输入一个参数,则判断该参数是否为矩阵:如是,则第一行的最大值,最小值、算术
平均值、中值,并将结果以行相量的形式返回。
如为标量则返回该标量
如输入为2个参数,则第一个参数为矩阵,第二个参数为行号,并求矩阵的相应行的最大值,最小值、算术平均值、中值,并将结果以行相量的形式返回。
主要问题是参数的判断与提取,多谢了
          
                              菜鸟师弟
#65
app21701appp2007-05-06 17:10
我的毕业设计是做:基于matlab的InSAR复图像强度图的生成
哪位师兄师姐能告诉我怎样用MATLAB读取二进制数和ASCII码啊?
以前老师上课讲MATLAB时只是讲得很浅,而我自己看了这方面的书,只是看了一些怎样导入数据。
另外,我那些数据是在EnviView那个软件上的。而那个软件我从来没用过,而且是英文的。
所以诸多难处真的是把我急死了。
哪位大哥大姐能给我一些指导或建议,诸如看MATLAB哪方面的书或是给我一个小例子,读一个二进制数据和ASCII码数据的例子。让我借鉴一下。谢谢你们了!!!!!!!!!!!!!!!我的毕业设计5月8号就要做出基本的程序了。所以我很急。谢谢你们了!!!!!!!!!!!!!!!!
#66
qqxutianchen2007-05-07 19:09
神经网络的matlab实现的相关程序会啊 ?/我急用啊 主要就是画图吧,一个城市的图,高手来帮帮我啊....
#67
fajessmi2007-05-08 21:55
有谁可以用帮我解一下这个积分quad('-40sinx',y,y+2pi),由于上下限也是有参数的,导师说只要用循环处理下,可还是不会编,麻烦帮我编一个计算过程,我的邮箱是fajessmi@yahoo.com.cn能发给那是我最好了,谢谢各位
#68
appleblue2007-05-10 14:15

搭建的电路为什么报错? 高手帮帮忙!!

只有本站会员才能查看附件,请 登录

#69
appleblue2007-05-10 14:19
matlab里面有没有集成运放? 为什么找不到?
#70
navyking2007-05-12 01:45
求助!!如何用simulink仿真自整角机的方位信息的传送?
#71
夜幽书生2007-05-12 10:45
各位高手:
  dx/dt=-0.15*y(t)
dy/dt=-0.1x(t)
x(0)=10000;y(0)=5000
这个微分方程用MATLAB怎样求解啊!

#72
lukaifengdd2007-05-12 16:26
请教
系统提示:
??? Error using ==> mpower
Matrix must be square.
那我问题到底是出在哪儿啊?
#73
wenzi11202007-05-13 17:24

有谁能够写出PCL-818数据采集卡的采集数据(温度)的原代码.
我急需!!!!

#74
果肉2007-05-15 00:24
问个比较简单的问题
∫t|e(t) |dt
现需要用MATLAB 实现上面方程的求解.
积分上下限为0~无穷

我现在编了
syms s a b c
x=ilaplace((s+c)/(s+a)/(s+b));
x=simplify(x)

请问后面要怎么遍?

导师说 其实e(t)就等于C(T)-1
C(T)就是传递函数C(S)的拉是逆变换
按照老师的意思我编了上述程序
后面怎么遍下去就不知道了
请高手指教

[此贴子已经被作者于2007-5-15 0:45:50编辑过]

#75
para2007-05-20 10:45

程序编出来了,调试怎么也过不了,希望哪位大虾帮忙调试一下!
u,v,ed,E矩阵需要预定义吗?怎么预定义?


tic;
clear;clc;
epsilon=8.9e-12; num=21; a=0.2/(num-1); t=(num+1)/2;cons=300;
dp1=0.1;pg=1.025;ki=2.2e-4;
v(t,t)=6e4;v1=u;N=0;
while max(max(abs(v-v1)))>cons;
N=N+1;
v1=v;
for i=2:num-1
for j=2:num-1
Eox=-(v(i-1,j)-v(i+1,j))/(2*a);
Eoy=-(v(i,j-1)-v(i,j+1))/(2*a);
E(i,j)=sqrt(Eox.^2+Eoy.^2);
u(i,j)=0.25*(v(i,j-1)+v(i,j+1)+v(i-1,j)+v(i+1,j)+a^2*ed(i,j)/epsilon);
ed(i,j)=-epsilon/(2*a)*(Eox+Eoy)-0.5*((epsilon/a*(Eox+Eoy))^2
+4*(epsilon/a*(Eox*ed(i+1,j)+Eoy*ed(i,j+1))))^0.5;
u(1,:)=0;u(num,:)=0;u(:,num)=0;u(t,t)=6e4;ed(t,t)=3e-5;
end;
end;
v=u;
end;
d=[0:0.2/(num-1):0.2];
[x,y]=meshgrid(d);
u=abs(u);ed=abs(ed);E=abs(E);
figure(1);surfl(x,y,abs(u));
figure(2);surfl(x,y,abs(ed));
figure(3);surfl(x,y,abs(E));
toc

[此贴子已经被作者于2007-5-20 15:57:07编辑过]

#76
app21701appp2007-05-20 14:42
clear all
fid=fopen('E:ASA_IMS_1PNPDE20040107_061300_000000162023_00120_09693_0012.dat','r')
sta=fseek(fid,1100*1100,-1 )
Array=fread(fid,10*10,'uint8')
上面的3,4行语句出错,能告诉我为什么吗?
谢谢了。我刚学MATLAB。
FSEEK函数倒底怎么用?
#77
2007-05-27 16:08
今天刚接到的课题```望达人帮忙解答下```
<用MATLAB实现头表地形图显示算法>
#78
dmxl19842007-05-31 15:17

我是一个新人,不怎么会编程,下面这个程序的错误在什么地方,帮忙看下
p(1)=0;
for t=2:512
p(t)=0.99*p(t-1)+sin^2*(800*t);
end
plot(p(t))
运行后提示是
??? Error using ==> sin
Not enough input arguments.

Error in ==> sin at 14
[varargout{1:nargout}] = builtin('sin', varargin{:});

Error in ==> l at 3
p(t)=0.99*p(t-1)+sin^2*(800*t);

#79
qinyanfeng002007-05-31 20:23
以下是引用woshiheshibo在2006-11-26 13:15:54的发言:

你好,能解释一下以下gui个部分的功能吗


function varargout = eggui3(varargin)
% EGGUI3 M-file for eggui3.fig
% EGGUI3, by itself, creates a new EGGUI3 or raises the existing
% singleton*.
%
% H = EGGUI3 returns the handle to a new EGGUI3 or the handle to
% the existing singleton*.
%
% EGGUI3('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in EGGUI3.M with the given input arguments.
%
% EGGUI3('Property','Value',...) creates a new EGGUI3 or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before eggui3_OpeningFunction gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to eggui3_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES

% Edit the above text to modify the response to help eggui3

% Last Modified by GUIDE v2.5 25-Oct-2006 20:33:17

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @eggui3_OpeningFcn, ...
'gui_OutputFcn', @eggui3_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end

if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT


% --- Executes just before eggui3 is made visible.
function eggui3_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to eggui3 (see VARARGIN)

% Choose default command line output for eggui3
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% UIWAIT makes eggui3 wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = eggui3_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)

% Get default command line output from handles structure
varargout{1} = handles.output;


% --- Executes during object creation, after setting all properties.
function figure1_CreateFcn(hObject, eventdata, handles)
% hObject handle to figure1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called

function zeta_edit_Callback(hObject, eventdata, handles)
% hObject handle to zeta_edit (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of zeta_edit as text
% str2double(get(hObject,'String')) returns contents of zeta_edit as a double
z=str2num(get(handles.zeta_edit,'string'));
t=0:0.1:15;
cla;
%plot(t,sin(t));
for k=1:length(z)
y(:,k)=step(1,[2*z(k),1],t);
line(t,y(:,k));
end

% --- Executes during object creation, after setting all properties.
function zeta_edit_CreateFcn(hObject, eventdata, handles)
% hObject handle to zeta_edit (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end


% --- Executes on button press in gridon_push.
function gridon_push_Callback(hObject, eventdata, handles)
% hObject handle to gridon_push (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
grid on;

% --- Executes on button press in gridoff_push.
function gridoff_push_Callback(hObject, eventdata, handles)
% hObject handle to gridoff_push (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
grid off;

% --------------------------------------------------------------------
function option_Callback(hObject, eventdata, handles)
% hObject handle to option (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)


% --------------------------------------------------------------------
function Box_on_Callback(hObject, eventdata, handles)
% hObject handle to Box_on (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)


% --------------------------------------------------------------------
function Untitled_2_Callback(hObject, eventdata, handles)
% hObject handle to Untitled_2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
box on;

% --------------------------------------------------------------------
function Box_off_Callback(hObject, eventdata, handles)
% hObject handle to Box_off (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
box off;

那些前面 % 的注释部分不要理会,没 % 的是你设计的GUI图形设计中控件的函数标题,在一个模块下,如:
% --------------------------------------------------------------------
function Box_off_Callback(hObject, eventdata, handles)
% hObject handle to Box_off (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
box off;
你可以编程 来实现这个控件的功能。
所有的控件都能实现的话,就可以组成一个GUI界面,就象VFP做的那样一个界面。
上学期,我们做DSP课设的时候就是我就是用这个GUI设计一个高、低、带通滤波器的界面,效果还不错。
#80
zz060032322007-05-31 20:59

各位高手.小弟我是新手入门了.想请教下.
如果我在matlab 编程中需要读取一些数据,,我要的只是数据,而不要用到前面的文件头啊..所以想在编程的时候在读取数据的时候把这些自动去掉...该怎么做啊~~(像这样的数据不只是一个..而是有很多个,所以手动删取是不可以的)这些数据是是如下:

#[FileName]
#FileName=D:\data\ww\0527\bk_40000.arc_data
#Map FileName=D:\data\zgp\526\2-2.arc_map
#Map Area Number=1
#TimeStamp=05-5-27 15:33:13
#[CCD]
#Interface=PCI Timer Card
#Controller=ST133
#CCDHead=HAMAMATSU 128 x 1024
#Mono_1_Present=true
#Mono_2_Present=false
#FilterWheel_1_Present=false
#FilterWheel_2_Present=false
#[CCD_Acquisition]
#TimeTaken=05-5-27 15:32:58
#DataCalc=ch1
#WaveUnits=nm
#Cycles=1
#DataSource=Area_Map_Profile
#IntegrationTime=40000
#HeadTemp=-34.5
#TempLocked=0
#AtoD=slow
#[CCD_Strips]
#Number_of_Strips=1
#[CCD_Strip_1]
#Startx=1
#Endx=1024
#Starty=1
#Endy=128
#[CCD_Mono3080650]
#Model=SP-308
#Serial=3080650
#WaveLength=5400000
#Grating=2
#Grating_Density=1200
#Grating_Blaze= 500 nm
#Entry=side
#EntryWidth=10
#Exit=front
#ExitWidth=0

# X Y
510.5755008 2842
510.6316334 2895
510.6877639 3008
510.743892 3035
510.800018 3030
510.8561416 3055
510.912263 3230
510.9683822 3211
511.024499 3273
511.0806137 3266
511.136726 3388
511.1928361 3307
511.248944 3328
511.3050495 3473
511.3611528 3411
511.4172539 3418
511.4733526 3396
511.5294491 3376

#81
jiangbao2007-06-01 11:59

我想求助一个matlab程序
是有关一种图像压缩编码的方法——算术编码的编码和解码过程的编程
各位高手帮帮忙呀
我后天就要交了,很急呀
拜托了

#82
zp59462007-06-04 10:29
我也是新手 可是怎么就是找不到学习语言的感觉啊 恳请各位大虾指点迷津
#83
weiyue5142007-06-08 23:44
请问怎么在离散的点图中显示按单位长度显示网格??
grid是默认显示无法按单位长度显示!!!
拜求答案!
急~
#84
leki2007-06-09 00:35
回复:(weiyue514)请问怎么在离散的点图中显示按单位...

看了一下grid函数
想到一个方法是可以自己设置显示的坐标,如下
t=0:pi/20:2*pi;
h1=plot(t,2*sin(t))
grid on

只有本站会员才能查看附件,请 登录

自己用下面的语句设置x轴坐标后
set(gca,'xtick',[1 3 5 7])

只有本站会员才能查看附件,请 登录


#85
leki2007-06-09 00:42
回复:(zz06003232)各位高手.小弟我是新手入门了.想...

内容格式

你好
1 11 111 1111
欢迎来到
2 22 222 2222
动力学与仿真控制论坛
3 33 333 3333
www.dytrol.com
4 44 444 4444
5 55 555 5555

读取方法

以下是由chinamaker编写的一种方法:
fidin=fopen('test2.txt'); % 打开test2.txt文件
fidout=fopen('mkmatlab.txt','w'); % 创建MKMATLAB.txt文件
while ~feof(fidin) % 判断是否为文件末尾
tline=fgetl(fidin); % 从文件读行
if double(tline(1))&gt;=48&&double(tline(1))&lt;=57 % 判断首字符是否是数值
fprintf(fidout,'%s\n\n',tline); % 如果是数字行,把此行数据写入文件MKMATLAB.txt
continue % 如果是非数字继续下一次循环
end
end
fclose(fidout);
MK=importdata('MKMATLAB.txt'); % 将生成的MKMATLAB.txt文件导入工作空间,变量名为MK,实际上它不显示出来
MK

MK =

1 11 111 1111
2 22 222 2222
3 33 333 3333
4 44 444 4444
5 55 555 5555

自己研究一下上面的内容吧

#86
leki2007-06-09 00:46
回复:(zp5946)我也是新手 可是怎么就是找不到学习语...

个人认为新手就是要多练习,我也是一样的才学不久
当然前提是一些基本的语法一定要知道

#87
zhfl13142007-06-09 23:23

如何用内插函数sinc函数对已经产生的序列进行从构函数啊??????

#88
wubotao10072007-06-13 02:56
我在毕业设计中 各位帮忙下
我程序是这样的
function f=fun(x)
global K L thetam ax alpha
for i=1:61
f=0
betae=atan(tan(alpha(i)./1-(K/L)*tan(alpha(i))));
A(i)=2*x(1).^2*sin(x(2)+alpha(i));
B(i)=2*K*x(1)-2*x(1).^2*cos(x(2)+alpha(i));
C(i)=2*x(1).^2-4*x(1).^2*(cos(x(2)).^2-4*K*x(1)*cos(x(2))-2*K*x(1)*cos(x(2)+alpha(i)));
theta3=2*acot((A(i)+sqrt(A(i).^2+B(i).^2-C(i).^2))/B(i)+C(i))
beta(i)=x(2)+theta(i0-pi)
if alpha(i)<=pi/18
f(i)=1.5*abs(beta(k)-betae(k));
elseif alpha(i)<=pi/180;f(i)=abs(beta(k)-betae(k));
else f(i)=0.5*abs(beta(k)-betae(k));
end
end
global K L thetamax alpha
K=input('输入两中心线间距(单位mm)K=');
L=input('输入轴距(单位mm)L = ');
thetamax=input('输入外转向轮最大转角(单位为度) ');
x0(1)=input('输入初始点的第1 个分量(臂长、mm ) ') ;
x0(2)=input('输入初始点的第2 个分量(底角、度) ') ;
thetamax=thetamax*pi/180;
x0(2)=x0(2)*pi/180;lb(1)=0.11*K;lb(2)=acot(K/(1.2*L));ub(1)=0.17*K;ub(2)=pi/2;
alpha=linspace(0,thetamax,61);
lb=[lb(1),lb(2)];
ub=[ub(1),ub(2)];x0=[x0(1),x0(2)];
options=opimset('To1Fun','le-10','To1-Con','le-6')
[x,resnonn]=Isqnonlin('fun',x0,lb,ub,options)
global K L thetam ax alpha
K= input('输入两主销中心线间距(单位mm )K= ') ;
L = input('输入轴距(单位mm)L = ');
thetamax=input ('输入外转向轮最大转角(单位为度) ') ;
x (1) =input ('输入优化解的第1 个分量(臂长, 单位为mm ) ') ;
x (2) =input ('输入优化解的第2 个分量(梯形底角, 单位为度') ;
thetamax=thetamax*pi/180;  
x(2) = x(2)*pi/180;
alpha= lin space(0, thetam ax, 61) ;  
for i= 1:61
betae= atan ( tan (alpha (i) )/(1- K/L ) *tan(alpha ( i)));
A ( i) = 2* (x (1) ).^ 2*sin (x (2) + alpha( i) ) ;   
B ( i) = 2* K* x(1) - 2* (x(1) ). ^2*co s (x(2) + alpha ( i) ) ;
C ( i) = 2* (x (1) ). ^ 2- 4* (x (1) ). ^ 2*(co s (x (2) ) ). ^2+ 4* K* x (1) * co s (x (2) ) -2*K* x (1) *co s (x (2) + alpha ( i) ) ;
theta3= 2*aco t (A ( i) )+ sqrt ( ( (A ( i) ). ^2+(B ( i) ). ^2- * (C ( i) ). ^2) )/B ( i) + C ( i) ) ) ;
beta ( i) = x (2) + theta3 ( i) - p i;
end
p lo t (alpha ( i) , betae ( i) ,‘r’) ; ho ld off
12