STM32/STM8技术论坛
直播中

硬汉Eric2013

9年用户 1035经验值
擅长:嵌入式技术
私信 关注
[资料]

【安富莱——DSP教程】第5章 Matlab 简易使用(三)

第5章  matlab 简易使用(三)
    本期教程主要是讲解Matlab的一些编程语句。
    5.1 条件控制——if,else,switch
    5.2 循环控制——for,while,continue,break
    5.3 总结

回帖(6)

硬汉Eric2013

2015-6-1 17:47:00
5.1  控制流
5.1.1  条件控制——if,else,switch
        下面我们通过三个简单的例子来说明这三个函数的使用。
l  if语句的使用
  1. % Generate a random number
  2. a = randi(100, 1);

  3. % If it is even, divide by 2
  4. if rem(a, 2) == 0
  5.     disp('a is even')
  6.     b = a/2;
  7. end
命令窗口输出结果如下:
    >>Untitled2
      a is even
l else语句的使用
  1. a = randi(100, 1);

  2. if a < 30
  3.     disp('small')
  4. elseif a < 80
  5.     disp('medium')
  6. else
  7.     disp('large')
  8. end
命令窗口输出结果如下:
           >> Untitled2
      large
l switch语句的使用
  1. [dayNum, dayString] = weekday(date, 'long', 'en_US');

  2. switch dayString
  3.    case 'Monday'
  4.       disp('Start of the work week')
  5.    case 'Tuesday'
  6.       disp('Day 2')
  7.    case 'Wednesday'
  8.       disp('Day 3')
  9.    case 'Thursday'
  10.       disp('Day 4')
  11.    case 'Friday'
  12.       disp('Last day of the work week')
  13.    otherwise
  14.       disp('Weekend!')
  15. end
命令窗口输出结果如下:
>> Untitled2
Weekend!
在这里顺便介绍一个类似于C语言中scanf的函数input并配合上面的if else实现一个小功能:
  1. yourNumber = input('Enter a number: ');

  2. if yourNumber < 0
  3.     disp('Negative')
  4. elseif yourNumber > 0
  5.     disp('Positive')
  6. else
  7.     disp('Zero')
  8. end
运行上面代码后,我们在命令窗口输入数字22,输出结果如下:
>> Untitled2
    Enter a number:22
    Positive

5.1.2  循环控制——for, while, continue, break
    这里我们也通过几个简单的例子来说明这几个函数的使用。
l for语句的使用
  1. for n = 3:32
  2.    r(n) = rank(magic(n));
  3. end
  4. r
命令窗口输出结果如下:
>> Untitled2
r =
      Columns 1 through 19
           0    0     3     3    5     5     7    3     9     7   11     3    13    9    15     3   17    11    19
      Columns 20 through 32
           3    21   13    23     3   25    15    27    3    29    17   31     3
l while语句的使用
  1. a = 0; fa = -Inf;
  2. b = 3; fb = Inf;
  3. while b-a > eps*b
  4.    x = (a+b)/2;
  5.    fx = x^3-2*x-5;
  6.    if sign(fx) == sign(fa)
  7.       a = x; fa = fx;
  8.    else
  9.       b = x; fb = fx;
  10.    end
  11. end
  12. x
命令窗口输出结果如下:
>> Untitled2
    x =
                 2.0946
l continue语句的使用
  1. fid = fopen('magic.m','r');
  2. count = 0;
  3. while ~feof(fid)
  4.     line = fgetl(fid);
  5.     if isempty(line) || strncmp(line,'%',1) || ~ischar(line)
  6.         continue
  7.     end
  8.     count = count + 1;
  9. end
  10. fprintf('%d linesn',count);
  11. fclose(fid);
命令窗口的输出结果如下:
>> Untitled2
    31 lines
l break语句的使用
  1. a = 0; fa = -Inf;
  2. b = 3; fb = Inf;
  3. while b-a > eps*b
  4.    x = (a+b)/2;
  5.    fx = x^3-2*x-5;
  6.    if fx == 0
  7.       break
  8.    elseif sign(fx) == sign(fa)
  9.       a = x; fa = fx;
  10.    else
  11.       b = x; fb = fx;
  12.    end
  13. end
  14. x
命令窗口输出结果如下:
>> Untitled2
     x =
                 2.0946

5.1.3  矢量化
    对于matlab而言,要想加快算法的执行速度可以通过算法的矢量化来实现,比如要实现如下的功能。
x = .01;
for k = 1:1001
       y(k) = log10(x);
       x = x + .01;
end
但是我们矢量化后,将更加方便和容易实现。
    x = .01:.01:10;
    y = log10(x);
但是有一点大家要特别注意,不是什么程序矢量化都能加快执行速度,要视具体情况而定。

举报

硬汉Eric2013

2015-6-1 17:48:21
5.2  Matlab中help的使用
    关于matlab入门方面的东西就跟大家将这么多,基本上有这些基础就够了,后面遇到什么问题在网上查找资料即可。也可以查看matlab本身的help帮助文档。点击这里就可以查看。
5.1.png

如果有不懂的函数,可以直接在命令窗口输入help 再加上函数即可,比如输入:
5.2.png


举报

硬汉Eric2013

2015-6-1 17:48:51
5.3  总结
    Matlab方面的教程就跟大家讲这么多,后面需要那方面知识的时候,我们再具体的补充。基本上面学会这些基本的操作就可以入门了。永远要记住,Matlab只是个工具,我们只需把它当个工具来用,没有必要花大量的时间去系统的学习,入门后用什么学什么即可。

举报

武力戡乱

2015-6-2 00:02:40
你家出DSP了????
举报

硬汉Eric2013

2015-6-2 10:23:33
引用: 武力戡乱 发表于 2015-6-2 00:02
你家出DSP了????

没有,这个是ARM官方的DSP库
举报

武力戡乱

2015-6-2 13:56:47
引用: 硬汉Eric2013 发表于 2015-6-2 10:23
没有,这个是ARM官方的DSP库

我看好你家啊!
举报

更多回帖

发帖
×
20
完善资料,
赚取积分