求MATLAB大神帮我看看程序哪有问题，谢谢！_百度知道
求MATLAB大神帮我看看程序哪有问题，谢谢！
1.04;((x+0.*exp(-4; index out of bounds because numel(y)=1;
g(n;end这个程序运行完屏幕上显示&gt，拜托大家帮忙;
r22=Rs&#47.5*O*(R/
%D的取值区间上限,m)=quad(f！function g=f()clear all.5; f at 31
g(n:c2,y);n=0.002*y1)^2+(R+2*r21)^2);
r11=Rs*(Rs+P21)&#47.5*(-(2*Rs+P21+R)+sqrt(D)); quad at 71if ~isfinite(y(7))Error in ==&
D=4*Rs^2+R^2+4*Rs*P21+2*P21*R+P21^2;for O=c1;t=0;b*2+1));syms y1.002*y1)^2+(R&#47:b:a
m=m+1.1;sqrt(D)*(1/
%O取值区间下限;2+r21)^2),(a/
n=n+1;sqrt(D)*(exp(s1*t)-exp(s2*t)).Error in ==&&gt? Attempted to access y(7):b;
s2=0,y1,100),0.5*(-(2*Rs+P21+R)-sqrt(D));c1=0,Z);
endendy=c1;a=2;d=0.1?;(4*(x+0,y?;R=0;(Rs+s2)*(exp(-Rs*t)-exp(s2*t))-1/2+r21)*(r11-r22)&#47:c2
if (O-fix(O)==0)
clear maplemex,m)=quad(f;
%D的步长.^2;;surf(x:a:d,0.4*10^(-8)*y1;
%O取值区间上限;b=0:d;x=-a;c2=10;
Rs=(O^2*(R+2*r21))/g=zeros(((c1-c2)/
%得到被积函数
f=@(y)subs(f.83*10^(-6)*y1^2)*r12;
m=0;P21=1;
f=2;r21=0;(Rs+s1)*(exp(-Rs*t)-exp(s1*t)))+exp(-Rs*t);
r12=0;d+1);
s1=0搞了一下午了;Z=g,100)
每循环一次，y就会被赋值一次，怎么会出现y长度不够的问题呢？
提问者采纳
surf(x. 下面的代码是没改的;&nbsp.*exp(-4, 不过你画出来应该是类似的;Z=g;%O取值区间下限;%得到被积函数f=@(y1)subs(f)+0.04;&nbsp:b.1;c1=0,Z):a,100);D=4*Rs^2+R^2+4*Rs*P21+2*P21*R+P21^2, 仅仅是更精细.a=2;r12=0我把步长改了改;for&%&y1.5*(-(2*Rs+P21+R)+sqrt(D));结果肯定不受影响;&&nbsp.r21=0;Rs=(O^2*(R+2*r21))/(4*(x+0.jpg" target="_blank" title="点击查看大图" class="ikqb_img_alink"><img class="ikqb_img" src="n=0.4*10^(-8)*y1,y;(O-fix(O)==0)clear&nbsp:d, 因为不然电脑运行太慢不出结果;syms&&nbsp,0;%O的步长.002*y1)^2+(R/2+r21)^2).^2;b=0;endendy=c1:am=m+1.jpg" esrc="/zhidao/wh%3D450%2C600/sign=ef33a879e3c0e/a08b87dfcb30e924b899f399;r11=Rs*(Rs+P21)/sqrt(D)*(1/(Rs+s2)*(exp(-Rs*t)-exp(s2*t))-1/(Rs+s1)*(exp(-Rs*t)-exp(s1*t)))+exp(-Rs*t);%D的取值区间上限;maplemex.*y1.1:c2,&g=zeros(((c1-c2)/d+1);d=0;就能运行了, 所以改长了.5.1.*r12;/zhidao/wh%3D600%2C800/sign=c77314a1cdea15ce41bbe80f863016cb/a08b87dfcb30e924b899f399.5*O*(R/2+r21)*(r11-r22)/((x+0;&%ρ11r22=Rs/sqrt(D)*(exp(s1*t)-exp(s2*t));这里加了一项0;&nbsp, 但是没有其他影响;我取的0;&nbsp,(a/b*2+1)).t=0;x=-a.002*y1)^2+(R+2*r21)^2);&nbsp,&nbsp,&f=2;%D的步长:<a href="&nbsp.83*10^(-6);图像是精简之后的.g(n;我取的0:b;&&c2=10:c2/zhidao/pic/item/a08b87dfcb30e924b899f399;&nbsp,m)=quad(f;&%O取值区间上限.R=0;s2=0;P21=1;s1=0.^2);endn=n+1;m=0;for&nbsp.baidu.5*(-(2*Rs+P21+R)-sqrt(D));&&nbsp,&nbsp.1://e;x=-a.hiphotos://e.baidu.5://e.*y1
确实可以啦！为什么必须要加“0.*y1”这一项呢？改之前的代码到底哪里有问题呢？
那个quad函数的确有问题, 包括quadl函数也一样, 被积函数只要是常数就有这问题, 你可以试试. 不知道为什么. 加上一项0就行了.
哦，那就是这个quad函数原本就有的bug了？万分感谢大神指点！
应该是, 因为常值函数不含自变量, 可能就出错了, 我猜是这样吧. quad函数可以调试, 但是我没分析里边的代码. 反正能出来结果, 不求甚解了.
提问者评价
其他类似问题
为您推荐：
其他1条回答
因是f是你的函数的名字，导致冲突，又是你函数里变量的名字。你把你的函数改个名字就可以了
不是这里的问题哦，已经解决啦，多谢！
你问题里的错误就是我说的原因导致的，另外quad没有问题，我也没有加什么0.*，只是把你程序里的f换成了ff就可以出来图，是这样的。回你只是因为觉得你这样写程序用和函数名相同的名字取变量名会给自己造成麻烦，以后要注意。
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁请问哪位大神知道下面的MATLAB程序哪里有问题？_百度知道
请问哪位大神知道下面的MATLAB程序哪里有问题？
;R=zeros(m;[n s]=size(r2),r2(,s):):m
for j=1,r2)[m n]=size(r1),j)=max(min(r1(i,:s
R(i;for i=1function R=js(r1,j)&#39;))
提问者采纳
，才能保证min(r1(i;)里面的两个向量同维：当r1是m*n维矩阵的话:)程序本身没有错误吧 ,r2(，r2应该是n*m维的，min要求同维？运行的时候r1与r2的输入应该满足,,j)&#39
提问者评价
已经解决了，不过还是要谢谢你！
其他类似问题
为您推荐：
matlab的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁查看: 332|回复: 0|关注: 0
求助：用matlab PTB进行MOT范式的显示和操作程序问题
本帖最后由 Yorty 于
23:47 编辑
求大神帮帮忙!!很紧急!!谢谢!!
我的是32位R2009a版本的matlab,需要用到psychtoolbox工具箱,是PTB 3.0版本的.
我的实验是关于MOT范式的,实验步骤简单解释如下:
实验开始→呈现8个静止的黑色小球，目标被红色圆圈框住→被试识记目标→点击键盘空白键→红圈消失，小球开始随机匀速运动→5s后小球停止运动→被试通过鼠标左键点击要求追踪的目标，所选小球变成蓝色（被试反应时间为10s）→点鼠标右键保存数据，进入下一次实验→实验全部完成后退出屏幕，统计各种条件下的正确率（命中数/目标数）。
我求助的问题是:1.请问怎么让小球运动5s后自动静止,以便等待被试按键操作呢?我试过用tic toc、pause、WaitSecs,但是不知道是不是语句放的位置不对,还是我用错了函数,以至于运行时小球要么一直静止不动,要么运动完全不受影响..下面的截图是我新建窗口画圆的语句,中间还有一些设置让小球匀速运动、碰壁反弹、两两相碰不重叠的程序,这里就不放出来了,反正小球运动的程序是可以正常运行的.
1.png (2.95 KB, 下载次数: 0)
23:46 上传
2.png (9.3 KB, 下载次数: 0)
23:46 上传
3.png (2.01 KB, 下载次数: 0)
23:46 上传
2.请问怎么让被试在小球静止后的画面上直接进行鼠标左键按键操作,以使所选黑色小球变成蓝色呢?我想到的是用GetMouse和FillOval语句,但是当我点击鼠标左键的时候,背景就自动变成了灰色,这就使我画蓝色圆的时候不能看到目标..语句如下,是紧接着上面那个循环语句后编写的,放在ShowCursor end之后:
4.png (8.3 KB, 下载次数: 0)
23:46 上传
我有想过是因为窗口设置的是一开始openwindow的那个灰色背景w,而不是小球静止的那个画面..不知道对不对..也不知道怎么解决..
3.请问怎么判断被试点击的小球是否是目标呢?我想到的是用GetClicks语句,判断条件为检测鼠标点击时那一点是否在目标小球的圆内,若是则正确选中目标数加一.大概想法如下图所示:
5.png (6.33 KB, 下载次数: 0)
23:46 上传
但是我不知道怎么获取目标小球的圆心坐标(x0,y0),因为我的实验中目标小球不止一个..也不知道怎么写判断条件中的语句才能达到效果..
Powered by求大神指导 程序哪个地方出问题了_matlab吧_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：88,102贴子：
求大神指导 程序哪个地方出问题了收藏
I1=imread('09.bmp');[thr, sorh, keepapp]=ddencmp('den','wv',I1);I2=wdencmp('gbl',I1,'sym4',2,thr,sorh,keepapp); %小波除噪I3=medfilt2(I2,[9 9]); %中值滤波I4=imresize(I3,0.2,'bicubic'); %图像大小I5=imrotate(I4,33,'crop');BW=edge(I5,'sobel');subplot(1,1,1);imshow(BW);[H,theta,rho]=hough(BW);imshow(H,[],'XData',theta,'YData',rho,...
'InitialMagnification','fit');xlabel('\theta'),ylabel('\rho');axis on,axis normal, P=houghpeaks(H,5,'threshold',ceil(0.3*max(H(:)))); %用houghpeaks查找H中的峰值x=theta(P(:,2));
%绘制峰值y=rho(P(:,1));plot(x,y,'s','color','white'); -————————————*******——————****————lines=houghlines(BW,theta,rho,P,'FillGap',5,'MinLength',7);%查找直线figure,imshow(I5),%创建图像，将直线段叠加到原始图像上max_len=0;for k=1:length(lines)
xy=[lines(k).lines(k).point2];
plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');
% plot beginnings and ends of lines plot(xy(1,1),xy(1,2),'LineWidth',2,'Color','yellow');
plot(xy(2,1),xy(2,2),'LineWidth',2,'Color','red');
% determine the endpoints of the longest line segment
len=norm(lines(k).point1-lines(k).point2);
if (len&max_len);
endend%highlight the longest line segmentplot(xy_long(:,1),xy_long(:,2),'Linewidth',2,'Color','cyan'); 在分界线之前都没有问题，但是程序运行出来说我引用了没有定义的点 求大神指导！！！！！！
上面是我处理的图片
登录百度帐号推荐应用
为兴趣而生，贴吧更懂你。或